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CORRECTION 


Canadian Journal of Physics, Vol. 36. On p. 803 the following illustration should appear 
above (c) and (d) in Fig. 2. 
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LOW TEMPERATURE BEHAVIOR OF 
DEBYE CHARACTERISTIC TEMPERATURES! 


G. K. Horton? AND H. SCHIFF 


ABSTRACT 

The low temperature behavior of the characteristic Debye temperature is in- 
vestigated using the six- and nine-direction Houston approximation and by way 
of illustration, the lattice dynamics of Born and Begbie (nearest neighbors only). 
We study the values of 69 and Bz for Cu, Au, Ag, Al, and Pb in the formula 
Or = 8 [1—B.(T/8»)*. . .]. The existence in Pb of a shallow maximum before 
the usual minimum, found earlier by Bhatia and Horton, is confirmed. The 
method is quite general and may be used for non-metals and with other theo- 
retical models. 


The low temperature behavior of the characteristic Debye temperature 
was worked out in detail by Bhatia and Horton (1955).* Their method 
determined the values of the a,’s in the expansion of the frequency distri- 
bution function (Blackman 1955), 


(1) G(w) = VY aw™ 


The equivalent Debye temperature 07 at a temperature T is then given by 
T 2n 
(2) Or = Oo», he ’ (T K ®o) 
n 0 


2 
(3) where Bo=1, B= — (20*/21) (ax/a)(# 0.) , etc. 


9No\ hi 
(4) and Qo = (2s) zs 


Qa) 


The values of the a,’s may be obtained by using Houston’s approximation 
(Houston 1948). This was done by Bhatia and Horton, who used Bhatia’s 
lattice dynamical model (Bhatia 1955) to find the value of B, for a number of 
face-centered and one-body-centered cubic metal. They showed that the 


1Manuscript received May 14, 1958. 

Contribution from the University of Alberta, Edmonton, Alberta, Canada. 

2Now at Solid State Section, Division of Pure Physics, National Research Council, Ottawa, 
Canada. 


*Similar results may be found in Salter (1956). 
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steepness of the drop of the 07 vs. T curve for a given crystal could be under- 
stood in this way, in fact that for some crystals the drop was much steeper 
than for others; indeed it was found that for Pb, using room-temperature 
values of the elastic constants since no 0° K values were available, the O7 
vs. T curve went through a shallow maximum before the usual minimum. 

The effect of using a different lattice dynamical model, again using Houston’s 
method, was investigated by Horton and Schiff (1956). They showed how the 
substitution of either the conventional two-neighbor central force model 
(Leighton 1948) or the nearest-neighbor non-central force model (Begbie and 
Born 1947) in place of the Bhatia model altered the conclusions of Bhatia 
and Horton; qualitatively they were unchanged. It seems certain that none 
of these models provides a particularly good description of any actual metal, 
but the broad similarity of many of the results obtained using different 
models suggests strongly that a further study of any one of them may throw 
light on the low-temperature behavior of crystals. 

In this paper we report the effect on these earlier low-temperature results 
of using the various refinements of Houston’s method that have been studied 
recently by Betts, Bhatia, and Wyman (1956) and by Betts (1958). We have 
chosen the Born and Begbie model for the present calculations, both because 
of its simplicity and because it appears more appropriate for the low-tempera- 
ture description of metals than the equally simple central force theory (Horton 
and Schiff 1958). The results of our calculation are summarized in Table I. 


TABLE I 


THE VALUES OF THE PARAMETERS @) AND By OCCURRING IN FORMULA (2) FOR THE LOW- 
TEMPERATURE DEBYE Or 











Metal —o = 2—(¢::—€12)/C44 Oo e, Oy» B,) B® B® 
Al 0.311* 427.71 428.14 428.13 17.35 15.77 15.82 
Au 1.307 155.43 162.77 162.19 6.327 4.480 5.155 
Ag 1.33¢ 218.50 228 .37 227 .62 5.545 9.347 9.445 
Cu Lote 329.72 345.64 344.43 1.950 8.817 8.684 
Pb 1. 50§ 82.98 89.98 89.38 —-19.03  —3.647 —3.925 








Note: The 0° K elastic constants were used in the following three references. 

*Sutton (1953). 

TNeighbours and Alers (1958). We would like to thank Dr. Alers for communicating his results to us before 
publication. 

tGaffney and Overton (1954). 

§The results for Pb were obtained using room-temperature values of the elastic constants. In order to make a 
rough estimate of the effect of 0° K constants we proceed as follows. Let us assume a linear temperature dependence, 
the same for all the c's. Further let the slope be chosen so as to yield the correct experimental 90 = 94.59 K (Keesom 
and Pearlman 1956), using the six-term Houston formula, which is known to give a result accurate to better than 
1%. It is then easy to show that the values of the B's are unchanged. 


The parameter o is a measure of the anisotropy of the metal, while 09, 
6o°®, Oo, and the corresponding B.’s were obtained using the Houston three-, 
six-, and nine-direction approximations respectively, which may be found 
in the papers referred to earlier. The characteristic temperatures Oo agree 
to 0.1% with those evaluated by the method of Budzinski and Schiff (1957). 
The correctness of the B2’s may be gauged from the convergence exhibited 
in Table I. (The B.y’s are made up of positive and negative contributions, 
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in contrast to 69, which consists of positive contributions only. Consequently 
their convergence is correspondingly slower.) The @o’s could also be obtained 
from the numerical tables compiled by de Launay (1956). These tables 
cannot be used to look up the B,’s, though tables for the B2’s have been 
announced by Marcus and Kennedy (1958) for certain lattice dynamical 
models. The low-temperature behavior for the central force model has also 
been discussed by Leighton (1948). 

The accuracy of the present approximations would seem to be sufficient 
to establish beyond reasonable doubt that the value of B, for Pb using this 
model and room-temperature elastic constants is indeed negative although 
it is now so small that it will be difficult to distinguish it experimentally from 
zero. A much larger positive value of Bz, has been reported by Dolecek (1955), 
but it does not seem to us that the experimental accuracy justifies a unique 
interpretation of his results. In particular Dolecek’s values for y and @o 
are higher than those estimated from other measurements (Keesom and Pearl- 
man 1956); it is perhaps significant that with a lower value of y, Dolecek’s 
results would imply a lower value of both @9 and Bz. 

When using the three-direction Houston formula, different models were 
found to give similar results. This may be modified when a more realistic 
lattice dynamical model becomes available. One might hope, however, that 
the present models may not be too different from a correct model and hence 
that the results even for the Born and Begbie nearest-neighbor model may 
have some physical significance. 

We would like to stress that our procedure is available for any model and 
may be applied to non-metals. The convergence seems sufficiently good to 
justify it for the analysis of the experimental results and to find the higher 
B’s. In many cases the six-direction formula will give good results when the 
crystal is reasonably isotropic. For non-cubic crystals, to which our procedure 
is equally applicable, one can replace the Houston method by those developed 
by Betts, Bhatia, and Horton (1956). 
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THE FREEZING POINTS OF HIGH PURITY METALS 
AS PRECISION TEMPERATURE STANDARDS 


IV. INDIUM: THERMAL ANALYSES ON THREE GRADES OF CADMIUM! 


E. H. McLAREN 


ABSTRACT 


An investigation of the freezing and melting temperatures of a sample of high 
purity indium (Cominco 99.999+ %) has been made. Plateaux of essentially 
constant (<=+ 0.0001° C) temperature with durations for over 1 hour are readily 
obtained on the cooling curves of induced freezes on this metal. The standard 
deviation of the plateau temperature (liquidus point) from a series of 24 induced 
freezes was +0.0001° C. The pressure effect on the freezing temperature of indium 
was found to be 0.0049° C for 1 atm. Alloy melting ranges were measured fol- 
lowing different types of freezing. 

An extensive intercomparison of liquidus points and alloy melting ranges has 
been made on a sample of 99.99°% cadmium and two samples of 99.999+ % 
cadmium from different suppliers. The liquidus points of the high purity samples 
were indistinguishable using precision resistance thermometry but one sample 
melted over a slightly smaller range of temperature than the other. Both these 
samples showed minor arrests on melting curves after induced freezing and 
detailed analyses of the melting contours after various types of freezing indicated 
some evidence of characteristic structure inside a range of 0.002° C. 


INTRODUCTION 
A study of the liquid—solid transformation temperatures in a sample of 
high purity indium and comparative thermal analyses on three grades of 
high purity cadmium have bee: made. This work is part of a continuing 
investigation of the liquidus temperatures and alloy melting ranges of high 
purity metals to determine the suitability of metal freezing points as precision 
temperature standards and the usefulness of precision resistance thermometry 
to discriminate among high purity metals. 
EXPERIMENTAL 
Indium 
The high purity indium was obtained from the Consolidated Mining and 
Smelting Company, Trail, B.C., Canada, and Table I shows a typical analysis 
of this material and an estimate of the expected depression of the liquidus 
point from that of the pure element. A 1.2-kg sample was used under a dry 
nitrogen atmosphere. 
(a) Thermal Analysis 
Following the techniques developed earlier to investigate the freezing and 
melting temperatures in high purity zinc, cadmium, and tin (McLaren 1957a, 
b, 1958), a similar series of thermal analyses was made on this indium sample. 


1Manuscript received May 26, 1958. ; 

Contribution from the Division of Applied Physics, National Research Council, Ottawa, 
Canada. 

Issued as N.R.C. No. 4817. 
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TABLE I 
TYPICAL IMPURITY CONTENT OF COMINCO* HIGH PURITY INDIUM 


Estimated liquidus depression 


Impurity content, 





Impurity p.p.m. (1G) Phase diagram 
Zn N.d.f — 
Cd 1 0.0001 Valentiner (1943) 
Pb 2 0.0001 Valentiner and Haberstroh 
(1938) 
Bi N.d. — 
Ag <1 <0.0005 Hellner (1951) 
Cu 1 0.0003 Weibke and Eggens (1934) 


Total <0.0010 


*See Turner, Mitchell, and Bryan (1957). 
+N.d.—not detected by routine spectrochemical analysis. 


In this paper normal freezing refers to undisturbed natural solidification 
while induced freezing relates to a technique that is used to shorten the 
recalescent period. The latter is effected by freezing a thin mantle of metal on 
the axial thermometer-wel!l after nucleation has commenced; this releases 
sufficient latent heat to raise the melt temperature to the liquidus temperature. 

Figure 1 shows freezing curves that were obtained using three different 
freezing techniques, curves I, II, and III arising from very fast normal, very 


MINUTES 





awn 3 ic HOURS T 
e Se 
ern. me °C 
I 
ee en T 
a" 


“Oy 
we NR ~0.01°C 


$ 





30 60 90 120 I 
MINUTES 


Fic. 1. Typical freezing curves on high purity indium: I, very fast normal; II, very slow 
normal; III, slow induced. 
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slow normal, and slow induced freezing respectively. The time scales on these 
plots begin at the end of the supercool which for the indium sample in the 
experimental arrangement used here amounted to about 0.07° C. These freezing 
curves are similar to those obtained with high purity zinc and cadmium; 
the plateaux developed on induced freezing vary less than 0.0001° C for periods 
depending on the gross rate of cooling of the furnace block; nucleation and 
recalescence, however, are measurably slower than in zinc and cadmium. 
Because of the slower recalescence in indium, temperature measurements may 
be made during the rise from supercool even on very fast freezes (see curve I, 
Fig. 1) and the sections of essentially constant temperature on the thermal 
curves of very slow normal freezes (curve II, Fig. 1) have very much shorter 
durations than are obtained on similar freezes of zinc and cadmium. To assist 
the melt to the liquidus point quickly on induced freezing it was found ex- 
pedient to nucleate the freeze at a more rapid cooling rate than was required 
subsequently to obtain a slow induced freeze. 

The tendency to a retarded recalescence in indium is most likely due to its 
smaller latent heat of fusion compared with those of zinc and cadmium. 

The plateaux of nearly constant temperature that are obtained on slow 
induced freezes on this high purity indium provide an excellent index for a 
temperature standard and, as do the plateaux on slow induced freezes of 
zinc and cadmium, probably arise from extended liquidus arrests caused by 
coring as the sample solidifies. 

Prior to determining the reproducibility of the plateau freezing temperature 
of, this sample, we made extensive measurements of the alloy melting ranges 
following various types of freezing to gauge the thermal effects of the dis- 
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Fic. 2. Typical melting curves on high purity indium after: I, very fast normal freeze; 
II, slow induced freeze; III, very slow normal freeze. 
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solved impurities. Figure 2 shows the thermal curves obtained on remelting 
the ingot immediately after the corresponding freezes of Fig. 1. The melting 
curves after the very fast and the slow induced freezes (curves I and II, Fig. 2) 
show melting occurring over a few thousandths of a degree only, but after 
slow normal freezing significant melting takes place over a range of about 
0.01° C (curve III, Fig. 2). The minor arrest on the remelting curve after 
induced freezing is probably due to a sharp segregation of impurities during 
the freeze; this segregation would account for the extremely flat major arrest 
on this melting curve. In contrast, the melting curve after the very slow normal 
freeze reveals that the impurities were received by the new solid during the 
slow recalescence over a broad range of temperatures and the small step near 
the liquidus break is the only evidence of heavy coring in the ingot. 

Figure 3 shows derived inverse plots of melting curves obtained after 
five different types of freezing; Table II outlines the transformation conditions 
that existed for the group. In order to obtain a large number of points along 
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Fic. 3. Inverse melting curves on high purity indium after freezes of type (a) very fast 
normal, (6) fast normal, (c) slow induced, (d) slow normal, and (e) very slow normal. 
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each melting curve a shortened measurement technique was used with NR 
resistance balances at 2-ma current, single daily measurements at the triple 


TABLE II 
TRANSFORMATION CONDITIONS FOR THE MELTING CURVES OF FIG. 3 


duration 


Freezing Solid Melting* Melting liquid to 
Type of Freezing duration anneal duration curve liquid 
freeze curve (minutes) (minutes) (minutes) (Fig. 3) (minutes) 
Very fast Fig. 1,1 7 66 147 (a) 220 
Fast normal _- 55 95 155 (bd) 305 
Slow induced Fig. 1, III 130 107 148 (c) 385 
Slow normal - 130 111 154 (d) 395 
Very slow 
normal — 310 45 145 (e) 500 





*Total melting times of Fig. 3. 


point of water, and approximate corrections for pressure over the melt.f For 
these reasons the difference indicated in Fig. 3 between the plateau tempera- 
ture (liquidus point) obtained on induced freezes and the melting tempera- 
tures may be in error by as much as 0.0005° C. On the other hand relative 
errors along a particular curve due to slowly changing gas pressures over the 
melt would be less than 0.0001° C. 

The histograms in Fig. 3 show a strong dependence on the thermal history 
of the sample. The width at the 1% levelf for the main arrests varies from 
about 0.0005° C after induced freezing (curve (c)) to over 0.004° C after very 
slow normal freezing (curve (e)). The wider alloy melting range after the very 
slow unassisted freeze shows that significant amounts of solid are laid down 
over a wider spread in impurity concentration that was due to the slow re- 
calescence and also probably due to some impurity diffusion in the solid 
during the long time interval in which the sample is held near its liquidus 
temperature. 

These differences in melting ranges following different rates and types of 
freezing emphasize the importance of using similar thermal treatments when 
comparing samples on the basis of their alloy melting ranges. 

(b) Pressure Effect on the Freezing Temperature 

The effects of external pressure on the plateau freezing temperature of 
indium were experimentally determined and the results are shown in Fig. 4. 
The increase in the freezing temperature of indium with l-atm gas pressure 
over the melt was 0.0049° C. An immersion of about 8 cm (of the mid-point 
of the sensing element of a Meyers thermometer) was required before the 
measured temperatures agreed within 0.0001° C of the variation deduced 


TNR balances mean the average of the resistance balances obtained with the thermometer 
connected in the Normal and Reversed positions on the Mueller-Wheatstone bridge. 

The width at the 1% level is an arbitrary melting range parameter and is simply the melting 
width at the 1% level on the histograms that are derived by determining the fraction of the 
total melting time spent in given temperature intervals from the direct time-temperature 
curves (McLaren 1958). 
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from the Clausius—Clapeyron relation. This latter variation is shown on the 
solid line on curve I and was calculated using the above measured gross 
pressure effect and a liquid density at the freezing point of 7.03 g/cc as de- 
termined by Williams and Miller (1950). 

(c) Reproducibility of Plateau Freezing Temperature 

Figure 5 shows results that were obtained on the reproducibility of the 
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Fic. 4. Pressure effects on the plateau freezing temperature of indium. 
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Fic. 5. Reproducibility of the plateau freezing temperature of high purity indium. 
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plateau freezing temperature (liquidus point) of this indium sample taken from 
a series of slow induced freezes carried out during a 3-week interval. Ry, is the 
resistance of the thermometer at the plateau freezing temperature corrected 
to l-atm gas pressure over the melt. The mid-point of the sensing element 
of the Meyers thermometer was immersed about 15 cm in the melt. Rrp is 
the average of the thermometer resistance at the triple point of water before 
and after the appropriate indium freeze. Both R;, and Rrp are zero current 
resistances corrected to unity bridge ratio and true bridge zero expressed in 
bridge units of resistance (McLaren 1957a). 

A single water triple point cell was used throughout the investigation; one 
ice mantle covered the work on freezes 1-8 inclusive and another mantle was 
used for the remainder of the work. Prior to these measurements the ther- 
mometer was annealed at about 450° C to remove as much as possible of the 
added resistance due to accidental cold-working of the platinum, and ad- 
ditional overnight anneals were made after freeze numbers 8, 12, 15, and 19 
respectively in order to maintain the thermometer in an annealed state. 

The standard deviation of the ratio Ry,/Rrp for the 24 induced freezes on 
this high purity indium was equivalent to +0.00009,° C. 


Cadmium 

Earlier work (McLaren 19570) had shown that the liquidus point (321.032°C) 
of a sample of New Jersey Super Pure cadmium had a reproducibility and 
stability comparable with that of any precision-fixed point on the International 
Temperature Scale. Two additional samples of high purity cadmium have 
been provided by the Consolidated Mining and Smelting Company, one a 
sample of special high purity and the other a sample of normal high grade 
Tadanac cadmium. Extensive thermal analyses have been made on the three 
samples. This work soon resolved into an attempt to distinguish between the 
two samples of highest purity on the basis of either alloy-melting ranges 
following various freezing techniques or liquidus temperatures measured on 
induced freezes. Table III lists the samples and gives the available information 
on their impurity contents as supplied by the sources. Estimated depressions 
of the liquidus points from that of the pure element are <0.0013° C and 


TABLE III 
IMPURITY CONTENTS OF CADMIUM SAMPLES 


Impurity content* (p.p.m.) 











Sample Source and grade Analysis Pb Cu Tl Ag Zn Fe Mg Al 

Cd.3 New Jersey Zinc Company, Typical N.d.V.f. — — N.d. N.d. V.f. V.f. 
Super Pure 

Cd.4 Consolidating Mining and Cd.4 1 <1 <I1INd -—- -—- —- — 


Smelting Company, High 
purity No. 9521 
Cd.5 Consolidated Mining and Typical 60 20 200 — >530 — —- — 
Smelting Company, 
Tadanac 


*N.d.—not detectable in qualitative spectrochemical analysis. V.f.—'‘may be very faintly visible in the spectro- 
gram but may represent electrode or atmospheric contamination rather than impurities in the metal’’; E. A. 
Anderson, New Jersey Zinc Company, Palmerton, Pa., U.S.A. 
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<0.046° C respectively for the impurity contents listed for samples Cd.4 and 
Cd.5. Sample Cd.3 was loaded in September, 1955, and has been used fre- 
quently for fixed point work since that time, while sample Cd.4 was loaded in 
November, 1957, especially for this comparison. 

(a) Thermal Analyses 

The three samples were allowed to transform under a variety of freezing 
conditions in order to make comparisons of the freezing curves and subsequent 
remelting curves. 

Figure 6 shows the results of very fast freezing on the three samples; these 
freezes were made by extracting the pyrex tubes holding the crucibles of cad- 
mium out of the furnaces and allowing the transformations to take place in the 
air. Although sample Cd.5 may be easily distinguished from the purer samples, 
it is immediately apparent that samples Cd.3 and Cd.4 have nearly the same 
impurity contents. 
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Fic. 6. Very fast normal freezes on three cadmium samples. 


Figure 7 shows freezing curves obtained from fast normal freezes on all 
samples; the freezes were made by allowing the samples to transform freely 
with the furnace power supplies turned off. The thermometer indicated a 
supercooling of about 0.03° C in samples Cd.3 and Cd.4 and less than 0.01° C 
in sample Cd.5. Although these freezes and additional similar freezes showed 
distinctive differences in the shapes of the curves during recalescence for Cd.3 
and Cd.4, it would be difficult to attribute this entirely to differences in 
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impurity contents of the two samples since the initial rate of nucleation of the 
solid phase is probably affected by the condition of the walls of the graphite 
crucibles and the degree of supercooling. This crucible factor might not be 
identical for both samples because sample Cd.3 had been used frequently 
during a 2-year interval prior to these measurements. 








f _ fel ge Pe ONO Oral OC apatag 4 
Py oe © =o—y 
~0.01°C Pr 00, 
VT f x Noon 
cules z PS 4Cd-3 oCd-4 * 





ge ee 
10 20 30 


° 
MINUTES 


Fic. 7. Fast normal freezes on three cadmium samples. 


For these reasons it has been more profitable to investigate melting ranges 
after various types of freezing to discriminate among samples of high purity 
metals, and Fig. 8 shows melting curves that were obtained immediately 
after the corresponding freezes of Fig. 7. Samples Cd.3 and Cd.4 melt over a 
temperature range of a few thousandths of a degree and there is evidence of 
structures on both their melting curves: melting can be detected over a tem- 
perature range greater than one degree for sample Cd.5. 

Figure 9 shows induced freezes on the three samples. The plateau freezing 
temperatures for samples Cd.3 and Cd.4 are indistinguishable on these freezes 
while a plateau of nearly constant temperature was not obtained for sample 
Cd.5. This failure to develop a plateau is probably a consequence of the larger 
impurity content of sample Cd.5; the slope of the liquidus for this alloy, 
over the range of compositions initially frozen out, must be large enough for 
freezing to take place over a measurable range of temperature even when the 
sample is assisted to the liquidus point by the induced freezing technique. 

Figure 10 shows derived inverse plots of the melting curves obtained after 
these induced freezes. It now becomes possible to make some distinction 
between sample Cd.3 and Cd.4. Although the melting ranges at the 1% level 
of the total melting time for the high purity cores of both samples are about 
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0.003° C it is clearly evident on comparing the histograms that substantial 
amounts of melting takes place over a larger temperature range in Cd.3 than 
in Cd.4. Both samples show minor arrests lying close to the main arrests. 
For completeness the inverse plot of the melting curve after the induced freeze 
on Cd.5 is also shown in Fig. 10. 
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Fic. 10. Inverse melting curves after the slow induced freezes of Fig. 9. 


Since the foregoing work showed that the thermal curves for Cd.3 and Cd.4 
were nearly identical under similar thermal treatments, further efforts were 
made to try to reveal distinctive differences in these samples by carefully 
examining melting contours after a variety of freezing conditions. Table IV 
summarizes transformation conditions and melting range estimates made 
on the widths at the 1% level on inverse plots of several series of direct melting 
curves; estimates of the melting ranges of the least pure sample, Cd.5, are also 
included in this table. Series 1 to 5 are similar to types of transformations 
already discussed but series 6 and 7 require further explanation. Curves II 
of Fig. 11 show composite inverse plots of the five melting curves in series 1 
to 5 for each sample which were fitted at the temperature intervals in which the 
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TABLE IV 


TRANSFORMATION CONDITIONS AND ALLOY-MELTING RANGES FOR THE THREE CADMIUM SAMPLES 
ESTIMATED FROM INVERSE PLOTS OF THE DIRECT MELTING CURVES DERIVED FROM THE FRAC- 
TION OF THE TOTAL MELTING TIME SPENT IN TEMPERATURE INTERVALS I, II, III, anp IV 


Melting range at the 1% level of the total melting 
time 
Prior freeze I II III IV 
—_—_—_—_—___—_————-__ Melting —- oo 
Duration duration Cd.3 Cd.4 Cd.3 Cd.4 Cd.5 Cd.5 








Series Type (minutes) (minutes) (°C) Ce) C7) (CE). AG) OE) 
l N cenit) 5 80 0.0038 0.0033 0.0028 0.0025 
2 2NR 5 120 0.0032 0.0024 0.0024 0.0020 0.20 0.13 
3 (Normat 40 65 0.0041 0.0034 0.0029 0.0026 
4 2NR 40 130 0.0044 0.0033 0.0029 0.0026 0.18 0.15 
5 (inane 160 90 0.0046* 0.0043* 0.0042* 0.0038* 
2NR 160 90 0.0032¢ 0.0033¢ 0.0028f 0.00247 0.30 0.18 
6 Composite 
( series 
1-5 —- - ~- 0.0028 0.0025 — — 
7 ae) 
2N 5 120 0.0038t 0.0029 0.0021§ 0.0021 — — 





*Complete samples. 

tHigh purity cores. 
{Temperature interval 0.00022° C. 
§Temperature interval 0.00011° C. 


greatest amount of melting occurred; the melting ranges from these average 
plots are listed as series 6 in Table IV. Curves I of Fig. 11 (series 7) were 
determined using 2N bridge balances only; resistance was estimated to the 
nearest micro-ohm (~0.000011° C).* 

The results in Table IV re-emphasize the importance of using similar 
thermal treatments and temperature analyses when comparing the alloy 
melting ranges of different samples with the resistance thermometer. 

Inspection of Fig. 11 shows that although the melting contours of the two 
samples are very similar sample Cd.4 appears to melt over a slightly narrower 
range of temperature than sample Cd.3. The close correlation of curve I 
to curve II fer a particular sample indicates that the 2N technique may be 
successfully used to compare melting contours of samples with narrow melting 
ranges provided the resistance thermometer is used in an experimental arrange- 
ment favoring a nearly constant temperature distribution over the connecting 
lead resistances. Table V gives a comparison of the incremented percentages 
of the total melting time spent in given temperature ranges from the plots of 


*The 2N technique allowed two measurements to be made each minute compared to about 
one measurement in 2 minutes using 2NR techniques. The 2N technique is open to criticism 
because any significant slow change in the lead resistance to the thermometer would tend to 
smear any structure on the melting contour; however, in practice it has been frequently 
observed that a Meyers standard thermometer immersed in a zinc ingot on a very slow induced 
freeze will maintain 2N balances (2R as well) to 10 micro-ohms for periods of at least 2 hours. 
Stability of this order in the lead resistances of a Meyers thermometer is undoubtedly due to 
the thermometer’s being deep-bedded in the heavily insulated furnace and to the use of heavy 
conductors connecting the thermometer leads to the Mueller bridge, all apparatus being 
maintained in an air-conditioned laboratory. 
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Fic. 11. Inverse melting curves: I, after very fast normal freezes, and II, fitted composite 
of five melting curves for samples Cd.3 and Cd.4. 


Fig. 11; the sensibly higher fraction of sample Cd.4 that melted in a range of 
0.001° C probably arises from a slightly different impurity content than is 
present in sample Cd.3. 

(b) Comparison of Liquidus Points 

A careful intercomparison of the plateau temperatures developed on slow 
induced freezes of the cadmium samples was carried out using two furnaces 
and the measuring sequence that was used earlier (McLaren 1958) to compare 
the liquidus points of samples of high purity zinc. The resistances at the 
cadmium plateau temperatures were corrected to l-atm nitrogen pressure over 
the melts. Measurements were made with the mid-point of the sensing coil of a 
single Meyers thermometer immersed about 15 cm in the cadmium. One ice 
mantle on a single water triple point cell was used to baseline all measurements. 
Fig. 12 gives results of comparing the plateau temperatures of samples Cd.3 
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TABLE V 


PERCENTAGE OF THE TOTAL TIME SPENT IN GIVEN MELTING RANGES DERIVED FROM THE 
HISTOGRAMS I AND II oF Fic. 11 





I (2N) II (2NR) 
Melting (a) (b) Melting (a) (b) 
range Cd.4 Cd.3 range Cd.4 Cd.3 
(10°*°C) %) (%) a/b (10-* °C) (%) (%) a/b 
24 17 12 1.4 1.4 17 10 1 Oy 
2.2 25 19 1.3 2.8 27 19 1.4 
3.3 32 26 1.2 4.2 35 26 1.3 
4.4 38 33 1:2 5.6 43 34 1.3 
5.4 47 40 1.2 6.9 54 39 1.4 
6.5 53 43 1.2 8.3 59 45 1.3 
7.6 62 48 1.3 9.7 63 53 | 
8.7 66 53 1.2 yt 66 58 et 
9.8 72 58 1.2 12.5 70 63 1.2 
10.9 74 65 1.1 13.9 73 66 EE 
12.0 79 69 1.1 15.3 76 68 ‘2 
13.1 81 70 Vt 16.7 79 72 j | 
14.1 82 72 tI 18.1 81 74 Ek 
15.2 83 75 Lok 19.4 83 76 Rok 
16.3 84 76 Lea 20.8 85 78 ri 
17.4 86 78 La 22.2 86 80 t) 
18.5 87 79 1-4 23.6 87 81 bo 
19.6 88 81 ed 25.0 88 82 a 
20.7 89 83 | 26.4 89 83 4 
| OS 89 84 Lo4 27.8 90 84 Ra 
22.8 90 85 14 
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Fic. 12. Intercomparison of plateau temperatures of samples Cd.3 and Cd.4. 


and Cd.4 with five determinations made on each sample. The intercom- 
parisons are summarized in Table VI. 
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TABLE VI 
INTERCOMPARISON OF CADMIUM PLATEAU TEMPERATURES 
; Reproducibility 
Ratio Rca/Rrp (equivalent 10-*°C) Plateau temp. 
No. of }=-—————_ _ ———__——_. rel. to Cd.3 
Sample freezes Mean 108 X Min. dev. Min. dev. Spread ce) 
Cd.3 5 2.21848667 +43 +1.2 2.9 
Cd.4 5 8627 +42 +1.2 3.7 —0.00011 
Cd.3, Cd.4 10 8647 +46 +1.3 4.0 
Cd.3 1 2.21848697 


Cd.5 1 0535 —0.023 








DISCUSSION 


The high reproducibility of the liquidus point of high purity indium is 
comparable with that of any fixed point on the International Temperature 
Scale and experimentally this freezing point is incomparably easier to realize 
than is the attainment of a precision Steam Point. For this reason the indium 
point (156.611° C Int. 1948) should receive serious consideration as a precision 
alternative to the Boiling Point of Water on the International Temperature 
Scale. 

The extreme difficulty encountered in attempting to distinguish between 
cadmium samples Cd.3 and Cd.4, even though one sample (Cd.3) had been 
used for over 2 years as a laboratory secondary standard prior to these measure- 
ments, illustrates the high stability and ready attainment of the liquidus 
point of high purity cadmium. 

Two distinct advantages of the liquidus points of high purity metals as 
precision temperature fixed points are: corrections can be easily and accurately 
made for pressure effects on the freezing temperatures, and alloy melting ranges 
can be occasionally measured to guard against decreases in purity and associ- 
ated changes in the liquidus points. For the 6 years during which studies 
have been made in this laboratory on the transformation temperatures of 
samples of high purity Zn, Cd, Sn, and In, no evidence of contamination of 
any melt has been observed; in fact the liquidus points of these metals, with 
the triple point of water, have provided a means for examining under high 
precision the stability of the temperature coefficients of the resistance thermo- 
meters. In addition the liquidus points of high purity metals are eminently 
suitable for use in comparing the temperature scales realized in different 
laboratories: it is no longer necessary to base such comparisons on transported 
resistance thermometers whose temperature coefficients are known to change 
as a result of cold-working during transit. 

Further work is contemplated on the freezing temperatures of high purity 
tin, lead, and bismuth. Earlier investigations on samples of 99.99% tin from 
two different sources had revealed highly reproducible plateau temperatures 
on particular samples (McLaren 19576) but a difference of about 0.008° C 
between the liquidus points of the samples; 99.999+ % tin has now been 
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procured and thermal analyses and annealing studies will be made on this 
metal. The highly reactive alkali metals with thallium and gallium are the 
remaining metals that have freezing temperatures in the range 0° C to 420° C; 
and although the freezing points of the former do not appear particularly 
promising as temperature standards, precision thermal analyses on high purity 
samples of these metals may provide useful metallurgical information. It is 
also hoped to examine zone-refined samples of metals since this might reveal 
the limit at which impurity differences can no longer be detected by precision 
thermal analysis. 

The structures that are developed on the melting contours of these high 
purity metals fall into two types. The minor arrests seen in Fig. 3 (c) and 
Fig. 10 (Cd.3, Cd.4) are almost certainly due to the segregation of impurities 
that are nearly insoluble in the solid solvent. The structures shown inside the 
melting contours of the high purity cores in Fig. 11 may be partly due to more 
soluble impurities, but what the “‘natural melting range’’ of the pure element 
would be and whether or not any isotope fractionation occurs for specific 
freezing conditions that could be detected on the melting contour is not known. 
Theoretical considerations by MacDonald and Lovejoy (1955) on the effects 
of differences in zero point energies on the melting points of the zinc isotopes 
suggested that significant fractionation on freezing for the heavy elements 
would be unlikely. However it would be desirable to have mass analyses made 
on sectioned ingots following different types of freezes to verify these theoretical 
expectations especially in view of the well-established fractionation of D 
and O!8 in normal water on freezing (Teis and Florensky 1941) and the very 
interesting results reported by Alieva (1958) on differences in isotopic com- 
position among samples of high purity zinc. 


CONCLUSIONS 


This investigation has shown that it is possible to develop sections of very 
steady temperature that are reproducible to the order of +0.0001°C on 
induced freezes of high purity indium. The forms of the thermal curves 
obtained on high purity indium are similar to those obtained with high purity 
zinc and cadmium. 

No measurable difference was found between the liquidus points of the 
two samples of cadmium with highest purity but close inspection of melting 
contours showed that one sample melted over a slightly smaller range of 
temperatures than the other. 
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CONDENSATION COEFFICIENTS OF SILVER, GOLD, AND 
COPPER IN VACUUM DEPOSITION! 


S. CHANDRA AND G. D. Scott 


ABSTRACT 


A beam of atoms striking a substrate will in general not condense completely, 
but some of the incident atoms will re-evaporate. The condensation coefficient— 
defined as the ratio of the number of atoms condensing to the total number in- 
cident on the substrate—has been determined for silver at substrate temperatures 
from — 30°C to 170° C, and for gold and copper at temperatures from 25° C to 
165° C. For silver at room temperature the condensation coefficient is found to 
be 0.936, which indicates 6.4% of the incident atoms are re-evaporated or 
diffusely reflected. By a simple graphical analysis of the results it is shown that 
the temperature below which complete condensation occurs is 200° K for silver, 
230° K for gold, and 180° K for copper. 


INTRODUCTION 


In 1915 Wood showed that mercury atoms condense completely onto a 
glass substrate only at temperatures below —140° C. Similarly, for zinc, 
mercury, and cadmium, Knudsen (1916) found that the temperatures below 
which complete condensation of the atoms takes place are between —183° C 
and —78° C. These authors also showed that in each case, above a certain 
substrate temperature, there is no condensation of the metal. In their experi- 
ments, the stream density, i.e., the number of atoms striking a unit area of 
a substrate normal to the forward direction in 1 second, was constant, and 
they did not attempt to find the effect of the stream density on the conden- 
sation temperature. However, Cockcroft (1928) has shown that the tempera- 
ture above which cadmium atoms fail to condense on the substrate of polished 
copper or silver disks depends on the stream density of the atoms. 

Frenkel (1923) has given a theoretical discussion related to these phenomena. 
For the condensation of metal atoms on a substrate, he assumed (a) that the 
metal atoms after striking the surface move over it as atoms of a two-dimen- 
sional gas for a short time, called the mean free time, and then re-evaporate, 
and (b) that these atoms while moving on the surface may collide with one 
another to form pairs which have a greater mean free time than a single atom. 
The pairs, or “doublet atoms,” act as nuclei of condensation for the atoms of 
the two-dimensional gas. Frenkel showed that if fo is the period of oscillation 
of an atom normal to the surface of a substrate about a point where its 
potential energy is #0, then its mean free time ¢ is given by 


t = to exp (uo/RT), 


where & is the Boltzmann constant and T is the absolute temperature of the 
substrate. The above expression shows that, as T increases, ¢ decreases, i.e., 
the probability of re-evaporation of an atom is greater from a hot substrate 


‘Manuscript received June 9, 1958. 7 
Contribution from the Department of Physics, University of Toronto, Toronto, Ontario. 
Can. J. Phys. Vol. 36 (1958) 
1148 
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than at lower temperatures. The condensation coefficient of a metal defined 
as the ratio of the number of atoms condensing on a substrate to the total 
number striking it will therefore decrease as the temperature of the substrate 
increases. 

Devienne (1953) made a study of the effect of the substrate temperature 
on the condensation coefficients of antimony, silver, and gold by the use of 
their radioactive isotopes. His tables for antimony show a decrease in the 
condensation coefficient of the metal with an increase in substrate temperature. 
For silver and gold Devienne was not able to obtain ‘‘very precise results”’ 
because of the low activity of the isotopes. However, repeated experiments 
on substrates of glass and copper at room temperature indicated that the 
condensation coefficient is not precisely equal to unity but fairly close to this 
value, ‘‘ranging between 0.9 and 1.”’ He also observes that when the tempera- 
ture of the condensing plate is raised, the coefficient decreases, but “‘it appears 
that the effect is only great enough to be noticeable at temperatures above 
300° C.” 

The re-evaporation or restitution of silver atoms incident on a substrate 
during vacuum deposition has also been established qualitatively by the 
“‘reflection’”’ experiments of Aziz and Scott (1956). In the present work the 
condensation coefficients were accurately determined for silver with substrate 
temperatures from —30° C to 175° C and for gold and copper with substrate 
temperatures from 25° C to 160° C. The measurements consisted of finding 
the thicknesses of the deposit on one substrate which received the atoms 
directly evaporated from the source and on a second substrate which collected 
the re-evaporated atoms from the first. The two substrates are referred to as 
the reflector and collector respectively. 


EXPERIMENTAL 


The evaporations were carried out in a 12-in. evaporating unit equipped 
with a liquid-air trap. Carefully cleaned microscope slides served as the 
reflector and the collector. The two slides were placed in a holder so that they 
were parallel, 1 cm apart, and inclined at 45° to the direction of the incident 
beam as in Fig. 1. The slides were mounted inside a double-walled cylindrical 
radiation shield which provided a constant-temperature chamber within 
the evaporating unit. A cooling coil of copper tubing and an electrical heating 
element inside the shield were used to obtain the desired temperature. For 
temperatures down to —30°C, alcohol cooled by an external bath of dry 
ice and alcohol was circulated through the coil. The substrate temperature 
was determined by two chromel—alumel thermocouples, one below and the 
other above the slides. 

The metals were evaporated directly from a ‘“‘boat-type’’ tungsten filament. 
The stream density of the atomic beam was kept approximately the same 
in all runs with a given metal by maintaining a constant voltage across the 
filament. The rate of deposition at the reflector for silver was about 45 A of 
thickness per second, which corresponds to a beam intensity of 6X10" atoms 


—2 


cm7? sec7!, 
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Fic. 1. Arrangement of the two microscope slides as reflector and collector. 


The thicknesses of the deposited films were measured at the center of the 
reflector and near the edge of the collector. The method employed for the 
thickness measurements was that of multiple beam interferometry as de- 
veloped by Scott, McLauchlan, and Sennett (1950). The typical variation 
of thickness of the film on the collector slide with distance from its edge is 
shown in Fig. 2. Because the reflector and the collector overlapped by 1 cm 
along the sides, the film was of nearly constant thickness for about 3 mm 
from the edge. It was this constant thickness which was used to determine 
the condensation coefficients. 

The main function of the glass slides was to provide mechanical support 
for the films. Under the conditions of these experiments a continuous film 
of the metal is formed on the reflector for a thickness of about 100 A. The total 
deposit on the reflector was usually in the range 1000 to 4000 A. Since the 
deposit on the collector was found to be proportional to the deposit on the 
reflector, it can be concluded that the restitution (i.e., re-evaporation) 
measured is primarily a characteristic of the metal film itself and is not 
influenced significantly by the nature of the substrate. 

A preliminary test was carried out to confirm the expected absence of 
specular reflection of the beam of evaporated atoms. The result showed that 
specular reflection, if present, must be less than 0.1%. 

The possibility that atoms “‘reflected” once might be “reflected’”’ a second 
time at the collector was considered. If this effect were present to a significant 
extent, a correction to the measured thicknesses would be required in order to 
calculate the condensation coefficient. No correction was deemed necessary 
since after the first reflection the atoms will lose a large fraction of their initial 
kinetic energy and the probability of a second reflection will be much smaller. 
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Fic. 2. Variation in thickness across collector slide. 


Geometrical Correction 

To calculate the condensation coefficient from the observed thicknesses of 
the deposits on the reflector and the collector, it is necessary to make a cor- 
rection for the particular geometry used. The required correction factor is the 
ratio of the deposit actually collected to that which would be collected if the 
reflector were an infinite plane uniformly ‘illuminated’? by the incident 
atomic beam. For the geometrical arrangement used (Fig. 1), the factor was 
calculated to be 0.773. 


OBSERVATIONS 


The results for silver, gold, and copper are given in Table I. The measured 
thickness, D, on the reflector, and the corrected thickness, d, on the collector, 
are listed at various substrate temperatures. The condensation coefficient is 
found from the ratio D/(D+d). It can be noted that in all three cases the 
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condensation coefficient decreases with increase of temperature as implied by 
the Frenkel theory: e.g., the values for silver at —28°C and 171°C are 
0.982 and 0.879 respectively. At room temperature the condensation coefficient 
is 0.936, which is within the range of 0.9 to 1.0 given by Devienne. 

The temperature, 7», at which complete condensation occurs was determined 
by trial in order to give a straight line plot between In(D/d) and 1/(T—Tp). 
The final plots for each of the metals are shown in Fig. 3. The values of J are 
200° K for silver, 230° K for gold, and 180° K for copper. 

Experiments were carried out with aluminum but no restitution was found. 
Even for a thickness of 5000 A on the reflector no deposit appeared on the 
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Fic. 3. Final plot to obtain the temperature T>—the temperature below which complete 
condensation occurs. 
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TABLE I 
CONDENSATION COEFFICIENTS OF SILVER, GOLD, AND COPPER 
Corrected 
Substrate Thickness D thickness d Condensation 
temperature, on reflector, on collector, coefficient, 

Metal *K D/(D + d) 

245 4990 91 0.982 

273 4920 233 0.955 

298 2090 142 0.936 

328 4230 311 0.931 
Silver 351 5496 558 0.908 

376 5360 615 0.898 

419 4843 598 0.890 

444 4945 677 0.879 

298 3387 98 0.972 
Gold 357 2940 119 0.961 

430 2817 137 0.954 

298 1020 115 0.899 
Copper 373 1100 233 0.825 

430 920 233 0.797 








collector with substrate temperatures as high as 175° C. Since deposits less 
than 5 A in average thickness can readily be detected, the restitution of 
aluminum atoms must be less than 0.1% at temperatures up to 175°C. 
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MEASUREMENTS OF ABSOLUTE PAIR PRODUCTION CROSS 
SECTIONS AT LOW ENERGIES! 


S. STANDIL AND V. SHKOLNIK? 


ABSTRACT 


Using an improved analysis in conjunction with a relative method for the 
measurement of absolute pair production cross sections, the pair production cross 
sections per atom of lead have been measured for the gamma rays emitted by 
Co® and Sb’, The Z-dependence of this cross section was measured for the Co 
radiation using targets of Al, Fe, Cu, Mo, Ag, Sn, and Pb. The results obtained 
are generally in agreement with other experimental determinations and are, 
along with those of other investigators, represented graphically. 


INTRODUCTION 


Several authors have measured the pair production cross section in lead 
for gamma-ray energies between the pair creation threshold at 1.022 Mev 
and at 3 Mev (Hahn ef al. 1952; Schmid and Huber 1955; Jenkins 1956; 
Standil and Moore 1956; Shkolnik and Standil 1957). In addition, exact 
theoretical calculations exist for lead at gamma-ray energies of 1.53 and 
2.66 Mev (Jaeger and Hulme 1936), and crude extrapolations based on these 
results have been made (Cleland, Townsend, and Hughes 1951; Hahn, 
Baldinger, and Huber 1952). The accuracy of such extrapolations between 
threshold and 1.5 Mev is ~ 20%, and experimental results are usually com- 
pared with the values obtained from the Bethe—Heitler formula (Bethe and 
Heitler 1934), which is not valid in these cases. Establishment of an experi- 
mental pair production cross section vs. energy curve over this energy range 
thus seemed desirable. 

Using three targets (Fe, Cd, and Pb), Hahn, Baldinger, and Huber (1952) 
measured the Z-dependence of the pair production cross section, and Dayton 
(1953) has published similar results obtained with four different targets (Al, 
Cu, Sn, and Pb). We have also made this kind of determination using our 
experimental method and seven different metal targets (Al, Fe, Cu, Mo, 
Ag, Sn, and Pb). 


THE EXPERIMENTAL METHOD AND TREATMENT OF DATA 
The experimental arrangement, shown in Fig. 1, and the general method 
have already been described (Shkolnik and Standil 1957) and in the present 
measurements only minor changes in the geometry were made. The pair 
production cross section per atom of target, ¢, was given by 


(1) oe (2-)(%4) (3) 


‘Manuscript received April 8, 1958. 

Contribution from the Department of Physics, University of Manitoba, Winnipeg, Mani- 
toba. This work was supported by the National Research Council of Canada. 

2On leave from the University of Belgrade, Yugoslavia. 
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Fic. 1. The experimental arrangement (not to scale). T—target; LB—lead; Q—Co® 
source (~150 mc); M—photomultiplier (DuMont 6364); C—sodium iodide crystal (3.5 in. 
diameter by 2.0 in.); CF—cathode follower; LA—linear amplifier; DD—differential dis- 


criminator; CC—coincidence circuit (resolving time 0.35 uwsecond); SC—scaler. 4 | “ 
where r = source-to-target distance, 
Nf = corrected coincidence counting rate (annihilation quanta) of 
the reference (noncalibrated) Na” source placed in the target 
position, 
Ni = corrected coincidence counting rate due to positron annihi- 
lations in the target, 
N = number of target atoms, 
S = number of reference source Na® positrons annihilated per 
second, and 
Q = number of gamma rays emitted per second by the source. 


99 


Nf was corrected for absorption of the annihilation quanta in the Na 
container, the corrections for one- and three-quantum annihilations being 
negligible; Ny2 was corrected for absorption of the primary beam and the 
annihilation quanta in the metal target and its envelope. Here too the effects 
of one- and three-quantum annihilations are negligible (Hulme and Jaeger 
1936; Bell and Graham 1953). 

The corrections for absorption in the lucite container of the Na” reference 
source and the graphite envelope of the target were determined by using 
known absorption coefficients and estimating the average paths in these 
absorbers. The correction for the effect of finite size was determined by 
measuring the probability of detecting annihilations as a function of the per- 
pendicular distance of their occurrence from the counter axis and averaging 
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this over the particular geometry. A point source of Na” was utilized in these 
auxilliary measurements. 

The optimum orientation of the target with respect to the incident beam 
and the counter axis, i.e., the position with minimum absorption, was found 
experimentally (a = 30°) by rotating the target about its vertical axis and 
detecting the coincidence rate for a given target. Neglecting the edge effects, 
this permits a more accurate correction procedure than that previously used. 
The coincidence counting rate for a given target, n>, may be written as 
(2) nis = Cme®", 
where m = target mass, 


B = an “effective’’ absorption coefficient per unit target mass, and 


‘a ini) 
@  ¢=($)., 


C is the coincidence counting rate due to positron annihilations per unit 
target mass, computed for an infinitely thin target, i.e., no absorption is 
involved. The experimental data for the coincidence counting rates and the 
corresponding target masses substituted into equation (2) and subjected to a 
least squares fit gave C and its error for all targets. The term (Ni2/N) in 
expression (1) is then given by 

Nw _ cA 


(4) y oR: 





where A is the mass number of the target material and No is Avogadro's 
number. The thickness-to-diameter ratio of the target did not exceed the 
value 0.1 in any case, hence the edge effects were in fact negligible. A typical 
least squares fit of equation (2) to the experimental data is shown in Fig. 2 
for the case of the Co® gamma rays and lead targets. 

The determination of the S/Q term of equation (1) was straightforward and 
as discussed in our previous paper. Because of the complex spectra of Co® 
and Sb™, these measurements were made using the highest-energy gamma 
rays in each case (1.33 Mev for Co® and 2.09 Mev for Sb"). 


THE RESULTS 


The sum of cross sections per atom of lead for the equal intensity (1.17 
and 1.33 Mev) Co® lines was found to be 


Dd o(k) exp = (0.274+40.025) barns. 


This is in better agreement with the experimental values of Hahn, Baldinger, 
and Huber (1952), (0.272+0.030) barns, and Schmid and Huber (1955), 
(0.274+0.008) barns, than with an extrapolation based on the Jaeger and 
Hulme (1936) calculation, which gives }°c(k) as (0.300+0.050) barns. 
Using the data recently published for the complex gamma-ray spectrum of 
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nb - number of coincidences per hour/10% 





Oo 1 2 3 
m - target mass in grams 


Fic. 2. Observed coincidence counting rate vs. target mass for Co™ on lead targets. The 
solid line is the least squares fit of equation (2) to these data. 


Sb" (Dzhelepov et al. 1957) the sum of the cross sections for all lines above 
the pair production threshold was found to be 


>, o(k) a. = (9.14+0.40) barn/atom, 
where the factor multiplying o() gives the intensity of the line with energy 
k Mev relative to the intensity of the 2.09 Mev line. The result based on 
extrapolation of the Jaeger and Hulme calculations gives for the above sum 
of weighted cross sections the value of (9.60+0.20) barn/atom. 

Our data previously published for the gamma radiation from Zn® on lead 
targets have been subjected to the treatment herein described. This gives 
a value for the cross section of (13.4+1.0) millibarns per atom, which is 3% 
lower than the result previously reported. The difference is due to an original 
overestimate of the average path in the target. 

The data for the gamma rays from Na* (Standil and Moore 1956) have 
been similarly treated and this yields a final value for the pair production cross 
section in Pb at 2.76 Mev of (3.4340.70) barn/atom, in agreement with the 
result of Schmid and Huber (1955). The large error of this result is inherent 
in the method then employed. 

Fig. 3 shows these results, together with those of other authors, joined by a 








CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 



















10 
5 
iy 
wt 
ze t 
1.0 fF 
0.5 F 
" r experimental 
= 
% 
5 
2 
£ [ 
Cc 
° 
3 
0,1 
= ail 
: I 
2 i 
° 
c 
° [ Bethe - Heitler 
re .05 
.01 
-005 
.001 
05 J 5 1.0 5 10 
(k- 1.022) - energy. above threshold in MeV 
Fic. 3. Pair production cross section in lead vs. energy above threshold. The explanation 


of symbols for experimental points is as follows: gj—Shkolnik and Standil (1957) (1.12 Mev, 
Zn®); W—Jenkins (1956) (1.12 Mev, Zn®*); O—Hahn et al. (1952), Schmid and Huber (1955), 
and this experiment (1.17 and 1.33 Mev from Co®); +—Jaeger and Hulme (1936) (calcu- 
lated); X—this experiment (1.70 and 2.09 Mev, Sb"); @—Hahn et al. (1952) (2.20, 2.62, 
3.20 Mev, Th (C+C”)); O—Schmid and Huber (1955) (2.76 Mev Na*™); &—Standil and 
Moore (1956) (2.76 Mev, Na**). 
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smooth curve and thus represents the variation of the experimental pair 
production cross section in lead with energy above the pair threshold. The 
values of o(k) for the complex radiation of Sb‘, Co®, and Th(C+C”) have 
been estimated from a knowledge of the sum of weighted cross sections and 
the general shape of the curve in Fig. 3 as determined by other points. The 
cross section for Z = 82, obtained from the Bethe—Heitler formula, is also 


plotted for comparison. 


Z-DEPENDENCE MEASUREMENTS 

The Z-dependence of the pair production cross section and radiation emitted 
by Co® was measured using thin targets (Al, Fe, Cu, Mo, Ag, Sn, and Pb) 
of equal diameter. For each element at least seven target masses were employed 
and the observed coincidence counting rate vs. mass curves least-squares 
fitted with equation (2). These measurements being of a relative nature, 
no corrections for absorption in the same target container were necessary. 
The results obtained are in agreement with a Z-dependence of o, given by 


(5) a ~aZ*+bZ', 
or 

2 Gd 2 
(6) az: 1+dZ ° 


A plot of ¢/aZ? vs. Z? is shown in Fig. 4. The slope of this line gives a value of 
d which is compared with the results of the data of other authors in Table I. 


TABLE I 
VALUES OF THE FACTOR d FOR THE RADIATION FROM Co? 


Author d.104 


Hahn, Baldinger, and Huber 


(1952) 1.66+0.30* 
Dayton (1953) 1.59+0. 13* 
Schmid and Huber (1954) 1.55+0.10 
West (1956) 1.81+0.30* 
Present experiment 1.53+0.03 





*The errors in these values have been estimate1 by the present authors. 


CONCLUSIONS 


An experimental pair production cross section vs. energy curve for Pb has 
been determined and the Z-dependence of this cross section for the Co 
radiation has been found. The present measurements are in good agreement 
with those of other authors with the exception of Jenkins (1956), as is seen 
from Fig. 3. We note too that most experimental values for these cross 
sections agree with those obtained by extrapolation of the Jaeger and Hulme 
calculations. The use of Figs. 3 and 4 enables the determination of the pair 
production cross section for the Co® radiation in any absorber. 
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Fic. 4. ¢/aZ? vs. Z* for Co® radiation. 
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ZEEMAN SPLITTING OF DONOR STATES IN GERMANIUM! 


R. R. HAERING? 


ABSTRACT 


The linear Zeeman effect of the 2p m = +1 donor states is calculated in the 
effective mass approximation. The resulting level splitting is independent of the 
longitudinal mass characterizing the ellipsoidal conduction band energy surfaces. 
This result is valid as long as the Zeeman splitting of the m = +1 states is 
smal! compared to the energy difference between the 2p m = 0 and the 2pm 
= +1 states. The Zeeman pattern to be expected in germanium is plotted as a 
function of the angle between the magnetic field and the (100) direction. 


1. INTRODUCTION 

The effective mass theory of impurity states has been very successful in 
explaining the nature of ‘‘shallow”’ impurity states in silicon and in germanium 
(Kohn 1957a).* This theory is expected to yield a good description of impurity 
states when the electron orbits are very large compared to the dimensions 
of a unit cell of the crystal. In this limit, the perturbation due to an isolated 
impurity is simply a coulomb potential modified by the static dielectric 
constant of the crystal in question (Kohn 19576, 1958). Such a description 
of the impurity perturbation breaks down in the cell in which the impurity ion 
is located (Luttinger and Kohn 1955), since the dielectric constant loses its 
meaning at such short distances. This correction modifies mainly the s-like 
states, which have a nonvanishing amplitude at the impurity center. We 
expect, however, that the effective mass theory will give an adequate de- 
scription of the p-like states in germanium.f 

In the following discussion, we shall calculate the Zeeman splitting of the 
2pm = + 1 states of donor impurities in germanium. We shall consider only 
the splitting due to the orbital angular momentum. Since the Zeeman effect is 
measured experimentally by means of optical (infrared) transitions, the 
splitting due to the electron spin will cancel in the observed energy differences. 

In germanium, the constant energy surfaces for the low-lying energies in 
the conduction band may be assumed to be a set of four ellipsoids oriented 
along the (111) directions. Consequently, the kinetic energy operator in the 
effective mass equation has ellipsoidal symmetry and the impurity eigenstates 
may be labelled by a magnetic quantum number which refers to the axis of 
symmetry. This quantum number corresponds to the m/ quantum number in 
the hydrogen atom. However, because of the lower degree of symmetry of the 
effective kinetic energy operator of the impurity atom, eigenstates which differ 
only in the value of the above quantum number m are no longer degenerate 


1Manuscript received April 22, 1958. 

Contribution from the Department of Mathematical Physics, University of Birmingham, 
Birmingham, England. 

2Now at the Department of Physics, Hamilton College, McMaster University, Hamilton, 
Ontario, Canada. 

*This review article also contains many other relevant references. 

+2p-like states are states which approach hydrogenic 2 functions in the limit (m,)/(m)—> > 
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in energy. For instance, the 2p m = 0 level has a lower energy than the 2p m 
= +1 level. The degeneracy between +m and —m remains, however. We 
shall see later that the above facts result in a considerable simplification of the 
calculation of the Zeeman effect. 

We shall find that the Zeeman splitting depends only on the transverse mass 
mt, which characterizes the ellipsoid in question. Furthermore, the level 
splitting depends only on the component of the magnetic field parallel to the 
axis of rotational symmetry (the z-principal axis) of the ellipsoid, and hence 
vanishes when the field is applied perpendicular to this direction. These 
results are a consequence of the nondegeneracy of the 2bm = 0 and 2pm 
= +1 states and may be expected to break down when the Zeeman splitting 
of the 2p m = +1 states becomes comparable to the energy difference between 
the unperturbed 2p m = 0 and 2pm = +1 states. 


2. FORMULATION OF THE PROBLEM 


In the one-electron approximation, the wavefunction of a bound donor 
electron in germanium satisfies the following equation: 


(1) ity vintue fe = E,® 


In equation (1) V(r) is the periodic lattice potential and U(r) is the additional 
potential due to the replacement of one germanium atom by a (singly ionized) 
donor ion. We shall assume this potential to be given by: 


2 
U(r) = -—-—, 
(r) = 
where « is the static dielectric constant of germanium («x ~ 16.0). 
The effective mass theory may be derived by expanding #' (7 denotes the 
different stationary states) in terms of the Bloch functions W,,: 


(2) @'= Sf d*kAd (kK) Vink. 


zone 


If the impurity state in question has a large orbit and if its energy lies much 
nearer to the minimum energy of the mth conduction band than to any other 
band, one may show that equation (2) becomes 


(3) b= | d°kAl(k)Vix 
all k 


oe i : . . 
and the coefficient A, satisfies the so-called effective mass equation: 


(4) \«=(0)— E,+ 2 (ki: tat) ath asc) 


_ 4me® > = 3u7 A i(k’) 
eee 4a ee eee 
K - " all k’ \k— —k’|** 


Here we have specialized to the case of an ellipsoidal energy surface whose 
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minimum is at k = 0. In germanium we actually have four ellipsoids at the 
zone boundaries oriented in the (111) directions. In the effective mass approxi- 
mation, the functions (3) belonging to different ellipsoids do not mix, however. 
Equation (4) is then valid for each of these ellipsoids if we measure k from 
the conduction band minima instead of from k = 0. We may transform 
equation (4) by introducing: 

. 
(5) Fi(r) = (24) °? J dk eS AL(k). 

all k 
The factor (27)? is convenient so that the functions F; are normalized when 
' is normalized. From here on, all functions F‘ and A‘ will refer to the mth 
band (in our case the lowest conduction band of germanium), so that we shall 
omit the index z labelling this band. With the aid of the relation (5), equation 
(4) may be transformed to become: 


) J ef d° -) he d° e*| i “f 
) T ~ 2m, dx dy saa a =e,F, 


where e; = £;—e,(0). 


3. THE ZEEMAN EFFECT 
When a magnetic field is applied, we must add to the Hamiltonian in 
equation (1) an extra term 


(7) H’ = —5-H- (exp). 


< 


Terms quadratic in HW can be ignored since we are calculating the linear 
Zeeman effect. In order to obtain the Zeeman splitting, we shall need the 
matrix elements of #” between the unperturbed states ®', i.e. we shall need 


(8) (!, Ho"), 


Here ©‘ is the total eigenfunction of the impurity state with energy £;,. It is 
convenient to express the matrix element (8) in terms of a matrix element 
between the modulating functions F‘ and F*. One finds: 


(9) (61 Ho") = 


ee, (» \11,( 2220») +11, 2:—208) + 4, (2P+_2Ps) te) , 
2c mM, ™;, mM, mM, Mm, MM, 


The proof of this result, which is rather plausible on the basis of the similarity 
between equations (1) and (6), is outlined in the Appendix. 

In equation (9), the x, y, and z axes are the principal axes of the ellipsoid 
in question. Let us now consider the 2)m = +1 states. The modulating 
functions for these states have the form 


(10) Ft ~ f(z", pet, 


where p = (x?+y?)'. The important feature of these functions is that they are 
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even in the variable z. On the other hand, the parts of the matrix element (9) 
proportional to H, and H, are odd functions of z. Hence (9) simplifies to: 


Cs prt dep — st 
2m .c (F ’ (xpy ypr) F ) 


2 Ataf ae +) 
= ee Pe 5" : 


Here 6 is the angle between the magnetic field and the 2-principal axis. We 


fe + 
(11) (o ; rt (exp)0*) 


— ee ; 
now note that the operator —7h rr has no off-diagonal matrix elements between 


the functions (10) and that 
sp aa +) (r, -a2 = 
¢ ; hag! = +h; ioe hag = —h. 


When these results are used, we obtain the following expression for the shift 
of the energy levels: 


(12) AEmasi = ¥ 


Note that this result is independent of the longitudinal mass m/ and that 
only the component of the magnetic field along the z-principal axis is effective 
in splitting the levels. This is a consequence of the nondegeneracy of the 
m = 0 level. If this state were degenerate with the m = +1 states (as it is 
for spherical E(k) surfaces), then the terms involving H, and H, would also 
contribute. Furthermore, the perturbation (9) would connect the degenerate 
states m = 0 and m = +1, so that the above simple treatment would not 
be valid. We may thus expect the above result to be reliable as long as the 
splitting AE is small compared to the energy difference between the unper- 
turbed 2pm = 0 and 2pm = +1 states. In germanium this means that 
(Kohn 1957a): 

AE K< 3X10 ev. 


We may estimate a characteristic splitting from equation (12) by putting: 
m, = 0.08m); H = 10‘ gauss; cos 8 = 1. One readily finds that 


AEvea10' ~ 7X10~‘ ev. 


It is also an easy matter to compute the expected Zeeman pattern for an 
arbitrary direction of the magnetic field. The result is plotted in Fig. 1 as a 
function of the angle 6 between the magnetic field and the (100) direction. The 
magnetic field is taken to lie in a (110) plane. In terms of this angle 6, the 
Zeeman pattern due to the four ellipsoids in germanium is given by: 
ehH ehH cos8 


sae RRO . aie aeeaaetadi siege 
i; cos(0+55°) and AE=+ oo." 


~ 
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Fic. 1. Zeeman splitting of{the 2p m = +1 donor states in germanium asa function of the 
angle between the magnetic field and the (100) direction. One of the three pairs of lines is 
doubly degenerate for all angles 0. 


4, DISCUSSION 


The calculation outlined above corresponds to the strong field Zeeman 
effect in the atomic case, since we have assumed the spin-orbit energy to be 
negligible. We may estimate this energy by using the analogy with the 
hydrogen atom. We find that the 2 level is shifted by an amount of the order 
of 

Ex» 
(137x)*’ 
which is completely negligible compared to a typical shift due to the magnetic 
field. Hence, if the above estimate is reliable, we are justified in neglecting the 
spin-orbit energy except in very weak fields. 
It is interesting to note that the Zeeman pattern in Fig. 1 does not depend 
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on whether the conduction band in germanium consists of four or of eight 
ellipsoids. Hence a measurement of the Zeeman effect will not vield additional 
information regarding this question. However, the Zeeman effect could yield 
an independent estimate of the transverse mass, if the effective mass theory 
provides a valid description of this effect. 
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APPENDIX 

Consider a typical term in the matrix element (9). 

(Al) M = (+ yp,"), 


On substitution of the relation (3) this becomes: 
(A2) M = J J d°kd®k' An '(k)A, (k’) + (Wn, vP Vn’). 


We may evaluate the matrix element of yp, between Bloch states by expanding 
the periodic parts of these functions to lowest order in k. 


ik-r r hk - Dain a 
(A3) Vik = €¢ Unk ™ - \ ES dea n'ofs 
h’#n 
where: 
9,\3 
Prin * z On) ty ol inv Jun dt 


and nn = €n'(0)—€,(0) is the energy difference between the n’ th and the 
nth band at k = 
If (A3) is substituted in (A2) one finds to lowest order: 
J thm 0 \ 


(A4) M = fave. A, (kK) +5) k, ag, Aj (k). 


Here x, y, and z refer to the principal axes of the ellipsoid. We have made use 
of the f-sum rule and of the relations 


9 


0°, = 6 he 

Oks Maa 
where: m,; = m,, = m,; mz; = m,. If we now introduce the functions F; 
defined by (5), equation (A4) becomes: 
(A5) M=— fare ihy= a4 pe 
A: M = — |drF,\— ety 

: Myre F \ ‘v9 ee , 
m 


F', yp,F"). 


mM, 
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Similar results are valid for the terms z2p,, yp., etc. When these results are 
added, equation (9) follows. 
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SOME GEOPHYSICAL APPLICATIONS OF THE ELASTICITY 
THEORY OF DISLOCATIONS! 


J. A. STEKETEE 


ABSTRACT 


In this paper some aspects of the elasticity theory of dislocations (E.T.D.) 
which may be of geophysical significance are discussed. The central position of 
A- and B-nuclei of strain in this theory is emphasized. Applications to cracks 
and seismological models are indicated. Some limitations of the E.T.D. for 
applications to geophysics appear in the discussion of strain energy. Finally the 
displacement field at the surface of a semi-infinite elastic medium which contains 
a particular B-nucleus at depth c is discussed. Although the considerations in this 
paper are in several respects incomplete and of a preliminary character it is ex- 
pected that the reader may recognize the analogy of the E.T.D. with potential 
theory and will appreciate the power of the theory to unify several isolated 
problems on cracks and seismological models and to suggest further applications 
which are of geophysical interest. 


1, INTRODUCTION 


It is the purpose of the following lines to discuss some aspects of the theory 
of dislocations which may be of significance for geophysical problems. The 
geophysical questions we have in mind are (i) problems connected with 
fault-plane studies of earthquakes and with such phenomena as the San- 
Andreas fault in California, (ii) problems connected with fracture zones in 
the crust and mantle which are believed to play a significant role in the 
structure of island arcs and certain mountain ranges as proposed in the 
orogenetic theories of J. T. Wilson (1954, 1957) and A. E. Scheidegger (1953, 
1958). 

It is clear that each one of the groups of problems mentioned is so vast 
and complex that one can hardly expect to obtain all the answers from one 
simple theory. Also it will be noticed that the considerations in this paper are 
in several respects incomplete and of a preliminary character. 

On the other hand it is found that the theory of dislocations is in so close 
agreement with some of the images and models used in seismology and in 
more qualitative geophysical discussions that it seems reasonable to expect 
some information concerning these problems from a study of dislocation 
theory. 

It should be made clear at this stage that the term dislocation theory may 
lead to confusion. While the author is unaware of theories of dislocations in 
certain branches of the medical sciences, it is convenient to distinguish in 
physics and mathematical physics between two different meanings. 

The first, which is at present no doubt the more significant one, indicates 
a field of solid state physics where one is concerned with certain types of 
irregularities in crystals. This field has developed considerably since 1934 

1Manuscript received May 20, 1958. 
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and a vast and rapidly expanding literature exists dealing with problems of 
slip in crystals, solidification, work-hardening, etc. (cf. Friedel 1956). More 
recently also attention is being paid to the connections with the phenomeno- 
logical theories of rheology (plasticity, non-linear elasticity, etc. (cf. Grammel 
1956)). 

The second, more restricted meaning indicates a branch of the mathematical 
theory of elasticity, which has been developed in the beginning of this century 
by a number of Italian mathematicians—scientists of whom especially 
V. Volterra (1907) and C. Somigliana (1908, 1914-15) should be mentioned. 
While the original investigators used the word “‘distorsioni’’, the term ‘‘dis- 
location” in this connection was introduced by A. E. H. Love (1944). The 
early investigators were chiefly concerned with states of stress in multiply 
connected bodies (cf. Volterra 1907) but it can be seen easily that this theory 
of dislocations performs for the Navier equations of elasticity the same task 
which the theory of surface distributions of poles and dipoles performs for the 
Laplace equation. To distinguish the second meaning from the solid state 
physics it may be convenient to call it the elasticity theory of dislocations 
(E.T.D.). 

While this separation in two theories seems useful for our purpose it should 
not be imagined that the theories were developed independently; on the con- 
trary, the progress in solid state physics was, and still is to some extent, de- 
pendent on the results of the E.T.D. 

In this paper we shall be concerned only with the E.T.D. There are, no doubt, 
in the dislocation theory of solid state physics several ideas and suggestions 
which may also find application in geophysics. However, a large number of the 
problems arise from explicit consideration of the crystal lattice, and this type 
of problem seems to be of little direct importance in geophysics, where the 
length scales are so large compared with the atomic distances in a crystal 
that the continuum assumptions need not be questioned. 

It is expected that the E.T.D. may be useful in geophysics for several 
reasons: (i) it will provide quantitative answers to some questions, (ii) it 
unifies certain methods and results already in use, (iii) it introduces the clarity 
which is characteristic of a mathematical model. 

To obtain this clarity, however, a price has to be paid; to make the situation 
amenable to mathematical treatment simplifications have to be made, and as 
is usual in mathematical physics only those features can be retained which 
are believed to be of major importance. For our case this means that we shall 
omit the curvature of the earth, its gravity, temperature, magnetism, and 
non-homogeneity and shall concern ourselves with an infinite or semi-infinite 
homogeneous and isotropic medium where the laws of the classical theory of 
elasticity hold. While it is no doubt true that in this way several of the in- 
teresting features are discarded it has to be admitted that we gain in simplicity 
and clarity so that we shall be able to investigate certain features with great 
precision. 

A further restriction we impose is that only static problems will be con- 
sidered. This implies that, for example, the waves and other transient features 
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associated with earthquakes will not be discussed. The formulae we use can 
be generalized to the dynamic case but the complexity of the mathematical 
expressions increases considerably and it seemed therefore reasonable to 
consider first the more elementary static case. It should be pointed out that the 
solutions in the dynamical case are in general composed of at least two parts, 
the P- and S-waves. In the static case these two parts merge and become 
undistinguishable. It is therefore necessary to take a certain amount of care in 
extrapolating results obtained here to the dynamic case. In particular it 
should be noticed that there are in the following pages several cases where 
results obtained for static situations are confronted with seismological material. 
The conclusions drawn in such cases can only be considered as tentative. 
Dynamical calculations will be necessary to confirm these conclusions. 

The further content of this paper is divided in eight sections. In Section 2 
we discuss the general theory in some detail and point out the significance of 
the A- and B-nuclei for the theory of dislocations. We further draw attention 
to the fact that the models used by seismologists fit in a natural way in the 
framework of the elasticity theory of dislocations. Section 3 deals with 
Galerkin vectors. We prove the equivalence of the static nuclei B and D 
and extend the use of the Galerkin vector slightly further. In Section 4 we 
recompute the displacement field and stress field of a Griffith crack to show 
that this crack may be considered as a special type of dislocation. In Section 
5 some properties of B- and C-nuclei are investigated. An additional argument 
is given in support of certain conclusions of Honda while a feature emerges 
which may throw some light on experiments of Press. 

In Section 6 the strain energy of the dislocation is discussed and the funda- 
mental theorem of Colonnetti is proved. Some of the restrictions of the theory 
of dislocations in its present form are pointed out. In Section 7 the displacement 
field at the surface of a semi-infinite medium which contains a certain B- 
nucleus at depth c is discussed. In Section 8 some general remarks on the 
elasticity theory of dislocations are made. 


2. FUNDAMENTAL CONSIDERATIONS 


Geologists and geophysicists distinguish faults usually in thrust faults, 
transcurrent faults, and normal faults (cf. Anderson 1942). The distinction 
depends on the orientation and the displacements of the fault faces with respect 
to the surface of the earth. The essential elements these faults all possess 
are the two faces of the fault which were created originally by fracture of the 
material. After the fracture, or while the fracture propagated, the fault faces 
suffered, in addition, a displacement with respect to each other. It is usually 
assumed that largely the same process takes place in the focus of an earth- 
quake (Byerly 1955; Scheidegger 1957). 

In the E.T.D. the following situation is considered; in an infinite or finite 
elastic body, which may be unstrained and at rest, the following largely 
fictitious process is performed: imagine an open surface = which may be 
situated entirely in the interior of the body. A cut is made over 2, and the 
two faces of the cut, which we denote as =+ and =~, are deformed in different 
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ways by applying some force distributions to them. If these force distributions 
are in static equilibrium the deformed state shall again be an equilibrium 
state and the situation created in this manner we call a dislocation over 2. 

If the body is already strained to begin with, the process here described 
has to be extended in this way that, while the cut is made, forces equal to the 
original forces acting across © have to be introduced on the two faces of the 
cut to maintain the original equilibrium state. 

It is clear that this process is largely identical with the geological process 
which is supposed to have taken place at a fault. In particular if the time 
in which the whole process took place is very small, so that inertia effects 
become of importance, it corresponds roughly with an earthquake. 

To make the description more exact we denote the components of the 
elastic displacement vector by u; (¢ = 1, 2, 3) and the components of the 
stress tensor by 7,; where 


(2.1) Tiy = NOG; Une tu(uij+uy;,i), 


\ and uw are Lamé’s constants, 6;; is the Kronecker symbol, and u;,; = 0u;/dx; 
where x; (7 = 1, 2, 3) is a Cartesian coordinate. 

Following Sokolnikoff and Specht (1946) we denote the force per unit 
area on a surface element as 

, 
(2.2) Ti; = Tips, 
where v,’s (j = 1, 2, 3) denote the direction cosines of the normal to the surface 
element. 

The dislocation is then determined by the shape of the surface = and by 
the discontinuity Au; in the components of the displacement vector across 
the cut, that is 
(2.3) Au, = u;*—4uz,, 


where u,;+ is the displacement vector for a point P+, originally being in P 
on 2 but now situated on 2+, and where u;~ is the displacement vector for the 
corresponding point P-, which is now on =~. It is clear that the edge o of 2 
shall be a singularity in general; Au; is not necessarily zero at the edge, and in 
that case the displacement components are not uniquely determined there. 

We shall restrict the deformations of 2+ and =~ to infinitesimal deforma- 
tions. The requirement that the final state of the body with a dislocation over 
> is an equilibrium state can then simply be written as 


(2.4) AE fT th 


where the forces are taken in the direction of the outside (or inside) normals 
to both 2+ and =~. It should be noted that the outside normal to =* is inside 
normal to =~ and vice versa. We shall, following some other authors (cf. 
Eshelby 1951-52), distinguish between two types of dislocations: 

(i) Volterra dislocations (also called dislocations of Volterra—Weingarten), 


(ii) Somigliana dislocations. 
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In the first case the discontinuity in the displacement components is pre- 
scribed to be 
Au, = ujt—uz = U;t+Qy x, 
(2.5) Qy = —MX4, 


where U, and Q,; are constants. The relation (2.5), which is the well-known 
Weingarten relation, states that the discontinuity Au, across = should be of 
the type of a rigid body displacement. In solid state physics one usually 
considers the case 2,,; = 0, and U; is then called the Burgers vector. 

In the case of a Somigliana dislocation Au; can have any form, provided 
the forces which maintain the dislocation do not violate the relations (2.4). 
For a Volterra dislocation, the relations (2.4) are automatically satisfied 
since a rigid body displacement does not generate any strains and stresses. 

In case the body was unstrained to begin with, it will be strained once 
the dislocation is made, as there is a force distribution acting upon =+ and =~ 
to keep the faces in their deformed position.* 

The formulae for the displacement field in an infinite elastic medium strained 
by a dislocation over a surface = was obtained by Volterra (1907). It may be 
useful to outline how this relation is obtained. 

We begin with the reciprocal theorem of Betti, which is a relation between 
two.possible but different displacement fields and stress fields for a particular 
elastic body. Both displacement fields are supposed to satisfy the equilibrium 
equations and the theorem can then be written in our notation as 


(2.6) J fue TP —u, TY YdS = 0, 
s 


where the superscripts inside parentheses refer to the two different displace- 
ment and stress fields. 

One of the two displacement fields is now taken to be the field described 
by the Somigliana tensor u‘(P,Q) (Steketee 1958).t This displacement field 
is generated in an infinite elastic medium by a force F = 82ye, in Q. If, as 
we have supposed, the medium is finite, the forces to be supplied on S to 
obtain the displacement field w4(P, Q) inside the body can be computed in a 
straightforward fashion. Furthermore it becomes necessary to distinguish 
the case where Q is outside or inside the body. In the last case a reciprocal 
theorem is obtained (Love 1944, pp. 245-246), which is, in our notation, 


(2.7) ux(Q) = oe J Ji u(P, Q) T,(P)—u(P) TXP, Q)}dS. 


*It should be remarked that the process of producing a dislocation has been intentionally 
simplified here. One finds in the literature several additional rules which make the operation 
physically more realistic; at the edge of 2, material may be removed if certain edges are 
in the way, while it is also usual to add that the faces are again welded together. While these 
rules may be useful for experiments, they are of no significance for the theory. For the theory 
it is important that the relations (2.4) are satisfied but whether this is done by a force distri- 
bution or by welding cannot be seen in the formulae. As the force distributions more clearly 
reveal that the body must be strained, the possibility of a welding operation has not been 
mentioned. 

tWe shall refer to this paper on several occasions and indicate it further as (St.58). If we 
refer to a particular formula this will be indicated behind the reference. 
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If we apply this theorem to an elastic body which is stressed by a dislocation 
over 2 while the outer surface S is left free from forces, it is clear that the 
integration has to be extended over 2+ and =~ in addition to the outer surface 
S of the body. If we reverse the normal to one of the two faces 2+ or =~ and 
keep in mind that (2.4) is satisfied we obtain easily the expression 


‘ es = P ook off hes 
(2.8) oe sat J Ju tas gaz J Jame Tha, 


which expresses the displacement components in a point Q in the body in the 
discontinuity Au; across = and the displacements of the outer surface S. 
The first integral vanishes if we assume that S recedes to infinity requiring 
at the same time that Q is not at infinity, and = has no points at infinity either. 
In this way we arrive at the formula of Volterra, 


. 
(2.9) u,(Q) = — J faucet (P, Q) »,d3. 

The discontinuity Au; has not been prescribed and the relation therefore holds 
for both Volterra and Somigliana dislocations. In the case of a Volterra 
dislocation it is usual to substitute (2.5) in (2.9) and write the expression as 
the sum of six “‘distorsions élémentaires” (St.58—2.7). In particular, if the 
3 Q4,’s vanish it can be shown (cf. Burgers 1939) that the displacement field 
does not depend on the actual shape of = but only on the shape of the edge 
a of 2. This theorem, which is analogous to a result for uniform dipole dis- 
tributions over a surface in potential theory, is of great importance in solid 
state physics. 

If we repeat the derivation once again but starting from an initial state in 
which the body is already strained, the same result (2.9) is obtained, provided 
no additional forces are applied on S when the dislocation is made. This 
indicates that the displacement field associated with a dislocation over = is 
independent of the initial state of stress and strain. This result, which is of 
interest for applications to geophysics, will be discussed further in Section 6. 

In the derivation of (2.9) r‘,(P, Q) was considered as the 7j-component of 
the stress tensor generated by a force 8xpe, in Q. It was shown in the preceding 
paper (St.58), and can be traced back to a paper of Love (1903-04), that 
ri(P, Q) may be considered for fixed i and j and up to a dimensional constant 
with dimension Au;/u as the k-component of the displacement vector in Q 
generated by a combination of strain nuclei in P; in the case 7 = 7 this is a 
combination of a center of compression and an additional double force without 
moment in the 7-direction; in the case 7 # 7 it isa combination of two coplanar, 
mutually perpendicular double forces with moments. If we denote by a sphere 
a center of compression or expansion which is a combination of three equal, 
mutually perpendicular double forces without moments, the two combina- 
tions may be schematically represented as they are in Fig. 1, (4) and (B). It 
will be convenient to refer to the two combinations as an A- and B-nucleus. 
It follows that the displacement component u,(Q) in (2.9) may be considered 
as the resultant effect of a distribution of A- and B-nuclei over 2. The magni- 
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tude of the contribution of a particular combination of nuclei in a point P on 
> to the displacement components in Q depends both on the local values of 


(A) i=j (B) i4g 
Fic. 1. 


Au; and the orientation of the surface element d= with respect to Q and the 
direction considered in Q. To clarify this further, consider a surface element 
d> of Y with normal in the x3-direction. Its contribution to the displacements 


in Q is 


(2.10) du,(Q) = 


I C . ‘J 
Se (Au, rist Aus rea t+Aus3 se) dd. 
MU 


It is clear that Aw, and Aus, which are in directions perpendicular to the 
normal, may be called the slip, while Au; represents the discontinuity in 
displacement of the two faces d=+ and d=~ in the direction of their normal. 
The effect of the slip can be described by nuclei of type B while the effect 
of the normal discontinuity is described by a nucleus 4. 

Owing to the transformation properties of 7',(P, Q) the same explanation 
holds from point to point on 2; the effect of a normal discontinuity is repre- 
sented by a type 4 nucleus, a slip discontinuity is represented by a type B 
nucleus. It is then also clear that the effect of a slip fault in our model is 
represented by a distribution of B-nuclei over the fault plane. 

If Au;, which varies already from point to point on 2, varies also with time, 
the effect is similar to changing the strength or intensity of the nuclei on 2 
with time. In that case we are no longzer justified in using the static nuclei 
derived from the Somigliana tensor but we have to use the dynamical strain 
nuclei discussed by Love (1903-94; 1944, pp. 304-307). With the introduction 
of these dynamical nuclei we also introduce the two types of elastic waves and 
the situation created in this way is essentially the model of an earthquake. 
It has to be pointed out that a detailed investigation of the dynamical case 
has still to be undertaken; while the generalizations mentioned here seem 
reasonable it is not impossible that in the general case also other terms appear 
which have to vanish in the static case as a result of the relation (2.4). 

It is, however, interesting to see how the theory of dislocations confirms 
the methods used for describing the mechanism in the focus of earthquakes. 
The focus of an earthquake is small compared with the dimensions of the 
earth so that the focus is usually taken to be a point or, if more detailed 
information is wanted, it is enclosed in a small sphere, the focal sphere. To 
describe what happens in the focus several mechanisms have been investigated, 
all built up from certain strain nuclei (cf. Nakano 1923). It is usually agreed 


= 
# 
« 
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that the three types in Fig. 2 approximate the phenomena observed most 
closely. 


ce ae 
a or 

B Cc D 
Fic. 2. 


It is clear that type B is the dynamical counterpart of nucleus B in Fig. 1, 
which is necessary to describe a slip fault. It was mentioned by Honda (1957) 
that the displacement fields of nucleus D and B in Fig. 2 are the same. A 
very simple proof for this equivalence, but valid only for the case of static 
nuclei, will be given in the next section. 

There has been some uncertainty among seismologists (Honda 1957; 
Hodgson 1957) whether the nuclei D and B or the nucleus C in Fig. 2 is the 
most significant. Some remarks about this problem can be found in Section 
5 of this paper. 

Among other theories for the mechanism of earthquakes there is a theory 
of Ishimoto (1932) (cf. also Byerly 1942) which suggests that the sudden 
motion of magma, ejected from a cavity into the rocks of least resistance, 
is the origin of an earthquake. It seems possible that the results obtained 
here might provide a test for this theory; if the magma suddenly forces itself 
into cracks and crevasses in the wall of the cavity, and if indeed this should 
be the major feature of the earthquake, it is clear that the mechanism is 
essentially described by dynamic A-nuclei. The seismograms should show 
the features of such a nucleus which can be easily investigated with the 
formulae of Love for a dynamic strain nucleus. 


3. GALERKIN VECTORS 
It was shown in the preceding paper (St.58) that the expressions we have 
to consider can be put in a very compact form if Galerkin vectors are used. 
The Galerkin vector F = I’; e; is defined as a vector from which a displacement 
field u = u,@;, satisfying the Navier equations, is obtained by differentiation 
according to a linear second-order differential operator, that is 


(3.1) u = (A—a grad div) I. 


In this expression A is the Laplace operator, grad and div the well-known 

vector operators, and a = (A+u)/(A+2yu). The operator in (3.1) is constructed 

by combining certain operators, which, each on their own, are invariant under 

rotations. It follows that the operator in (3.1) is also invariant under rotations. 
It was previously shown (St.58) that each component I, is a biharmonic 

function, that is 

(3.2) AAT, = 0. 


Further, it is trivially obvious that two displacement fields will be the same 








1176 CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 


if the components of their Galerkin vectors are the same. Because of the 
rotational invariance of the operator in (3.1), this theorem still holds if the 
components of the Galerkin vectors become identical only after a rotation 
of the coordinate system. 

The converse of this theorem, which states that the same displacement 
fields have the same Galerkin vectors, is, however, not true, as was shown 
by Westergaard (1952). 

We shall use the preceding theorem to show that the displacement fields 
represented by the nuclei D and B in Fig. 2 are identical, by showing that 
their Galerkin vectors are the same. 

We have seen in the preceding paper (St.58) that the displacement field 
generated in an infinite elastic medium by a force F = 87ye, at O is obtained 
from the Galerkin vector 


(3.3) r = rex, 


where 7 is the distance to O. 

The displacement fields of double forces are obtained by differentiation of 
the single-force displacement field. It is then clear that the Galerkin vector 
for the static nucleus B in Fig. 2 (or Fig. 1) is given by 


(3.4) r= “79 or = r 1 Ge. 


Rotation of the axes over an angle of 45° around the x;-axis gives us the 
transformation formulae 


(3.5) xy = — (%14+-%2), e: = —5 (e,*—e:*), 
i (x2—%1), 5. ~~ (e,*+e,*). 


Transforming the vector in (3.4) with (3.5) gives us on the new set of axes 
(3.6) r= —r,e,*+r2e.*, 


where the subscripts stand for differentiations with respect to x,* and x,*. This 
Galerkin vector is recognized as a combination of two mutually perpendicular 
double forces without moments, which is represented by nucleus D in Fig. 2. 
(It is easily checked that also the sign of the double forces is correct; the 
double force along e.* compressive in the x2 direction, etc.). The Galerkin 
vectors (3.4) and (3.6) which are the Galerkin vectors for the nuclei B and D 
respectively are therefore the same and it follows from the theorem stated 
before that the displacement fields also are the same. 

In the preceding paper (St.58) we have also indicated that the Somiglana 
tensor u’(P,Q) can be obtained from the Galerkin tensor 


(3.7) M(P, Q) = dir = bur, 


where r is the distance PQ. It was shown further that the displacement field 
r(P,Q) with i and j fixed can be obtained from the Galerkin tensor 
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(3.8) P(P, Q) = Nbi5 r tude rj t+b 5 ¥ 3), 


where 7 and j have to be kept fixed. 

It may be convenient to extend this process a step further and to consider 
also the Galerkin vector I,(Q) of the displacement field generated by a dis- 
location over 2. If the differentiations under the integration sign may all be 
performed, and in general this will be the case, we can write 


1 
(3.9) T,(Q) = sae J Jam I',(P, Q) »,d. 


This relation shows a formal analogy with the expression (2.9). Differentiation 
of (3.9) according to (3.1) will give in general the displacement field (2.9). 
While the Galerkin vector has little interest for its own sake, the relation (3.9) 
may be useful as the integrations in (3.9) will be usually simpler than the 
integrations in (2.9). 


4. THE GRIFFITH CRACK 

In this section we consider a two-dimensional problem, which has received 
a certain amount of attention in the literature. Its direct geophysical signifi- 
cance is only small. (It is mentioned in Anderson (1942) in relation to the 
fracture of rocks.) 

Some thirty years ago Griffith (1921-24) put forward a theory of fracture 
for materials like glass based on the assumption that the surface of the material 
contains a large number of cracks. To compute the stresses and strains associ- 
ated with such a crack, Griffith used an analysis of Inglis (1913) in which is 
considered the two-dimensional problem of the equilibrium of an elastic 
medium containing a thin elliptic crack which is being opened by a uniform 
pressure. The solution of the problem is in this case obtained by using the 
confocal coordinates of Lamé, in which the crack is the limiting ellipse. The 
same case, with some generalizations, has been recomputed twice by I. N. 
Sneddon, the first time (Sneddon 1946) by using a stress function due to 
Westergaard, the second time (Sneddon and Elliott 1946) by using a Fourier 
transform method. Green and Zerna (1954) solve the problem with special 
complex variable methods. 

Here we shall show that the Griffith crack can also be considered as a type 
of Somigliana dislocation. It is our main concern to show that this crack 
problem, and therefore also other crack problems (Sack 1946; Starr 1928) can 
be put in terms of dislocation theory and we shall restrict ourselves, therefore, 
to a reconstruction of the formulae for the displacement and stress field. For 
further details and discussions on the physical significance the reader may be 
referred to the literature which has been mentioned. 

The problem under consideration is two-dimensional and one may therefore 
begin with the computation of the two-dimensional Somigliana tensor. It 
turns out to be slightly easier to consider the two-dimensional crack as the 
limiting case of a three-dimensional crack which stretches in one direction 
to infinity. In this way we define the Griffith crack as a dislocation for which = 
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is the strip |x1|<c x2. = 0 = |x3|< © while the discontinuity in displacement is 
given by 


(4.1) Au, = Au; = 0 Aus = Po. /(c*—x?), 
ba 


It will be found that fo in (4.1) is the pressure in the crack. 

The normal to > is clearly in the x2-direction and the general expression 
(2.9) from which the displacement field generated by the crack is determined 
reduces to 


2 ogi taees 
(4.2) u,.(Q) jai a 


e+e 


*+00 
dx | dx3 Aur 759(P, Q), 


—a 


where P denotes a point (x1, 0, x3) on 2 and Q a point (y1, ye, 0). The choice 
y3 = 0 leaves the generality of the problem unimpaired. 

The stress 732(P,Q), which can also be considered as the displacement 
vector generated by a nucleus of type 4 (Fig. 1) has to be obtained from the 
Somigliana tensor (St.68—2.1). We find 


(4.3) r22(P, Q) = 2d(1 aoil . )a-+tuba( : iti 

where’ r = V{(yi—%1)? +y22+-;?} while we further remember that 
or or i 

(4.4) r,= i — iy. =-—r?r 


It is clear from (4.2) and (4.3) that the calculation of the displacement field 
requires the evaluation of six integrals. Two of these integrals vanish, however, 
as their integrands contain derivatives to x3. It follows that u3(Q) = 0 and 
the problem is then clearly a case of plain strain. The four remaining integrals 


are 
ato ate ——- 1 
Dei j ax, | tM e=a(!).. 
+00 +e -—3—-—-5- 1 
Io = J dx3 dream 7) 
- = = F 
(4.5) = 
®+a e+e -—5——_—5- 
foe J dx, | dxynJ/ (c’ —x1")r, 221, 
®*+a atc —3——~5— 
I; = J ax, | dxyvf (=x Jr. 


With this notation the two remaining displacement components can be 
written 
Po 
ui(Q) = ~-—3 {A(1—a) Ii — nal}, 
; trap 
(4.6) 


w2(Q) = Ps {u(Qat1)Ie— nals}. 
4 rap 
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The integrals in (4.5) can be reduced to some different forms if we use the 
well-known expression 


. ‘ ae dx 9 2 
a —- = —? te. fa 
(4.7) lim J has a me 2 2 Inv/ (x1 +x.) 


One obtains then easily 
h=-2) dnV@-m)* aS” 
pe eo ee ee 
I, = af dei (x1) “3, 
a+e =e ai +c ; 
~—* 2 | dxin/(C =x) * i af avec —x}) 
I, = sf dxin/(C° —x’)- 42 f apes 


Substitution of (4.8) into (4.6) and some further elementary reductions 
enable us to put (4.6) in the form 


(4.8) 


ui(Q) = Be {—(1—a) ryi+[y2"+ (l—a@) (c°—y1")] Ts 


2 mya 

(4.9) +(1—a)yi(r>—c") Ie}, 
ux(Q) = <P { —(1-a) xy2— (2—a) yxy2 Ts + [(1—a)r*-+e"lye Toh, 

21 ua 

where 
po x1dx1 

= J. Pym tr)Ve—m)’ 

(4.10) ome v1 ; 


ete 
dx, 
I, = ee ee 
J. (X11 —2yixitr’)V/(c —x1°) 
and r? denotes now 7? = y,;2+y.". To evaluate the two remaining integrals 
we introduce the complex variables 


(4.11) 2 = ytiyo = re e Z = yi—iVve = re 


It is then easily found that 


os ef. yy _ = mg) ‘a 
(4.12) ow lVe—8)~ VE  ~ y INT EHS 

4 = teh} See it 

6 Qy2 (V/(2?—e’) Jt || anes i 


Using the results of (4.12) we can finally put the displacement components 


in the form: 


0 inti LV = ‘= - \ 
ui(Q) = 2 \— (1-a)yt 52 1 BYES B ws ale Ate, 
(4.13) 7 : . <3" 


0 - — an 
u2(Q) = : \—-(-a)yt45 y2 BY—"S* B +54 {> 


2ua 
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where 
hai iad eae pede cos i ee 2 
(4.14) Vv (2°—c") + me Ve —c’) + VE —c’)’ 
a asi B= 3-5 - 
V(e-c) VRC)’ V(e—c) Ve 


Once the displacement field is known the strains and stresses can be found by 
straightforward differentiations. It is easily computed that 


711(Q)+722(Q) = po(Bt—2), 
(4.15) 711(Q) —T22(Q) —ipoy2c2 C, 
712(Q) = 3 poc® y2 Ct 


ll 


where C+ and CW stand for 
mee _ = 1 1 
4. 1 + =, C 6 SS eae 
( 6) C afte = iy + VJ (F—c ey VJ (s°—c’*)° V (#"—c’)* 


To compare these relations with the formulae of Sneddon we put 


- rt) 8 
(4.17) z-c=ne', etc =roe"". 





Fic. 3. 


With these notations and (4.11) the relations (4.16) can be put in the form 


711(Q)+722(Q) = 2peh 7, 008 (0— 36,— 40. saat : 


3/2 
(4.18) 711(Q) —T22(Q) — 2po- (2) (0:+62), 


ei 
p 0" ES sin@ cos 3 (+0. »). 


11" 


712(Q) 


It is easily seen that the relations (4.18) are the same formulae which appear on 
p.428 of Sneddon (1951). Sneddon has taken the crack along the x-axis, 
and the angles 6;, 6, and 6. are measured from the normal to the crack in his 
case. It is immediately seen that for points on the crack we have 


(4.19) 711(Q)—722(Q) = 0, = 711(Q) +722(Q) = —2po, 


which shows that fo is the pressure in the crack. 
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5. DOUBLE FORCES WITH MOMENTS 


In the preceding section we discussed a particular distribution of nuclei of 
type A. Here we shall be concerned with the nuclei B and C. It has been re- 
marked at the end of Section 1 that there is some uncertainty among seis- 
mologists whether the mechanism in an earthquake is best approximated by 
the equivalent nuclei B and D in Fig. 2, or whether type C is the most proper 
representation (cf. Honda 1957); while most of the Japanese investigators 
favor mechanisms represented by B- or D-nuclei, there seem to be several 
indications from Russian work which favor C-nuclei (Hodgson 1957; 
Scheidegger 1957). 

If the nuclei B and C are compared, a significant dynamical difference is 
noticed; while the nucleus B is in equilibrium, the nucleus C is not. If we 
compute the displacement and stress field associated with a static B nucleus 
and, moreover, calculate the resulting force and moment on a part of the 
elastic medium containing the nucleus, it is found that the resulting force 
and moment both vanish. If the same is done for a static C nucleus it is found 
that equilibrium of the medium is maintained only by a moment on the 
surface of the medium containing the C nucleus. 

In the case of dynamic nuclei the situation is more complicated and we 
cannot give such a simple argument. Firstly the two couples of the B nucleus 
may vary with time in different ways. In this case it may be preferable to speak 
of two different and mutually perpendicular C nuclei, but it is clear that this 
phase difference will leave, in general, a resulting moment which disturbs the 
equilibrium. Secondly, if we consider an elastic body which is different from 
a sphere with the nucleus at the center the P- and S-waves reach different 
points on the surface at different times. If we integrate the forces over the 
surface at a particular instant it is clear that there will be no equilibrium in 
general; parts of the surface may still be undisturbed while other parts have 
already been reached and strained by the P- and S-waves. 

To associate this result with the mechanics of earthquakes we recall that 
the information, from which conclusions concerning earthquakes are drawn, 
is mainly the sign and time of arrival of the first P-, S-, and other phases 
which supposedly left the focus at the initial moment of the collapse. This in- 
formation may be supplemented with visual observations and with results 
from triangulations before and after the earthquake. 

From the idealized situation in the theory of dislocations we have found 
that the equilibrium state of a slip fault is described by a distribution of 
static B-nuclei. It must be considered improbable that the equilibrium state 
can be represented by a distribution of static C-nuclei as an additional 
moment will be required which has to be supplied at the surface of the medium 
(in our case at the earth’s surface) to maintain equilibrium. The statement 
has been qualified as improbable since one can imagine special configurations 
in which some C-nuclei will not disturb the equilibrium, while it must also 
be remembered that the earth is not at rest but performs complicated motions 
which may accelerate or decelerate slightly as a result of some unbalanced 
moments inside the earth. 
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If we omit, however, these special configurations and neglect also the motion 
of the earth we conclude that the slip fault in the focus of the earthquake 
can only be described by a distribution of B-nuclei once equilibrium has been 
established. The information obtained from triangulations is therefore expected 
to favor the model with the B-nucleus. 

The first arrivals of P- and S-waves give information which concerns a 
non-equilibrium state, for the occurrence of the earthquake is essentially 
the collapse of equilibrium. So far as this information is concerned there is no 
a priori reason to prefer either the B- or C-nucleus. If, however, the information 
of first arrivals would indicate the C-nuclei as the proper mechanism, it is 
still true that in the final equilibrium position there must be B-nuclei, which 
means that between the initial collapse and the final equilibrium state a second 
couple has to be supplied at the focus. The most reasonable way to generate 
this second couple is a rotation of the fault plane, giving compressions in 
certain quadrants and dilatations in other quadrants, all together resulting 
in the required couple. 

We arrive, therefore, at the following suggestion: if the information from 
seismic waves favors a mechanism with C-nuclei it must be assumed that in 
addition to slipping of the fault faces, a rotation of the fault plane takes place 
before equilibrium is re-established. If the information from seismic waves 
favors a mechanism with B-nuclei it must be assumed that slipping of the fault 
faces over each other will suffice to reach a new equilibrium state. It follows 
also that it must be possible to obtain a mechanism which favors two mutually 
perpendicular C-nuclei of unequal strength, which reach equality as a result of 
rotation of the fault plane at the time that equilibrium is re-established. 

To illustrate some of the differences between the B- and C-nuclei we shall 
consider some aspects connected with the two nuclei. We begin with a detailed 
analysis of the displacement and stress field associated with the B- and C- 
nuclei. Starting from the Somigliana tensor u (P,O) (St.58—2.1) we obtain 
easily for the displacement field of a double force with moment 


us (P, O) = - wpa ehecioa 
_ r 
(5:1) 


a1, 0) = - atu (2). terase 


Expanding these expressions gives us for the C-nucleus 


9 Xo eg .. 
u(P,O) = (2-a) “3+3a > ae 
9 + x Xe" 
(5.2) us (P,O) = — a-3+3a —*3* , 
r : 4 
2, NX 1XoN: 
us (P, O) = 3a s & 


and for the B-nucleus 








es 
a 
ze 
& 
S 
& 
& 

2 
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1,2 2,1 , X21 2 x1 X2 
uy +u; = 2(1—a) —3+6a—-, 
r r 
ans 1,2 2,1 - Ke. Sie 
(5.3) us’ +uz) = 2(1—a) pt ba can 
2, 2, . Xion: 
us +uz) = 6a. 
Substituting cylindrical polars by putting 
(5.4) x1 = pcos, X2 =p sin8, X%3 = X3, 
we obtain easily for the C-nucleus 
) {ae . 
up, = )(1l—a) 3t+za 3 sin20 
p \ ( ) r? 9 r ’ 
rae p ‘ 
(5.5) ue = — pt (l-a) ‘ cos2é, 


us = ~a—-s sin20. 
r 


For the B-nucleus we obtain in the same way 


3 
up = 1201 —a)Pa+3aes t sin20, 


(5.6) uy = 2(1 —a) Py cos26, 


X3 . , 
Uz; = 3a p =~ sin26. 
r 


It is observed that u, reaches for both nuclei extreme values for @ = 7/4, 
32/4, etc. In case of the B-nucleus, u» vanishes for these angles. In case of the 
C-nucleus, u is composed of two terms, the first one independent of 6, the 
second similar to u for the B-nucleus. The first term in us» is negative and is 
therefore a displacement in the clockwise direction, corresponding with the 
sense of the generating couple. For the B-nucleus this constant term has 
vanished since the second couple contributes a similar term in anticlockwise 
direction. For @ = 0 we find for the C- and B-nucleus respectively 


(5.7) ug = — ats ; ug = 2(1—a), 4 


We see that the displacements on the x;-axis are in opposite directions, in 
particular with \ = uw we have a = 2/3 and the displacement components 
are in opposite directions and of equal magnitude. 

We compute also the forces and moments with respect to the origin, acting 
upon a sphere which has O as center. We therefore compute first the stress 


. 
tensor and then the components 7, of the force, keeping in mind that 


(5.8) vy = x,/r. 
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For the C-nucleus we find 


3 
\X1%4 X, Xe 
—6u(1—a)= $4 30an* #2, 





sian 
ra = 6u(1ta)™5?— Bay, 
ra = 6u(L— a)? 30a, 
page Guar — 9g Zt 

(5.9) ra = —Gou(1—a)=¥*— 30a! 
Baad 2u(1—a)y—6 2%-+-Gay 
Ty = —4u(1—a)™}—24ay ize 
T: = 2u(1-+2a)"4— ag 
Ty; = Pha EF 

r 


v 
Because of the odd powers of x, in 7; the resulting forces on the surface of 


the sphere vanish, that is 
(5.10) J T,dS = 0. 


For the same reason two components of the moment vanish. There remains 
only the component of the moment vector along the x3-axis, that is 


(5.11) re (x,T. eigen i= Si 3 2u(1+2a) dpkg a dS. 


Substitution of the spherical polars 


9 


(5.12) x= ~r7 sind cosd?, x.=~r7 sind sing, dS = Pr sind dé dd 


reduces (5.11) to the form 


(5.13) Jao f dé sin{2u(1+2a)sin’6 cos’¢+4u(1—a)sin’6 sin?¢}, 
0 0 


which gives the answer 8zu. It is this moment on the surface which maintains 
the equilibrium of the sphere with the single couple in the center. 
Repeating the calculations for the B-nucleus we obtain 
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x ‘S. 
Ti = 12ya- “5 — 60am 


X 1X2 x1x2° 
T22 = 12ya- nae 60ay~ 7 , 


x1 x X20X 
733 = 12u(1—a) “3 — 60a ou 2. : 
X1X% "x 
T23 = 6u(2a—1) ty — 60a - a3 
f 7 Nox; Hos ot 
(5.14) T31 = 6u(2a—1)- 3 G4 Xt x28 . 
r r 
x 24 ; xX re 
T12 = 4u(1—a) ,+6u(2a—1) 5 — 0a —, 
Pr x1 Xe 


7, = 2u(4a— 1)" —48 que , 


qT, = 2a (de 1)2}— 48 aE 


v 
ops XX 23 
I, = —A8 au" 5 
The last three terms represent again the force components acting upon a 
sphere with the nucleus in its center. It is easily checked in this case that 
both the resulting force and the resulting moment vanish. 
We finally introduce in the expressions for the B-nucleus spherical polars 


by means of the relations 


(5.15) x, = rcos#i, %x*.=rsin@cosd, x3; = r sin sing, 


v , 


¥ 
and compute the components 7,, 7», T4 of the force acting upon a sphere 


v 
with the nucleus at the center. The expressions 7;, etc., are obtained from 


Es T, cos6+ ( ie cos¢+ Ts sing) sind, 


(5.16) Ty = —T,sind-+(T2cosé+T3 sin$)cos@, 
y = T; cos¢@— Ts: sin ¢. 
We obtain then from (5.14) and (5.16) 


T, = —(200-+1) SP TOse 


(5.17) Ts = Galt) cose rose 


: cosé sin p_ 


T = —2y(4e—1) —3- 
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It has been suggested (Honda and Masatsuka 1952) that the mechanism 
of a deep focus earthquake is properly described by an impulsive radial force 
distribution in a small spherical cavity of radius a, which is determined by 


(5.18) (F,)rma = Fsin20cos® (Fe)rma = (Fs)rma = 0, 


while the mechanism of a shallow earthquake is properly described by a 
B-nucleus. In a more recent paper, Honda (1957) shows that the wave patterns 
of the two mechanisms are essentially the same. Inspection of the formulae 
(5.17), although obtained for only a static case, indicates the same result. 
For the case \ = yu, we have a = 2/3 and one finds 


(5.19) u(20a+1) = 438/3u, 2yu(4a—1) = 10/8u, 


which shows that the radial force exerted on a sphere with the B-nucleus at 
its center is more than four times as large as the tangential forces; thus the 
effect of the tangential forces may be easily overlooked in practice and result 
in formulae of the type (5.18). 

We shall now consider a two-dimensional case and begin with the compu- 
tation of the two-dimensional Somigliana tensor. There are several methods 
by which this result can be obtained (Love 1944, pp. 208-214). We distribute 
forces in the x2-direction along the x3-axis from — © <x; <+© in sucha 
way that the force per unit length is 8rues. Remembering that 


a+L 
~ . . 1X3 2 2 
e .20 | J een a —~———_—~— = —2 In x1 Xo . 
© oe ~~ V1 +%X2 +%3° Vv‘ pe) 


one easily obtains 


tx 
“i, = | dy3(—a@r.o1) = 2a(x2 In p) 1, 
J —~ 


wr 
bo 
— 
= 
= 
rr 


QO+o 9 
| ay? ra) = —41n p+2a(x2 In p) 2, 
U3 = j | dy3(—a Yr 93) = 0, 


where P= V xP? + (x3—ys)"} p= V (x1? +2"). 


In a similar way we find ui and us and the two-dimensional Somigliana 
tensor. To keep the formulae in this case as simple as possible we normalize 
them to a force of 42 instead of 8rp per unit length. The two-dimensional 
Somigliana tensor then becomes 


(5.22) u\(P,O) = —26% In pta(x;, In p), ;. 
Computing the displacement field of a C- and B-nucleus gives us: 


1,2 v NS. ¢ X1 Xe 
te = =a) -a-+-2Ze 

p p 
1.2 xy X 1X2” 
¢ =—azs 


3+2a 4 
p p 








TO AAT NC CSTE 


IED ACARI A TREE MOTE MELA TRATES OT ONEROUS Pe 


EAE RPI PI TALON PRI STOLEN Om eV 


eae ete 


SPUR ORES aN NE 





I 
fe 
E 


Sm pER 
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a bes = 2(1—a) atte Xi 2 
(5.24) 9 
1X: 
i 


us’ +uz = 2(1-a) * Fs die — 
p 





Substituting polar coordinates x,=pcos6 x2.=psin@ gives us from (5.23) 


b., 
u, = — sin20, 
p 


(25) 
l cos26 
up = —-+(1—a) ————. 
p p 
From (5.24) we obtain 
a. 
yp = ~sim2e 
(5.26) 
cos26 


ue = 2(1—a) — 


It is clear that the displacement fields display the same features we have 
already found for the three-dimensional C- and B-nuclei. 

We distribute now the two types of two-dimensional nuclei in a uniform 
fashion along the x-axis over the interval —1 < y; < +1 ye = 0 and study 
the displacement field so obtained. The displacement field generated by the 
C-nuclei is determined from 


aaa U, ¢* 2 
(5.27) u(P)=— | dy, ux?(P, Q). 
dr =; es 
The displacement field generated by the B-nuclei is determined from 


(5.28) u(P) = uf" dy,{uk*(P, Q)-+u? (P, Q)}, 

which is a special form of the general relation (2.9). The situation considered 
here may in fact be considered as the expression for the displacement field 
of a two-dimensional Volterra dislocation with Au, = U;, Aus = Au; = 0, 
while 2 is the strip |yi/<1 yo=0 = [ys|\< @. The second term in (5.28) 
can be integrated immediately. If we denote with Q* and Q~ the points (1, 0) 
and (—1, 0) (5.28) can be written 


(5.29)  ,(P) = uss uz(P, Q*)—uz(P, O- +f dy, us"*(P, at. 


Physically this result is clear; the couples u71(P,Q) all cancel each other 
inside the interval —1 < y; <+1 in such a way that only two forces at the 
end-points remain. It is then also clear that in this case the nuclei u{?(P, Q) 
are responsible for the discontinuity in the displacement field across the 
x,-axis in the interval |y,! < 1, while the second couples are necessary only 
to maintain the equilibrium. Evaluation of (5.27) using (5.23) gives us the 
displacement field 
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_ Uif atl aml a ‘u-_sct)} 
u,(P) = a ytan = tan . +5 *2 Re Re /f? 
ae Ry od Sy 
uo(P) = on ona, +9” Re Re (- 
Adding the displacement field 
ux(P, Qt) —ux(P, Q) 
to (5.30) results in the displacement field (5.29). That is 
as, : -1_x,+1)\ 
uP) = ~i Ye tan ET tax {e )t 


2a 





(5.30) 








(5.31) 
u(P) = — a4 (1—a) In 4 = Eten: (ah) {- 
In these expressions 

Ri = Vi (x; — 1)?-+-2%}, Ry = Vi (%1+ 1)?+-x27}. 


It is easily checked that u,(P) suffers a discontinuity of magnitude U; when 
crossing the x-axis in the interval |x,;| < 1. This result holds for both dis- 
placement fields (5.30) and (5.31). 

If the stress field is computed it is found that the displacement field (5.30) 
is maintained only by a moment acting along the circumference of any segment 
of elastic material which contains the section |x,| < 1 of the x-axis. In the 
case of (5.31) there is no resulting force or moment acting upon a segment 
which contains the section |x,| < 1 of the x-axis. It is clear that the moment 
which keeps the displacement field (5.30) in equilibrium has the magnitude 
2uU);. 

If we consider the displacement component u2(P) on the x-axis we find 
from (5.30) and (5.31) respectively 


Lie jg Ri 

R,’ 
It is again clear from these expressions that the displacement components 
on the x;-axis are in opposite directions. In particular for a = 2/3 their mag- 
nitudes are equal. While the displacement of the C-distribution corresponds 
with the clockwise sense of rotation of the C-nuclei, the displacement com- 
ponent of the B-distribution is in the opposite direction. Although the results 
here obtained are derived for the statical case only it is not impossible that 
they may throw some light on certain dynamical experiments performed by 
Press (1957) which were reported at the recent I.U.G.G. Congress in Toronto. 
The experiments were performed with a thin disk, in which there was a slit 
along part of a diameter (Fig. 4). If the faces of the slit were moved in the 
way as indicated by the horizontal arrows the first movement of a point P 
on the rim was in the direction indicated by the vertical arrow. Assuming 
that it is again permissible to extrapolate the results from the static case to the 
dynamic problem it is seen that the slit in the experiments should be considered 


(5.32) ux(P) = cba) in 


hh) = = NR,” 





} 
F 





Ir PPT 
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as a distribution of B-nuclei and not C-nuclei. This is borne out by the fact 
that the rim was free in the experiments while the slit was held in a fixed 
position, no rotation of the disk as a whole being allowed. There are, in ad- 
dition to analogies, some differences also between the case of Press and the 


Fic. 4. 


situation we considered here. The thin disk is in general considered as a case 
of plane stress while our case is essentially a case of plane strain. It is, how- 
ever, well known that the one solution can be converted into the other by a 
suitable change of elastic constants. Secondly, our calculations are essentially 
performed for an infinite elastic medium while the disk is definitely finite. 
To counter this argument it may be pointed out that the experiments of 
Press were concerned with only the first arrival of the waves from the slit at 
the rim. Up to that time the waves were still unaware of the existence of a 
boundary and therefore behaved in the way in which they would in an in- 
finite medium. Finally, we once again recall that the result will be definite 
only after a recalculation of the same case with the dynamic formulae. 


6. THE STRAIN ENERGY OF THE DISLOCATION 


We have seen that in a body which is unstressed, stresses are generated if a 
dislocation is produced. To obtain an expression for the strain energy we 
recall that the elastic potential or strain energy per unit volume is given by: 


(6.1) W= 3 Ty Cy. 


For the total strain energy of a body with volume V we have then immediately 


(6.2) W = { wav. 
Vv 


The integral (6.2) can be transformed into the form 


J rumsaV = +f (14s), dV, 
Vv Vv 


where we used the expression ty wy = 0, with wy = 3(uj,;—ui,), and the 
equilibrium equation. The last integral can be replaced by two surface integrals 
if the divergence theorem is used. We obtain 


(6.4) w=3f futas+s ff autaz. 
e Ss z 


(6.3) W= 


tol 
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The first integral vanishes if no forces are applied on S when the dislocation 


is made, and we obtain for the strain energy 


(6.5) 


W = bf f aut. az. 


The integral (6.5) is in most cases an improper integral as a result of sin- 


v 
gularities of 7; on 
evaluation of the strain energy using (6.5) presents therefore a number of 


> and especially along the edge of 2. The straightforward 


analytical difficulties which will require further study. For the time being, 
however, it is assumed that the analytical tools which are used, for example 
in aerodynamics (leading edge singularities, vorticity singularities, etc.), will 
be sufficient to cope with this situation. 

We suppose now that the body is already strained before the dislocation is 
made and let e7;, 77), w* define the initial equilibrium stress distribution in 
the body. If e,;, rij, and w denote as before the strains, stresses, etc., due to 
the dislocation we have for the elastic potential w’ in the final state the 


expression 


(6.6) 


’ 1 * * 
2 a (Tis t+7is) (Cis Hei;), 


* re. * 
ee tw o(rijeis +7 ijeij)- 


The expression between parentheses in (6.6) may be written in the forms 


(6.7) 


1 * * * * 
2 (7 jig +7 ij€i;) = Tijlij = Tijlij- 


The total strain energy in the body in the final state is then given by 


(6.8) 


W' = W+W+ Jiriawsav = Wt w+ | rut saV, 


where we have denoted with W’, W, and W* the volume integrals of w’, w, 
and w* respectively, while we also used the expressions 


(6.9) 


* * 
T ijMig F TijMij = 0. 


The volume integrals in (6.8) can be transformed into surface integrals. 
Considering that 7,; and 77; each on their own satisfy the equilibrium equations 
we obtain in turn 


(6.10) 


W’ 


wtw+f fut as+{ | Au,T* d3. 


wwf f uttas+ | J AurT’; d2. 
s z= 


We suppose that to make the dislocation, forces were applied over = while 


the forces on S were left unchanged. In that case 7; = 0 on S. Furthermore, 
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Aut; = 0 on & as the initial displacement field is continuous throughout the 
body. From (6.10) we then deduce the two relations 


W’ = W'+W, 


jj uT rast] { aut d= = 0. 


This is the Theorem of Colonnetti (1915). 

The last relation states that the work performed by the initial forces over 
> while making the dislocation is equal and opposite to the work performed 
by the initial forces on S when the dislocation is made. 

The first equation in (6.11) shows that the strain energy in a body with a 
dislocation is equal to the sum of two strain energies: (i) the strain energy of 
the initial stress distribution, (ii) the strain energy accumulated in the body 
when we produce a dislocation in the unstressed body, leaving the surface 
S free from forces. 

The Theorem of Colonnetti shows that the strain energy associated with 
the dislocation is independent of the initial state of stress or strain which 
may exist in the body. This result is in agreement with the statement that the 
displacement field in an elastic body that is due to a dislocation over © is also 
independent of the initial state of stress or strain in the body, which was 
mentioned in Section 2. 

If one considers problems of fracture in geophysics, the Theorem of 
Colonnetti seems puzzling, as it is felt intuitively that a fracture should relax 
the stress distribution while it is found that the strain energy is increased 
when a dislocation is made. To put this difficulty in the proper light we recall 
that so far we have considered mainly cases where an elastic body was strained 
exclusively by a dislocation. If boundary conditions were considered we 
prescribed that the surface was essentially free when the dislocation was 
produced, no additional forces being applied on S. With these boundary 
conditions the Theorem of Colonnetti inevitably follows. 

It is, however, possible to select other boundary conditions. One can assume 
that the outer surface S is held fixed while the dislocation is made. It is clear 
that in that case also forces on S have to be applied while the dislocation 
is made, and the displacement and stress field generated in the body will be 
different from the displacement and stress field obtained under the former 
assumptions, even if © and Az, will be the same. It is easily found that with 
the conditions of a fixed boundary the relation (2.8) has to be replaced by 
the expression 


; 1 iam 1 . 
(6.12) “:(Q0) ==— JJ u\T ,dS+=— Sf Au,T% d=. 
Sau Ss Simp z 
where 7; in the first integral is now the force to be applied on S to keep the 
surface fixed. This same relation is also found if the body is already in a state 
of initial stress, but the surface S is held fixed in that initial position when 


(6.11) 
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> dislocation is made. e boundary con ‘ition u; = 0 or is applied i 
the dislocation is made. If the boundary con ‘ition 2 0 on S is lied in 


the formulae (6.10) we obtain 
wwe f Au,T; d2, 


wtw+f at? a5. 
Ss 


In these formulae W can be written in the form (6.5) as the first integral in 
(6.4) vanishes also if uw; = 0 on S. The integral in (6.13) can have several 
different values and W’ is not necessarily greater than W*. The value of W’ 


v 


W’ 


(6.13) 


ll 


will be determined if we know Au; or 7; on Y. In connection with fracture 


G. I. Taylor (cf. Starr 1928) has suggested that over the fractured surface > 
complete relaxation of stresses will take place. We have then in the final 


position for points on ¥ the condition 


(6.14) PAT, at 


With this condition it follows from (6.10) that W’ = W*-—W, 


(6.15) [ faut d= = JJ “iT, dS, 
e z Ss 


where W is given by the expression (6.5). 

It is noticed that with these conditions the formula (6.12) loses a con- 
siderable amount of its usefulness for the calculation of the displacement 
field as Au; is not known over = but has to be determined first. Under these 
circumstances a different formulation of the problem will be necessary. 

With these new boundary conditions many new problems arise. A body 
can be held fixed over part of its surface while other parts are free, or again 
parts of the surface can be subjected to given loads while the dislocation is 
made. Also on = the discontinuity Au; can be prescribed or certain conditions 
of partial or complete stress relaxation can be imposed. Depending on these 
conditions the strain energy will increase or, starting from a given initial 
state, the strain energy may decrease, as for example in (6.15). 

It is clear that for geophysical problems only such boundary conditions 
will be physically realistic which will cause a reduction in strain energy when 
the dislocation is produced. On the other hand we have to recognize that the 
surface of the earth is essentially free and if a dislocation is made under 
those circumstances, Colonnetti’s Theorem shows that the strain energy 
can only increase. 

In this dilemma one may perhaps see an additional argument for the 
generally accepted idea that earthquakes are not world-wide phenomena but 
that their cause lies in local conditions; if one considers the whole earth it is 
difficult to see any restriction on the displacement of its surface; on the other 
hand if a particular region is considered it is easy to imagine that certain 
parts may obstruct or prevent the displacement of others and give rise in this 
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way to boundary conditions of the type in which parts of the surface are no 
longer free, creating in that way the possibility of escaping Colonnetti’s 
Theorem. 

It is clear that the problem of applying dislocation theory to geophysical 
features becomes very complex if the strain energy is required to decrease 
when the fracture is made. In the cases which have been considered so far 
the situation was usually simplified. Housner (1955) applied dislocation theory 
to earthquakes and introduced certain physical and statistical assumptions 
which enabled him to obtain amongst other things an estimate for the relation 
(6.5) and to relate his results to certain formulae obtained by Gutenberg 
and Richter. Rochester (1956) obtained approximate expressions for the 
displacements and stresses in a semi-infinite elastic medium which was stressed 
by a Volterra dislocation over a rectangular surface. The strain energy is 
determined from (6.5). It is then assumed that this is also the strain energy 
which will be released if an earthquake produces the given Volterra dislocation. 
This seems the most reasonable and most simple assumption. As no initial 
stress was given, the condition (6.14) could not be used. The only remark 
which can be brought against this method is that it may be better to impose 
the boundary conditions that the displacements have to vanish far away than 
to assume that there are no forces applied at infinity. It is, however, also clear 
that this will complicate the situation. If the strain energy is computed for 
the examples discussed in this paper and for other cases where strain nuclei 
are distributed in an elastic medium we shall obtain expressions similar to 
Rochester’s procedure, as it is assumed that these strain nuclei are the only 
cause of stress in the body. The paper of Starr (1928) is more satisfactory 
from the point of strain energy but his analysis is not yet entirely correct; 
while the discontinuity Au; across the crack should be computed on the 
assumptions of a shear at infinity and complete stress relaxation over the 
crack, the crack is prescribed as an ellipse in his special confocal coordinate 
system. The analysis of Starr has been taken up again recently by Knopoff 
(1957). For completeness sake it may be mentioned that there is also some 
Russian work on dislocations in geophysics (Vvedenskaya 1956), but the 
author had so far no opportunity to see this paper. 

It should be pointed out finally that the approximations on which our 
theorems are based in the case of initial stress are suitable only if the initial 
stresses are sufficiently small. The stresses in the earth are usually so large 
that these approximations are not justified. The next step is then to consider 
the dislocations as problems in initial stress (Love 1944, pp. 109-111; Jeffreys 
1929, pp. 161-163). 


7. THE DISPLACEMENT FIELD w (P,Q) AT THE SURFACE 
OF THE SEMI-INFINITE MEDIUM 
The considerations in the preceding sections were concerned mainly with 
an infinite elastic medium. If we wish to consider dislocations in the earth 
we have to supplement these considerations with certain suitable boundary 
conditions. We have seen in Section 6 that a proper choice of boundary 
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conditions in addition to a suitable initial stress distribution is necessary 
to obtain a physically satisfactory picture of the strain energy. 

Following other authors (cf. for example Lamb 1904; Pekeris and Lifson 
1957) we approximate the earth by a semi-infinite medium with a stress-free 
surface. The medium is unstressed to begin with and is then stressed by 
making a dislocation while no additional forces are supplied at infinity. We 
have seen in Section 6 that under those circumstances no satisfactory picture 
of the strain energy is obtained, but this disadvantage is considered to be 
compensated for by the simplicity of this model. 

In the preceding paper (St. 58) a Green’s function method was developed 
to deal with this problem. To obtain the displacement field, generated by a 


dislocation, in that case, we have the relation 


(4,1) u,(Q) = JSS. Auw'(P, Q)v, dd," 


which differs only from equation (2.9) in its kernel; the functions riy(P, Q) 
have been replaced by w,(P, Q). A method was indicated for the construction 
of the wi(P, Q) functions and the set wix(P, Q) was explicitly constructed.t 

In this section we wish to consider in some more detail the displacement 
field wh(P, Q) at the surface of the semi-infinite medium, which is given by 
(St. 58—7.20). 

The B-nucleus lies in a plane parallel with the free surface (x; = 0) and 
at a depth c. The displacement components were written in the general notation 
with two fixed subscripts and one free superscript but we shall denote them 
here simply as u; (¢ = 1, 2,3) while we also replace y; by x;. Omitting the 
elastic constant the displacement components at x3; = 0 can then be written 


in the form 





uy = - pare a 7 
(7 2) uo = . 
X 1X2 
us = —7C, 
r 
where 
ean ¢ - 
(73) inti, 
pp 


2 4 6 


B = 6+16°—48-3+-38—12, 
pp p p 


*In (St. 58) we wrote out all the expressions for the case of a Volterra dislocation, but one 
easily verifies that the procedure applies also to the more general Somigliana dislocations. 

tlt may be good to indicate that in Section 7 of (St. 58) some confusion may arise from the 
use of r*,(P,Q) throughout. In particular, in the heading of the section and in the final 
formulae (7.18) and (7.20), r4(P, Q) should be replaced by wh( P,Q) to maintain agreement 
with the notation of the preceding sections in that paper. 
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3 5 
Cc »€ c 
C = 4—18°+26%;—125 , 
p p p 
r = x1 +x", 
p° =¢+r’*, 


It has to be remembered that u; is positive in the direction into the medium. 

To investigate the above expressions further we substitute polar coordinates 
x, = rcosé, x2 = rsin@, and obtain the radial- and tangential-displacement 
components u“, and us in the plane x3; = 0. We find 


u, = ritA+3B} sin 26, 
(7.4) ue = m3 A cos 26, 


U3 = 51 C sin 26. 
For a given depth of focus c and radius p, only the trigonometric functions vary 
in (7.4) and one easily notices the distribution in quadrants; u, and u3; change 
sign on the coordinate axes while m» reaches extreme values, ue changes sign 
at 45°, 135°, etc., and u, and uw; then reach extreme values. 

To study further the displacement as functions of the distance to the 
focus we express p and r in terms of c and the angle ut (Fig. 5). If we substitute 


ct 


Fic. 5. 


p = cos t = C/p it is easily checked that 4+3B, 4, and C have p = lasa 
double root. After simplification, the expressions (7.4) can then be put in the 


form 
p°\/(1—p’) (5+ 14p+p*?—12p*—6p') sin 20 
4, = B P SS 
r (1+ p)° a 
p42 5 ee > 
ass _ 2p V/(1—p') cos 20 
Che) ug = (+p) -* 
aid p’(1—p)(2—5p—12p"—6p*) sin 20 
sali i+p e 


We notice that the displacements at the surface decrease with increasing 
depth. This decrease is fairly rapid as the square of the depth appears in the 
formulae. To obtain an idea of the dependence on p we plotted in Fig. 6 the 
relations (7.5) when c*, cos26, and sin2@ respectively were taken equal to one. 
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Fic. 6. 


The curves for u, and us have few interesting features; the radial displace- 
ment reaches an extreme value near p = 0, 7 and so does up. The displacements 
u, and us keep the same sign for all angles y. The curve for u3 reaches an ex- 
treme value at a value of p which is slightly larger than in the case of u, and 
us, but in addition us changes sign. The expression 2—5p—12p?—6p* = 0 
has a root which is approximately p = 0, 242. The angle ¢ is then slightly less 
than 76°. As tan 76° is about 4 we see that the vertical surface elevations change 
sign on a circle with a radius which is about 4 times the depth. To illustrate 
this behavior the distribution of uw; in the plane x3; = 0 is indicated in Fig. 7. 
The property of uw; which has been found here raises the intriguing question 





Fic. 7. 


whether the same behavior will occur if the dynamic case is studied. It might 
throw some light on the difference between first arrivals at Berkeley from 
shallow and deep shocks in Southern California (Byerly 1955). Further in- 
vestigations will no doubt be necessary here. 
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In several papers Tsuboi has reported the change in the surface of the land 
caused by an earthquake. In particular the horizontal displacements given for 
the Idu earthquake of 1930, which was associated with the Tana fault, show 
a pattern which is not unlike the pattern which we obtain for our resulting 
horizontal displacements (Tsuboi 1932 and 1939). Tsuboi also supplies 
information concerning dilatations, rotations, and several more features. To 
compare this with our model, we shall have to make further calculations. 


8. FINAL REMARKS AND ACKNOWLEDGMENTS 

In the preceding sections some simple cases in the theory of dislocations 
were discussed and an attempt was made to relate them to certain geophysical 
features. While the treatment is fairly systematic from the point of view of 
the theory, the geophysical aspects discussed show a fair amount of diversity. 
This aspect of the theory to order different features into one consistent 
scheme seems especially attractive in a field as geophysics where there is 
such a large scope for speculation. 

The work is of a preliminary character and the conclusions reached are in 
several cases of a very tentative nature as only static cases were considered. 
Further work will no doubt be necessary, to confirm several of these points 


or to prove them wrong. 

It is again a pleasure to thank Professor J. T. Wilson of the University 
of Toronto for his invitation to work during the summer of 1955 on geophysical 
problems. The author is grateful for financial assistance received during the 
summer of 1955 from the California Standard Oil Company and in the summer 
of 1957 from the National Research Council of Canada. 
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THE THEORY OF FORBIDDEN 8-DECAY' 
G. E. LEE-WHITING 


ABSTRACT 


The angular distributions of electrons and of neutrinos and the polarization 
of the electrons emitted by an oriented source in 8-decay of arbitrary forbidden- 
ness are calculated in terms of various angular-momentum coupling constants 
and of certain generalizations of the well-known combinations of electron radial 
wave-functions, Ly, M,, N,,.... The interaction Hamiltonian used encompasses 
all five invariants; neither parity conservation nor time-reversal invariance is 
assumed. A particularly simple procedure for calculating the longitudinal 
polarization of electrons emitted by an oriented source is pointed out. The 
polarization is equal to +v/c, to good approximation, if the two-component 
neutrino theory holds. Explicit results for various allowed and first-forbidden 
effects are quoted. 
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1. INTRODUCTION 

The work done in recent years on the properties of the various angular- 
momentum coupling coefficients—see, for example, Rose (1957) or Sharp 
(1957)—makes possible the formulation of a general procedure for calculating 
the consequences of the Fermi (1934) theory of 8-decay in a transition of arbi- 
trary forbiddenness. The revival of interest in 8-decay experiments induced 
by the suggestion (Lee and Yang 1956) and the subsequent experimental 
discovery (Wu et al. 1957) of the nonconservation of parity increases the need for 
such a formulation. It is desirable to use the parity-nonconserving interaction 
Hamiltonian of Lee and Yang and to include the effects of the possible failure 
of time-reversal invariance. Some work along these lines has already been 
done by Alder et al. (1957) and by Dolginov (1958). Results applicable to 
first-forbidden transitions have been given by Curtis and Lewis (1957) and 
by Berestetsky et a/. (1958) as well as by the former authors. Our purpose is 


1Manuscript received June 10, 1958. 
Contribution from Atomic Energy of Canada Ltd., Chalk River, Ontario. 
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to derive, in a manner both simple and general, results analogous in form to 
the well-known formulae for correction-factors (Greuling 1942 and Pursey 
1951). 

Because the succeeding sections of the paper abound in calculational detail, 
it is thought desirable to give at this stage a nonmathematical description of 
the path to be followed. At the same time the opportunity is taken to point 
out, for the benefit of the reader who is not particularly interested in the 
derivations, where formulae corresponding to different experiments may be 
found. An expansion in spherical waves of the wave-function of an emitted 
electron associated with a given nuclear transition and with the absorption 
of an antineutrino from a given state is calculated by time-dependent pertur- 
bation theory; the perturbation is the interaction Hamiltonian of Lee and 
Yang (1956). The form of perturbation theory used is identical except for 
notation to an adaptation to the relativistic case by Greuling and Meeks 
(1951) of some earlier work by Bethe (1930). 

The asymptotic expansion of the electron wave-function contains matrix 
elements between the initial and final nuclear states of operators composed 
of Dirac operators and products of lepton wave-functions. The derivation of 
a comprehensive formula for these matrix elements appropriate to a {-inter- 
action containing all five invariants is greatly facilitated by the use of a 
notation, devised by Banerjee and Saha (1954), for labelling the eight Dirac 
operators: J, 8, 6, 86, Ys, Bys, a, Ba. Each operator is put into one-to-one 
correspondence with one of the eight sets of values of three two-valued para- 
meters, viz. £, 7, and ¢. The value of & is 0 for the scalars and 1 for the vectors 
(in three-space); 7 takes the value 0 for the ‘“‘nonrelativistic’’ operators (the 
first four in the list), and 1 for the others; ¢ is 0 for those containing an explicit 
8, 1 for the rest. More fundamental properties of 7 and ¢ are discussed later in 
Section 1. It is found necessary to introduce a particular representation of the 
Dirac operators by means of 4 X 4 matrices. With the aid of the parameters 
(é, n, ¢) one can obtain a simple formula for any element of any of the represen- 
tative matrices. The use of this formula enables one to write the interaction 
Hamiltonian in manageable form. The lepton covariant in the matrix element 
for the transition can be reduced to a sum of products of functions of r and of 
electron energy with spherical harmonics. The dependence upon electron 
energy can be approximately factored out of the matrix element. Next, follow- 
ing Spiers and Blin-Stoyle (1952a), we express each product of a heavy- 
particle Dirac operator and a spherical harmonic as a sum of irreducible tensor 
operators (Rose 1957). Each term in the sum can be expressed (the Wigner- 
Eckart theorem) as the product of a vector-addition coefficient and a second 
factor, called the reduced matrix element, which is independent of magnetic 
quantum numbers. The reduced matrix elements, though independent both 
of electron energy and of magnetic quantum numbers, are dependent upon 
the structure of the nucleus; because of our ignorance they are usually treated 
as arbitrary parameters. There remains, however, under the integral sign a 
factor of the form (r/X)'*", X being the Compton radian wave-length of the 
electron and / and /’ being effective orbital angular-momentum quantum 
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numbers. Since the de Broglie wave-lengths of the electron and of the neutrino 
are rarely as small as 104, and since 4 is very much larger than any nuclear 
radius, the most important terms are those with the smallest values of /+/’. 
The degree-of-forbiddenness approximation of 8-decay theory is the rejection 
of all terms in the expansion of each type of matrix element except those 
having the minimum value of /+/' consistent with conservation of angular 
momentum. Resulting from the steps which we have discussed is a general 
formula for the transition matrix element, and hence for the emitted electron 
wave-function, containing coupling constants, reduced nuclear matrix ele- 
ments, energy-dependent lepton radial wave-functions, and assorted angular- 
momentum coupling coefficients; the angles specifying the direction of emission 
of the neutrino enter as the arguments of a spherical harmonic. 

The angular distribution of electrons is obtained by calculating the outward 
radial component of the particle current using the asymptotic wave-function 
discussed in the preceding paragraph. The angles describing the direction of 
neutrino emission remain as parameters. Of course contributions from dif- 
ferent orientations of the resultant nucleus and of the neutrino spin must be 
summed, and the appropriate average over orientations of the initial nucleus 
must be performed. The polarization of the electrons is found by calculating 
the statistical-mechanical average of the expectation value of the operator 
é in the electron rest-frame. Taking advantage of the similarity in form of 
the expressions for the angular distribution and for the spherical components 
of the polarization vector we introduce functions S“ such that So° gives the 
angular distribution and S,“ (with « = +1,0) gives the polarization. The 
suMs over magnetic quantum numbers in the expression for S/ are carried 
out with the aid of the properties of various angular momentum recoupling 
coefficients. In the resulting formula, (48), for S# the angular co-ordinates 
of the directions of emission of electrons and neutrinos appear as arguments 
of spherical harmonics of orders g and g’ respectively. Information about the 
orientation of the source is contained in the factor Z(k) (see (49)), which is 
essentially the ‘‘degree of orientation of order k’’ discussed by de Groot and 
Tolhoek (see Siegbahn 1955). Except for trivial factors the remainder of the 
quantities upon which S/ depends have been collected into a set of coefficients 
R(gg'kk’t). In the general formula (58), for R(gg’kk’t) we have grouped the 
products of electron radial wave-functions into certain convenient com- 
binations which we call L,,(¢or), etc.; these functions depend upon the 
electron energy and the nuclear charge as well as upon the nonmagnetic 
quantum numbers of two electron waves; they are, of course, generalizations 
of the L,, etc. defined by Greuling (1941). Simple formulae are derived for 
the generalized functions for the approximation in which the nuclear radius 
is negligibly small compared to the electron de Broglie wave-length. 

The remainder of the paper is devoted to expressing the results corres- 
ponding to certain types of experiment in terms of the coefficients R(gg’kk’t), 
and to evaluating these coefficients for allowed and first-forbidden transitions. 
In Section 4 the formulae for the shapes of spectra of unoriented sources are 
examined, and the matrix elements of our spherical tensors are related to 
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those of the traditional Cartesian tensors. The formulae for longitudinal 
polarization of electrons emitted by such a source are shown in Section 5 to 
be very similar in form to the shape-factors of Section 4. We show how to 
calculate the polarization by means of a simple set of modifications of pub- 
lished shape-factors. The angular distribution of electrons emitted by an 
oriented source is given by the formula (83) of Section 6; expressions (87) 
and (88) can be used to calculate the polarization. The angular correlation of 
the directions of emission of electrons and neutrinos is to be found in (94) 
of Section 7; electron polarizations can be calculated from (95) and (96). 
Explicit formulae for many of the pertinent coefficients R(gg’kk’t) for allowed 
and first-forbidden transitions are to be found in Sections 8 and 9. 

Beta-decay is more complicated than many other somewhat similar prob- 
lems of nuclear physics, because two particles are emitted simultaneously, 
and because both these particles must be treated relativistically. The emission 
of a pair of particles introduces no great difficulty in principle, though it does 
increase the complexity of the results considerably. One of the more noticeable 
effects of the two sources of difficulty is the large number of different types 
of nuclear matrix elements which must be handled. 

In treating the motion of an electron in a general, central, electrostatic field 
it seems to be necessary, or at least more convenient, to introduce a particular 
representation of the Dirac operators. We shall use a notation in which the 
Dirac equation for the four-component wave-function, ¥, of an electron of 
energy E is 


(1) [—a-p—B+V(r)]¥ = Ey. 


V(r) is the potential energy of the electron in the electrostatic field of the 
nucleus. In (1) we have used the customary units in which the unit of mass 
is the electron mass, m, the unit of length is im~'c~', the unit of momentum is 
mc, and the unit of energy is mc?. The traditional representation used in p- 
decay theory is that given by 


; ae _ fo €). 
(2) e=(( E «-(9 oi 


the 2 X 2 matrices J’ and @’ are the unit-matrix and the Pauli matrices—our 
6’ is identical with the ¢ of Schiff (1955, p. 232). If ys is defined by the equation 
Ys = —ila,aya,, then it can be shown that the operator for the spin of a Dirac 
particle is 5 6, with 6 = y;@ = ays. In the representation (2) we have 


Ga) 69) 
ae ae ~ <a ee 


The form of the Hamiltonian density for the interaction responsible for 6- 
decay, 4%, which we use is that given by Lee and Yang (1956, Appendix). 
Because we will need to discuss selection rules for transitions between nuclear 
states of definite parities, it is convenient to express Ws in terms of operators 
with simple behaviors under spatial reflections. In other words we wish to 
express the products of the matrices y, used by Lee and Yang in terms of J, 


(3) Y5 


pn RT eRe nc ae “TOA SEEN EBSD 
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8, 6, a, and v5. Let us suppose that the y-matrices are so defined that the time- 
dependent wave-equation for a free particle is 


ly-V —tys(0/dt)+J]y = 0. 


If V is replaced by 7p and 0/dt by —7E, this wave-equation can be trans- 
formed into the free-particle form of (1) if 


(4) 71 = —B, vi = —2Ba, for k = 1, 2, 3. 


Elimination of the y-matrices from the Hamiltonian density of Lee and 
Yang by means of the relations (4) leads to the formula (5) for #%. 


(5) Me = (W*B®)-[(Y*B(CstCs'y5)¢] 
+ (W*d) -[Y*(Cy+Cy'y5)o] — (U*a®)- [Y*a(Cy+Cy'ys)¢] 
+ (W*B 6) - [y*86(Crt+ Cr'y5)o]+ (Y*Ba®) - [y*Ba(Cr+Cr’ys)¢] 
+ (W* 6d) - (Y*6(C44+Ca’y5)b] — (WF ys) - [Y*y5(CatCa’ys)o] 
+ (W*By5®) -[W*Bys(Cp+ Cpr’ys)¢]. 


The symbols ¥, ©, ¥, and ¢ stand for the Dirac wave-functions for the proton, 
neutron, electron, and neutrino respectively; ¥* is the Hermitian adjoint of 
W. When the primed coupling constants are equated to zero, expression (5) 
reduces to the customary 6-decay interaction—see the chapter by Konopinski 
in Beta- and Gamma-Ray Spectroscopy (Siegbahn (1955)). The operator 
Ys = Yi1v2v37¥4 defined by Lee and Yang is identical with the operator ys 
= —la,a,a, which we have already defined; hence the coupling constants of 
the two forms of .#g are precisely equivalent. 

A notation which enables us to write (5) in a convenient, shorter form has 
been devised by Banerjee and Saha (1954). The eight Dirac operators appear- 
ing in (5) are put into correspondence with the eight different combinations 
of the values of three two-valued parameters £, 7, and ¢. The definitions of 
these parameters depend upon the properties of covariants of the form 
(y*Ky), K being one of the eight Dirac operators, Y being a solution of an 
equation of the form (1). The variable é is 0 or 1 as the covariant is a scalar 
or a vector; 7 is 0 or | as this covariant is invariant or changes sign under a 
spatial reflection; ¢ is 0 or 1 as the covariant is invariant or changes sign 
under a total reflection. The operators for which 7 = 0 are those sometimes 
described as ‘‘nonrelativistic’’, viz. J, 8, 6, and 66. Those operators in which 
8 appears explicitly have ¢ = 0. The parameters 7 and ¢ also specify the com- 
mutation properties of the operator A with 6 and ys respectively, i.e. 


(6) BK = (-—1)"KB 
and 
(7) 5K = (—1)*'Kys. 


It is worth noting that the operator formed by multiplying the operator with 
parameters (, 7, ¢) on the right with ys has parameters (, n+1, ¢), » = 2 
being here defined as equivalent to 7 = 0. 
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In order to write (5) in the desired simplified form we need the spherical 
components 14 of a three-vector A; according to the usual definition 


(8) Al = —2-4(4,+iA,), A°=A,, A! =2-4(4,-iA,). 


The spherical components of r, for example, are (42/3)!rY\4. With this 
notation the scalar product of two vectors A and B becomes 


A-B = >> (-—1)*4“*B". 


It is also convenient to extend the definition of spherical components to a 
scalar, A, with A+! = 0, A® = A. We are now able to write the interaction 
Hamiltonian density in the form 


8 1 
(9) He = Ye (-1)" De (1) "(HK *#)[WKn' (Cot Cr'r5) 4. 
n= _ 
The coupling constants C, and C,’ are equal to the Cy and Cy’ for the inter- 
action to which K,, belongs. The parameters £, 7, and ¢ are assumed to be func- 
tions of ». In what follows no significance is attached to the exact corres- 
pondence between the K,, and the eight integers spanned by 2. 

Observable quantities connected with single 8-decay depend upon the 
coupling constants only through certain quadratic combinations. To facilitate 
our discussion in later sections let us introduce abbreviations for these com- 
binations now. 


(10a) Dxy= Dyx = Re(CyCy*+C,y’Cy"*), 
(10d) Axy = —Ayx = Im(CyCy*+Cy’Cy"*), 
(10c) Gxy = Gyx = Re(Cx’Cy*+CxCy"*), 
(10d) Tey = —lyx = Im(Cy’Cy*+CxCy"*). 


If Pauli’s (1957) more general interaction is to be used, then our formulae 
are to be changed by replacing Dyy and Ayy by the real and imaginary parts 
of Pauli’s invariant Kyy, and by replacing Gyy and Ixy by the real and 
imaginary parts respectively of the invariant Lyy (Liiders 1958). 

The practicability of working out the consequences of a general .#3 of the 
form (9) for arbitrarily forbidden transitions stems from the discovery of a 
simple analytical formula for the elements of the matrices A,“. Following 
Spiers and Blin-Stoyle (1952a) we write each of the eight matrices appearing 
in (5) as the Kronecker product, as defined by Banerjee and Saha, of two 2 X 2 
matrices. By inspection one learns that each of the scalars is of the form 
6 X I’ and that each of the vectors is of the form 6 X ¢’, 6 being a 2 X 2 
matrix whose elements are functions of 7 and ¢. The four components of a 
Dirac wave-function are to be labelled by two two-valued indices, 8 and o. 
The two upper elements have 6 = 1, the lower two 8 = —1; the upper element 
of each pair has o = 3, the lower o = —}3. Similarly we write A,,3. for an 
element of a 4 X 4 matrix A. In the factorization into Kronecker products 










: 
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which we have been discussing, the elements of 6 are designated by the indices 
8 and 8’, those of the matrices J’ and @’ by the indices o and oa’. In fact 


f+1 
bag’ = 68 6g, (1-29) "+ 
Next note that the matrix whose oo’-element is 
(—1)**' 24(3 $0 —olé w) 


reduces to J’ when — = O and to o# when £ = 1. We use the notation 
(j:jzmim2|jm) for the vector-addition coefficient (j:j2mmz2|j1j2jm) of Condon 
and Shortley (1953). Thus we have 


(11) (Kn") 308"0 7 BE+155 (1299 (— 1)*-#' 23 (44 o —o'|Eu). 


Note that we have been able to write the dependence on ¢,o’ in a form con- 
venient for carrying out sums over magnetic quantum numbers. 

The representation defined by (2) and the sign conventions of (1) are 
believed to be those most commonly used in 6-decay. The author wishes to 
express his thanks to Drs. Konopinski, Greuling, and Pursey for confirmatory 
correspondence on this point.: Though Spiers and Blin-Stoyle do not specify 
a representation explicitly, one can deduce from their Table 2 that the repre- 
sentation used is equivalent to ours with the signs of @ and ys changed; there 
appears to be a misprint of sign in this table in the case of the elements of Ba. 
Likewise one can deduce from their expression (29) that Banerjee and Saha 
used the same representation as Spiers and Blin-Stoyle. Since it is easily 
shown, using theorems given by Pauli (1936), that quantities of the form 
(¥*Kq@) are invariant under change of representation, the expressions for 
the various observable quantities worked out with the different representations 
must be precisely the same, even though the nuclear matrix elements appear 
to depend upon the K,. 

So far, we have discussed negative-electron decay only. The Hamiltonian 
density for positron decay is of course the complex conjugate of (5), and may 
be put into a form analogous to (5) by a transformation of the coupling 
constants. Let Cy+ and Cx’*+ be coupling constants in the positron-emitting 
part of interaction after it has been brought into the form (5). It is easily 


shown that the new coupling constants must satisfy Cy’+/Cx+ = —Cy’*/Cx*. 
The definition of the transformation may be completed by putting Cx* 
= Cx* when X stands for V or T and Cyt = —Cx* when X stands for S, A, 


or P. It is also necessary to take the appropriate sign for V(r) in equation (1), 
which defines y. In the remainder of this paper we shall deal with negatron 
decay only; results for positron decay can be obtained by making a suitable 
transformation of the coupling constants and by changing the sign of Z. 


2. ELECTRON WAVE-FUNCTION 
In this section we calculate the probability amplitude for electron position 
at a large distance from the transforming nucleus. As is customary the absorp- 
tion of an antineutrino is considered instead of the emission of a neutrino; 
plane-wave eigenfunctions are used for the antineutrino. The nucleus is 
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supposed to make a transition between states with definite magnetic quantum 
numbers. Thus the calculated electron amplitude is a function of the nuclear 
quantum numbers as well as of the momentum and spin of the antineutrino. 

It is convenient to follow the simple treatment of Greuling and Meeks 
(1951), who use ordinary time-dependent perturbation theory, in order to 
find the amplitude of the emitted electron wave. The field-theoretic Hamil- 
tonian density of Section 1 is replaced by the Hamiltonian Hg; in place of 
expression (9) we now consider 


(12) Hs = it 2. (—1)"**K, *(H)K,"(L) (Cat Cy'¥s)tu715(Xu—Xz). 
n » 


The indices H and L refer to heavy and to light particles; r_ and ry are oper- 
ators transforming neutrons and antineutrinos into protons and electrons 
respectively; Xq and X, are position vectors for the heavy and light particles. 
In writing (12) a summation over all neutrons in the nucleus has been omitted 
for the sake of simplicity. 

The initial wave-function of the system is the product of a nuclear wave- 
function, ®, and the antineutrino wave-function. Wave-functions of a free 
particle in our representation are well known; see Schiff (1955, page 327). 
By inspection of these expressions one can show without difficulty that the 
components of the wave-function of a particle of zero mass, of energy —4q, 
and of momentum —q can be written as 


sa O3,' = (2x)'(e  s—o ol} s)¥s(6", oe, 


with 


B. 


tole 
Nix 


ée= 


Two orthogonal solutions are specified by the two values, +}, of the para- 
meter s; the normalization ¢*¢ = 1 has been used. For a particle of nonzero 
mass, s is an eigenvalue of a, in the nonrelativistic limit; since we shall always 
sum over all values of s, we need not enquire further about its physical sig- 
nificance. The angles 6’ and ¢’ are the spherical polar co-ordinates of q, the 
neutrino momentum vector, with respect to the z-axis; Y,“ is a normalized 
spherical harmonic equivalent to that defined by Condon and Shortley. 

The wave-function of the final state of the system is the product of a nuclear 
wave-function and an electron wave-function, ¥. We expand y in a complete 
set of wave-functions, ¥/"*. The quantum numbers j and m refer to angular 
momentum and its z-component; the two types of solution obtained from 
the two values, +1, of a@ are analogous to the two solutions with different 
orbital angular momentum for a given value of 7 in a nonrelativistic theory. 
If the relativistic parity operator is defined as ~8 X (a spatial inversion), then 
its eigenvalues are (—1)**4, 


(14) je = r'FY(p, r)(L4 m—o o|j m)Y7""(6, o), 
with 
l = j+4ap. 
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The wave-function (14) is identical with that given by Rose (1955), who uses 
the same representation, if the following relations between the radial parts 
are satisfied: 


(15) r'Fy = —if., r'F,* = Se x = —a(j+}). 


Differential equations satisfied by f, and g, are given by Rose. The quantity 
p is the magnitude of the linear momentum at a large distance from the 
nucleus, and is related to the total energy, E, by p = (E?—1)?. The normaliza- 
tion of the radial wave-functions is fixed by insisting that they have the 
asymptotic expansions 


(16a) f,e = —r(1—E-")! sin (pr+é), 
r(1+E-!)! cos(pr+s). 


(16d) &« 


The normalization (16) is usually used in 6-decay (Konopinski and Uhlen- 
beck 1940; Rose et a/. 1953); it gives rise to the relation 


(17) sume, ry?” (EE v)dr = rpE"5(E—E’)6;)5mm'Saa’- 


The phase-shift 6 depends upon « (i.e. upon 7 and a); if V(r) falls off as r~, 
6 also has a logarithmic dependence upon r (Rose 1937). If the potential energy, 
V(r), is nonsingular at r = 0, then the factor F%"* of the radial wave-function 
approaches a finite, nonzero value as r approaches zero. For the values of E 
usually encountered in 8-decay, F?* is practically constant throughout the 
nucleus. Our relations (15) differ in sign from the corresponding equations of 
Spiers and Blin-Stoyle and of Banerjee and Saha, because these authors use 
a different representation of a@. 

The application of the time-dependent perturbation theory is straight- 
forward except for the apparent presence of ingoing as well as outgoing waves 
in the asymptotic expansion of the electron amplitude. Using a method 
originated by Bethe (1930), Greuling and Meeks showed that the ingoing 
waves remove themselves by destructive interference; note, however, that 
because of an apparent misprint of sign their expression (9) appears to lead 
to a cancellation of outgoing rather than of ingoing waves. Dropping the 
time-dependent factors, the operator tq, and a unimodular factor, we may 
write the asymptotic expansion of the outgoing wave as 


(18) vee ~ [E(E—B)]*(pr)7'e”” Se? f|Hg|i) (Am —oo| jm) ¥"~“ (6, ¢), 


jma 


with 
A = j+ 3a8, 


(19) (f\Hpi) — = +. (— 1 yam \ (w*K,, *b)[yrr*K,"(Cat+Ca'ys)o" lar. 


6(j, a) is the phase-shift, as in (16), associated with ¥’"*. Relation (18) gives 
a formula for calculating a component, gs, of the outgoing electron wave. 
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The various results to be obtained in the rest of the paper all flow from the 
use of the formula (18). But before we can proceed much further, we must, 
of course, evaluate the matrix element defined in (19). 

Expression (19) consists of a sum of matrix elements, between two nuclear 
states, of operators of the form Q, 


Q = (—1)*K,*(y""* Ko"). 


In order to relate the quantum numbers of the electron—neutrino emission 
to the nuclear angular momentum quantum numbers, we next express Q as a 
sum of spherical tensor operators. The lepton covariant in Q is a function of 
r; it must therefore be expressed as a sum of spherical harmonics. The plane 
wave in the neutrino wave-function is expanded in terms of spherical waves 
in the usual way (Blatt and Weisskopf 1952, Appendix A). It is convenient 
to define a function G,(gr) = (qgr)~‘j.(gr), j:.(z) being the spherical Bessel 
function; G,(gr) is very nearly constant throughout the nucleus. Using the 
analytic formulae for wave-functions and for the K, one may write 


Q = 84°"K,* > BP '83.0—2ne" De (—i)"g'r*" Fr (p, r)Gy (qr) 


me ee (—1)""*"'(1 3 m—o alj m)(e 3 s—o’ a’ \4s)($ 40 —o' lw) 


x ¥i™""(G, 6) Ve" (0, 6) Ve" (6, 6) Vu * 6’, ¢'). 


It is well known (e.g., see Blatt and Weisskopf, Appendix A) that products 
of the form Y,”"(6,¢) VY" (9,¢) can be expressed as linear combinations of 
Yro'(6,¢) with L ranging from l'—/ to /+l'. We next express the product 
r’ VY ¥>*K,~* asa linear combination of the tensor operators 7), (J, L), defined 
by (20). 


(20) 7,%(J,L) = (4x)'i*U"'(J/,L,8) D) (LE M+y —plJ M)r* VR, *, 


U(J,L,&) = (J(L+é) YE (LY) [2L4+1)!2°-7}' for J > 0. 
(24 U(000) = U(O11) = 1. 


Recall that € is a function of n. The presence of the factor 7” in the definition 
20) leads to a convenient behavior of reduced matrix elements under time 
reversal, as will be seen later; U~'(J, ZL, &) is included to bring about an 
equality between certain products of our reduced matrix elements and the 
corresponding sums of products of components of Cartesian tensors found 

the work of Greuling and of Pursey. It should now be clear that the left 


side of (19) can be written as a linear combination of terms of the form 
. 
on Ja ‘ l¢+U—L - M 
(22 jv Fy* Gyr’ 7,” (J, L)@ae. 
« 


Jn order to separate the energy-dependent factors from the nuclear matrix 
elements we apply the mean-value theorem to (22); i.e. we take /)“*Gy out- 
side the integration and evaluate it at p, an appropriate mean-value of r 


roughly equal to the nuclear radius, Because /%“*G, does not vary much 
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within the nucleus (in other words, because the electron and neutrino 
de Broglie wave-lengths are much larger than the nuclear radius), it is not 
necessary to use the true value of p. With electron wave-functions pertaining 
to a point-charge it has been customary to put p equal to the nuclear radius; 
with the wave-functions for a distributed charge it would probably be better 
to use half the nuclear radius. Although one would expect, in general, to find 
a different value of p for each different matrix element, we shall follow the 
usual approximation of using the same value for all. 

Next we apply the Wigner-Eckart theorem (Rose 1957) to the modified 
form of (22). The quantity r'+"-* 7,™(J, L) is, of course, a tensor operator 
of order J. Let J;,.; and J;,M, be the quantum numbers of the nuclear 
spin and of its z-component in the initial and final nuclear states. Then the 
Wigner-Eckart theorem says that, if we write 


(23) (flr*Ty-™ (J, L)\i) = (Kat JL; k)(—D 7+" IM,M IM), 


the reduced matrix element (K,;JL;k) is independent of the magnetic 
quantum numbers; the notation for the reduced matrix elements employed 
is meaningful only for tensor operators which differ from those defined in 
(20) by at most a power of r. This definition of reduced matrix element differs 
from those of Rose (1957) and of Sharp (1957) by trivial factors. 

After the right side of (19) has been expressed in terms of reduced matrix 
elements, all summations over the magnetic quantum numbers wy’, ¢, and o’ 
can be carried out with the aid of recoupling theorems. The properties of 
the Racah coefficient W(adcd, ef) are now well known (Biedenharn e¢ al. 
1952), and will be used freely in the next sections; relevant properties of the 
X-coefficient are summarized in the Appendix. By the operations we have 
described the matrix element of (19) can be brought into the form 


(f\Hali) = 4(—1) DE (- 1) ED B81 C8s, 293° + OSs, ens] 


x YE (JIM MIM) (Ka; JL; + U—L) US, £8) 


JLU 


24) 4 i »: Shee enn Xe 
\ x > (4h sm—M-—s!|j’m—M)(j i’ m M—m') J MY," (0, @ 


i woe 


x i O'R 1+ ( pers a, 


— 
wer 
™ 


l= j+}a8, M=M,—-M, ¢ = 4-38 


a ope 


( ay ee 


is a modified V-coefhicient and is detined in the Appendix: it is zero unless 
[+1'+L is even. Because A,“ys is obtained trom A,* by changing the value 


of yn, the coethcient of C,o in (24) is the same as the coethcient of C, except tor 
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a change in the value of n. In deriving (24) use was made of the identity (25) 
containing a Racah coefficient with special values for some of the parameters; 
a proof is given in the Appendix. 


(25) [2(2e+1) (21’+1)]*(el’00|h0) W (e3/'7’, 3h) = 3(1+(—1)*"*"]. 


Conservation laws can be used to eliminate reduced matrix elements speci- 
fied by certain combinations of J, Z, &, and n from the expansion (24). Note 
that J is the quantum number of the total angular momentum carried off 
by the electron—neutrino pair; L and & represent the corresponding orbital 
and spin angular momenta. As may be seen from the first vector-addition 
coefficient in (24), conservation of angular momentum requires that 


we shall call AJ the change in nuclear spin. Let the nuclear parity change 
Am be defined to be +1 when the initial and final nuclear states have the 
same parity and —1 when they have opposite parity. Then the nuclear matrix 
element is zero unless 


(27) An = (—1)"*", 


The reduced matrix elements in the expansion (24) of (f|H, 7) contain a 
factor r'+"; / and I’ play a role here analogous to that of the orbital angular 
momenta ina nonrelativistic theory. Because the value of r (in units of hm~'c—!) 
is of the order of magnitude of 0.02, the dominant terms in (24) are those with 
the smallest value of /+/’. The approximation usually made in 6-decay theory 
is that of discarding all terms in the expansion of (f| Hg|7) except those with 
the minimum value of the degree-of-forbiddenness parameter, f = /+/'+n, 
permitted by the nuclear quantum numbers. The term 7 is included in the 
definition of f because matrix elements containing a K," for which 7 = 1 
are expected to be smaller than those for » = 0 by a factor roughly equal to 
some mean value of v/c for the nuclear states in question. With this definition, 
matrix elements satisfying the nuclear selection-rules (26) and (27) are of 
approximately the same magnitude if they have the same value of f. Note 
that there are no allowed, i.e. f = 0, matrix elements with 7 = 1. The minimum 
value of /+/’ is L. In discussing the dependence of ZL upon nuclear parameters 
it is convenient to divide the nuclear transitions into two classes dependent 
upon the relation between AJ and Ar. 







I. If Ar = (—1)4/, then the minimum value of f is AJ, and three types 
of matrix element must be considered. 











(i) » = 0, f= Dior 1, aay J eA 
(11) =I, =a b= J—i, SD sep: 
(iii) » = O, fia bie 1. J = AJ+1. 






When AJ = 0, case (ii) and the — = 1 alternative of case (i) are to be 
omitted; case (iii) does not satisfy (26) if either J; or J; is 0. 






I]. If Aw = (—1)4/*!, then the minimum value of f is AJ—1, for AJ # 0, 





SMELT CD TI AIT OREN REE RRO 


Set PEPE NET HO SAPNA EB aera ES 


SPEAR RRR NRE RIN RENEE 


PRT EP REECE EstENERRETENTEN 


NA TTT NSN TC 
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and only terms with 7 = 0 and é = 1 contribute. These matrix elements 
are the same as those considered in case (iii) of I. If J; = J; = 0, then 
the only types of matrix element which can occur are: 


(i) » =0, g = 1, L=1, J =0; 
(ii) » = 1, t= 0, L=0, J=0. 


Both these types are first-forbidden, i.e. have f = 1. If J; = J; ¥ 0, 
matrix elements of types (i) to (iii) of I also occur, with AJ put equal to 1. 


If the 8-interaction is completely lacking in either A or 7, then the shapes 
of spectra of Class II are completely calculable—i.e. they do not depend 
upon unknown ratios of coupling constants or matrix elements. Such transitions 
are often said to be unique. We shall extend the definition slightly and term 
the transitions of II the U class; similarly transitions of I will be called the 
non-U class. The degree of forbiddenness of a given transition is defined to 
be the minimum value of f corresponding to a matrix element in the relevant 
expansion of the form (24). 

In treating the effect of .the P-interaction in first-forbidden (AJ = 0) 
transitions we have not followed those authors who have assumed that 
{¥*Bysédr = (875; 00) is so small that the contributions from P are negligible 
for any reasonable value of Cp. A brief discussion of this question with refer- 
ences to the literature is given by Konopinski (Siegbahn 1955). 

Although it is only in the most careful measurements of spectral shapes 
that the presence of matrix elements other than those of minimum forbidden- 
ness could be detected, we shall include all terms in our general formulae 
of the next section. With the formalism employed this inclusion requires no 
extra effort. Only dominant terms will be given, however, in the particular 
formulae discussed in the remaining sections. The degree of forbiddenness 
of the most important higher terms is always two greater than that of the 
dominant terms. Formulae for the second-degree corrections to allowed 
spectra have been given by Zweifel (1954). 

Let us assume with Longmire and Messiah (1951) that the action of the 
Wigner time-reversal operator, 4 on nuclear wave-functions is described by 
(28), 

(28) NGF = Os, 
the phase ¢ being dependent on the nuclear state. .¥ is the product of a 
unitary matrix, U, and an operator which takes complex conjugates. In our 
representation we may put LU’ = toy, or 

Uses'e = 539'(—1)"*e,-2", 
whence it follows that 
(29) UIT, MUR Ut =: (=1)""" TF Ue). 


(The elements of the matrix [A*]” are the complex conjugates of the elements 
of A.) Using (29) one can show without difficulty that 


(30) (Ki; Jil; ki)*(Ko; Jobo: ko) = (Ki; Jil; ki)(Ko; JoLo; k2)*. 
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3. GENERAL ANGULAR CORRELATIONS 

The angular distribution of electrons emitted when the absorption of an 
antineutrino of given momentum and spin causes a particular nuclear transi- 
tion may be obtained by calculating the particle-current for the wave-function 
(18). For y’s satisfying the wave-equation (1) the current is — (~*ay); the 
angular distribution function for the electrons is found by multiplying the 
radial component of this current by r?. To obtain results comparable with 
experiment we must sum over all possible values of .W, (because of the ortho- 
gonality of final nuclear wave-functions, cross terms between y's corresponding 
to different values of 1, do not occur) and average over all possible values of 
M,. We shall write P(M;,) for the probability that the nucleus is initially in 
a state with magnetic quantum number /,, the normalization condition 


> P(M,) =1 

M; 
being understood; the form of the dependence of P upon M; is determined by 
the mode of producing the oriented nuclei. Summations over s and over the 
states in an appropriate part of the neutrino momentum space are also to be 
performed. Let the probability that an electron is emitted into an element 
of solid angle dQ at (6, ) with its energy in an interval dE at E while a 
neutrino is emitted into dQ’ at (6’, ¢’) be W(6, 0’, ¢—¢’)dQdQ'dE. Then 
(31) W(0, 0’, 6—¢') = — D> P(M,)(20)~q°r?(Y*ar). 

Mi Mys 

Note that the y in (31) depends implicitly upon W;, M,, and s, all quantum 
numbers related to orientation of angular momentum vectors. The symbol 
a, in (31) represents the component of the vector e in the outward radial 
direction. It is customary in the theory of 8-decay to quote results in the form 
of a correction factor, So(6, 6’, ¢—#’) as defined in (32), which relates the 
statistical distribution to the true distribution. 


(32) W(6, 0, 6—$') = (20°)-'Epg?Fo(E, Z)So(8, 0’, 6-9’). 

F\(E, Z) is the Fermi function as defined by Greuling (1942); extensive tables 
of a closely related function have been published by the National Bureau of 
Standards (Fano 1952). 


Simplification of the right side of (31) is facilitated by the use of a more 
convenient form of the current. Equation (1) may be rewritten in the form 


(33) vy = (V—E)by—Ba-py. 


By taking the arithmetic mean of the two formulae found by substituting 
expression (33) for each of the y’s in — (¥*ay) in turn, one can show that 


(34) — (Wray) = 2{[(py)*-a]bay+y*a8(a-p)y}. 
Now to the first order in r~' the y of (18) satisfies 


(35) py = pny; 
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here p is the linear momentum operator, n is a unit outward radial vector, 
and p is a number equal to (E?—1)?. Substitution of (35) in (34) with the 
use of the anticommutation relations for (6, a) leads to 


(36) — (Yay) = —p(y*By). 
Similar to the derivation of (34) is the proof that 
(37) vy = (V—E) (¥*6y) — 31 (Y*Ba- py) + (py)*- ay]. 


Since V(r) approaches zero as r goes to infinity, we can, with the aid of (35) 
and of the anticommutation relations, transform (37) into 


(38) y*y = —E(y*6y), valid for large values of r. 


The application of (38) to (36) gives rise to the asymptotic formula (39) for 
the outward radial current. 


(39) — (Yaw) = pE“"(y*y). 


Since in the units being used pE~ is equal to the magnitude of the velocity 
of the emitted electrons, expression (39) simply says that the asymptotic 
current is the product of velocity and particle density. It is of interest to 
note that the relation (35), and hence the equations (36) and (38), are exactly 
true for the momentum eigenfunctions of a free particle (Bethe 1948). Com- 
bination of (31), (82), and (39) leads to the formula (40) for the correction- 
factor. 

(40) So(0, 0’, 6— ¢") = [4E°F(E, Z)]* SS P(M)r’ (Wy). 

ss 


M;M 


The question of a suitable definition of a vector polarization, {, for a Dirac 
electron has been discussed by Tolhoek (1956), who has pointed out that the 
direction of { is that of y*(J—8) éy. If, as is desired, the magnitude of { is 
to be equal to the degree of polarization, then { should be defined to be the 
appropriate statistical-mechanical average of the expectation value of the 
operator (1+£~-!)-!(I—8) 6. Thus 
(41) ¢ Do P(Mi)(¥*¥) = Do P(M)A+E")'l¥* (I-8) oy). 

MiMyss MiMyss 
Because the efficiencies of the various methods for measuring polarizations 
described by Tolhoek are given in terms of %, this vector provides a con- 
venient way of expressing polarization results. 

In deriving formulae for %, it is more convenient, of course, to work with 
the spherical components, ¢#. Let us define a set of four quantities S$“; 4 = 0 
when ¢ = Oand wp = +1 or 0 when? = 1. 

(42) Si = [4E*Fo(E, Z))" DY P(M dr’ (Gy). 


Mi; Mss 
For ¢ = 0 we put @ = J, and So° of (42) is the correction-factor So of (40). 
For t= 1 we put @ = (1+£7')-'(I—8)o*; the spherical components of 
the vector polarization are then given by (43). 


(43) ¢= Si4/So. 








1214 CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 


Thus all the required information about the distributions can be obtained 
from a knowledge of the quantities S/. 

Next we turn our attention to the evaluation of (¥*@/W), using the y of 
(18). By a slight adjustment of formula (11) for the matrix elements of K,*, 
one can show that 


1 


(Of) poser = (1 +6 E-!)—8g9- (1 —6 18) (— 1)-223(3 3. —o'\t yp). 
In carrying out the required contraction, sums of the form (44) are encountered. 
(4) B= 2) (1-BE™')(1—duB) 2 (—1)" "4 be —o"lt w) 
xX Ar 3 m\—o olf my) (Xz 3 mM2—a' a" |j2 m2) Yxi—* (6, ¢) gael (0, o). 
hi = jit 3a 6. 


By expanding the product of the two spherical harmonics and by recoupling 
the two vector-addition coefficients, expression (44) can be transformed into 
(45). 


(45) B= D> [(4e)(Qji +1) (Qje+1) (2t4+-1) 2k’ +1) (2g+1) "IF? 
k’.g 
x (jije —m m2|k’ mMs.—m,)(k' t M2—My, bg M2—m+ph) Peis | ¢) 


> (—1)2*™™ (1 — BE~*) (1-8 nB)[(241 +1) (242+1)]' 


jie M 
(A14200|g0)X| j2 2 Ae 
B. ie 


Because \; and )2 are functions of 8, it would appear at first sight that the 
summation over 8 in (45) could not be carried out explicitly. However, we 
show in the Appendix that 


ja & Ms 
(46)  [(2d1 +1) (2A2+1)]'(A1200|g0)X| jo 4 De 
k’ t g 


= (= 1) 2° (1) (2-41) (2k +L) (hr fe 2 — 41’ 0) 
x {{~ 1)??? (12/0020) — 5 1A joy, na (R's g)] 


with 
(47) Ayo(R’, g) = 2A (A+ 1] (1k — 1/90) [(u1) + (- 1) +1). 
Now, recalling that \; = j7;+4a,8, one sees that 
(—1)7-"* = — fas 
and that 


(gy a (1)? “Fees, 








EEE 


ETON 


a 


a 


a 
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Substituting (46) in (45) and using 
A fem 


we can now carry out the summation over 8 in the last line of (45) without 
difficulty. The presence of the factor 


[1+ (—1)**2-"a,a,] 


shows that the order of the spherical harmonic describing the electron emission 
is even or odd as the parities of the electron states 


” aii and e 
are the same or opposite; this fact is used extensively in the manipulations 
which follow. 

Once the spinor contraction has been performed one can bring S/ of (42) 
into the form (48) using the definitions of the two recoupling coefficients, the 
orthonormality relations for the vector-addition coefficients, and the identity 
(A.15) of the Appendix. 

(48) S*£= YS (4n)7(—1)"[(2t+1) (2g’+1)] 22 (k) R(gg’kk’t) 


kk’ gg’ 


xX dS (kg’Om|k'm) (k’ gm w—mlt uw) ¥,"(0’, ’) ¥,""(0, ¢). 


™ 
The angles (6’, ¢’) and (@, @) are, it may be recalled, the spherical polar co- 
ordinates with respect to the axis of nuclear orientation of the directions of 
emission of the neutrino and electron respectively. Information about the 
degree of nuclear orientation is contained in the functions Z(k) defined by (49). 


(49) Zk) = dD (-1)72I,4:1)'(Si J, Mi — Mik 0)P(M)). 


Mj 


It follows from the properties of the vector-addition coefficients that Z(k) 
is zero if k is greater than 2/;. Note that 


(50) Z(0) =1 and Z(1) = 3{J(J,4+1]}° Viz), 


(J iz) being the mean value of the component of nuclear spin in the direction 
of the axis of orientation in the initial state. If the decaying nuclei are aligned 
but not polarized, i.e. if P(—.M,) = P(M,), then Z(k) is zero for odd values 
of k. When the initial nuclei are unoriented, i.e. when P(.W,;) = (2J;4+1)7', 
then (49) reduces to 


(51) Z(k) = dxo. 


A discussion of earlier uses of coefficients similar to Z(R) is given by de Groot 
and Tolhoek in Beta- and Gamma-Ray Spectroscopy (Siegbahn 1955). 

The outstanding coefficient in (48), R(gg’kk’t), contains combinations of 
coupling constants, reduced nuclear matrix elements, electron phase-shifts, 
electron and neutrino radial wave-functions evaluated at r = p, and various 
functions of “angular momentum” quantum numbers. 
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(52) R(gg’kk't) = (—1)°*7/-* (8 Fop®)"[ (2: +1) (2k +1) (2k’+1)7'}? 
z, (—1 ymbitaatet sites pee ee 
nin2 Ji J2Li L211’ lr’ 


USL )U GLOW SAAS) FF (—* 


Jij2jr’ jo’ 
[(2jr +1) (2jo#1) Qj’ +1) jo’ +1) ]*Gr jo 4 — $1’ 0) (jr'je’s — 43g’ 0) 
jh jv Ji ane 
A je ae Se > ghee Knog 4, (gt Mga] 
k’ g’ k aia2 
[(tk’00|g0) —a25 1A joy, ny (R’, 8)] p> Bi ee rT re 
152 


ly 2 ji Le 3 je 
Adh’ 4 ji’ Aah’ 3 it 
Ly £1 Ji L. £5 J» 


ees: JiLy; h+d;’—L,)* (X.? 
Jala; latl’—La) Do By'Bs'[1+81'Ba(—1) 8] 


* rk , 
(Ca §61,—2n81' + Cn 831, (2m—181' [C25 ge, (1—202)82' + Cre 832, (2n2—1)82'] 


By inspection of the last lines of (52) one can see that the two parity con- 
serving terms (Lee and Yang 1956) C,,*C,, and C,,’*C,,’ have the same co- 
efficient. Similarly the coefficients of the two parity nonconserving terms 
C,,*C,,' and C,,’*C,, are the same. In the coefficient of the product of two 
coupling constants of which v are primed the nonvanishing of the Kronecker 
deltas demand that 


Bi'Bo’ = (—1)"*"*"B Bo. 
Now from the nuclear parity selection-rule (27) we can show that 
(—1)"*™ as (jy) 
and from the properties of the A o-coefficient we recall that 
(-1)"*" = (-1)", i= 1,2. 
Thus 


sat(-1" tol (—1)"t**8*49,8, 
ee (—1)""" +h-ntemi-yitg 2, 


(— 1 ) vto'+ii—J2 1d. = (— ays, 


Hence we have shown, as expected, that for the parity-conserving terms, 
v = 0 or 2, the sum of the orders of the spherical harmonics describing the 
electrons and the neutrinos must be even; for the parity nonconserving terms, 
v = 1, this sum must be odd. 








E 
: 
j 
f 
i 
; 
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To simplify the rewriting of (52) and to conform to custom, we replace the 
summations over j; and j2 by summations over the non-negative integral 
parameters up = j2— 3 and vy = j,—}. We have shown in the preceding para- 
graph how the summations over 8,’ and 2’ in (52) may be performed. The 
summations over di, @2, 81, and 82 result in the appearance of generalizations 
of the functions L,, M,, and .V, defined by Greuling (1942). The way in which 
this step is carried out depends upon the relative sign of a; and dz; it is there- 
fore convenient to eliminate the sum over a, in favor of one over a new para- 
meter o, also taking the values +1, defined by a; = oa». In order to arrange 
that the combinations of electron wave-functions that we are about to define 
be real numbers we make the expansion 


etss.02)-851-60] _ 4718.9 COS Qy» (G2, 7) +5,15iN Q,» (ae, o)] 


r=0 
wherein 


(53) Qu(a, 0) = 5(u+4, a) —d(v+}, oa). 


With each of the four choices 8; = +a;, By = +as the summation over a2 
in (52) gives one of the three types of function defined in (54) 


(54a) L,»(tor) = i" (2p? Fy) pe a F;* 10% ets va =n..(a, a), 
(54b) M,,(for) = i***(2p?Fo)* > (—a) FFAS * Fee i" S0,,(a, o), 


(54c) Ny(tor) =i * (29°F) Dd) al Fei’ Fo Sn, (a, @) 
a 

Note that each choice also corresponds to a definite relationship between 
J; and j; and between /: and j.. For example, when 6; = —a; and Bp: = —az 
we have /; = j1 — 3} and /, = j2—43, and summation over dz yields L,,(for); 
there is also a term containing L,,({+1, ¢7) when ¢ = 1. The index 7 in (54) 
takes values 0 and 1, the former corresponding to cos and the latter to sin 
in the $Y ambiguity. The parameter ¢ in (54) is an abbreviation for ¢:+¢2. 
Since ¢ appears only as an exponent to the quantity +a, and since the dummy 
a takes the values +1 only, it is clear that all even values of ¢ are equivalent 
to 0 and that all odd values are equivalent to 1. As has already been stated, 
a has the two values +1, corresponding to a; = +d. Thus the two possible 
values for each of the three parameters {, a, and 7 inserted in the L-, M-, 
and .V-functions of (54) give rise to 24 different functions of » and p (and of 
the nuclear charge distribution), each dependent upon the two integral 
indices uw and », in terms of which the general electron—neutrino polarization 
distribution can be expressed. For the Z- and M-functions different choices 
of u and v do not always lead to independent functions. From the definitions 
(15) and from the reality of f, and g,, it follows that 


(55) Fy** = —BF;". 
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Using (55) one can show without difficulty that 


(56a) Ly(gor) = (—1)'o**'L,, (or), 

(560) M,, (tor) = (—1)'o'*'M,,(¢o7). 

From (56) follow the useful results 

(57a) L,,(¢11) = 0 = M,,(¢11), for all ¢, 

and 

(576) L,,.(¢-—1r) = 0 = M,,(¢—11), if ¢+7 is even. 


Though V,,(¢11) = 0, no relation analogous to (576) holds for V,,(¢—1r). 
Equation (55) can also be used to show that each of the functions defined in 
(54) is real. The definitions (56) have been so devised that Z,,(010), 17,,(010), 
and .V,,(010) are identical with the Z,, W,, and .V, of Greuling (1942). Similar- 
ly L,,(110), 7,,(110), and V,,(110) are equal to the P,, Q,, and R, defined 
by Pursey (1951); these same three functions are called L,~, M,-, and NV,~ 
by Banerjee and Saha (1954); they occur only if the B-interaction is drawn 
from both S7’P and VA. 

Before rewriting R(gg’kR’t) in its final form we shall introduce one more 
simplification. Because m; and m, both run through the same eight values, 
the complex conjugate of the product of each pair of different matrix-elements 
occurring in (52) also occurs. Since, by (30), the complex conjugate of such 
a product is equal to the original product, we can see that each product of 
different matrix elements occurs twice. By re-naming summation dummies 
in such a way that the symbols bearing the subscripts | and 2 are interchanged, 
and by the use of the symmetry properties of the various angular-momentum 
coefficients, one can show that the two coefficients of a twice-occurring 
matrix-element product differ by the sign (—1)*+#’+7 and by the complex 
conjugation of the coupling constants. The incorporation into (52) of this 
simplification and of the changes described in the preceding paragraphs gives 
rise to (58). 


(58) R(gg’kk’t) = 4[(2k+1)(2k’+1) 7}? DS (1+6n) 


nigne 


YD (Ht 2I, 41)? WS Tet Jy RI) U(IiLibs)U (Sooke) 


JiLliJd2L2 


a" (CPC +C., 76,7) + (— 1" (CarOus* +Cn,'C,.” 11 (> l id 


Tt 


ve [a Cat Gn'Gn + (— 1)°(G,,'G,*+G,.G.™ oer fm I ery 
DS ot SY [r+ 1) ie +1) it +1) H+] 
o By fa’ 53°41" Ts’ 


di ht OT; 
(jij — 2/R'0) (jr jo’s — 3 \¢'0) Al jo je’ Je 
Be og 


/ ji— jr’ > Lo—le’—Lith’_ li’+l2’ + jo 
(—3yr"* jj le-k ns se *(1+(—1)° ji ” gf] 
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(Kno; Sele; w+ls—Le)(Kny; Jili; »+hi—L1)* 


wood js v 4 fi 
Ad 4) AD AH 
he te Bin. Sy dj 


X [(tk’00|g0) Ly, (for) —6,, Ay Ly (F+1 o 7)] 
+(K,,; Jobo; e+ 1 +13—L2)(Kny; Jil; yvt+1+i-—L, e 


~ 


utl & js ert th 
Ag 1. 3 je’ Ad ly > ji’ 
le tee Ly & Ji 


X [(tk’00|g0) MM, (for) —6;, Au My. (¢+1 ¢ 7)] 
+i(—1)**(Kys; JoL2; wt+1+ls’—Le)(Kni} Jily; v+/,'—L,)* 


etl 4 = je v > 2 
Ag 1,’ i je A, Ly’ 5 ji 
L, & Je hha I 


X [ (tk’00|g0) Nu. (fo7) —b Ay Nur (E+1er)] 
—i(—1) aK, ; Sole; ute! —L2)(Koi; JiLi; v+-14+h!—Li)* 


Le t je vt+l 2 ft 
Aj 1,’ } je Ao L,’ 2 ji 
bs E> J» I & Ji 


X [(tk’00|g0)N (for) —cb Ay N y(E+1er) |} 
S= mbitnletitsetmtmtinstéis wo = g’tktitr Ay = Ay(k’, g) 

c= bits je = ut+3 i = v+3 52 = Bn ino) 71 72911 29 pv) 11° ta! 
There are several tests that one may apply to (58). First, the expression 
must be invariant under interchange of all the quantities with subscript 1 
with the corresponding quantities with subscript 2, since the ordering of the 
operators AK, can have no physical significance. That such an invariance holds 
may be demonstrated with the use of (30) and (56) and of well-known 
properties of the various angular momentum coefficients. From the required 
reality of So° and of each Cartesian component of ¢ it follows that S“* = 

(—1)4S,;-*, which condition in turn requires that 


(59) R* (gg’kk’t) = (—1)*+9+9'+ *R(gg’kk’t). 


It is not difficult to show that the formula (58) satisfies (59). Other tests of 
the correctness of the complicated manipulations leading to (58) are the 
agreements with formulae previously derived for various special cases as 
discussed in succeeding sections. 

In order to make use of the formula (58) one must have some estimate of 
the magnitudes of the combinations of electron wave-functions defined in 
(54). A convenient rough approximation, with which we shall be contented 
for the present, is to use a point nuclear charge and to retain only the first 
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term in the expansions of the radial wave-functions in powers of (pp); this 
procedure will be termed the zeroth-order approximation. Note that no 
expansion in powers of (aZ) is implied. Retention of more than the first term 
in the (pp)-expansion is probably not justified in a calculation in which the 
finite distribution of nuclear charge is ignored. Series expansions for f, and 
g,. are given by Rose et al. (1953); the phase-shift is given by Rose (1937). 
The presentation of the zeroth-order approximations obtained with the aid 
of these expansions is facilitated by the use of the following definitions: 


(60a) L,,.(fo0) = 15 *’T,T,[L,»(¢o0) cos (€¢,—€,) +Ly»(¢ol) sin (e,—e,)], 
(606) L,,(¢ol) = 49” ’T,7T,[Lyr(fo1) cos (€,—€,) —Ly»(to0) sin (€,—€,)] 
vy, = [~+1)*-a°Z"]} yy =aZEp™ 


(2pp)7-? E@rt1) | T(y,+iy) 
ai r (27,+1 ) | T'(yot1y) | 


= 3 ryt arg I'(y,+7y). 


OF) TRE 


T, 


Z and a@ are the atomic number and the fine-structure constant. The M- 
and .\-functions are to be expressed in forms analogous to (60) except for 
factors p~? and p~! to be included on the right side for the former and latter 
respectively. In the zeroth order the results are: 


(61a) Ly»(010) = [(u+1+y,)(v+1+y,)+a2Z?] 
L,(110) = E“[(u+1+y,) o+1+7,)-—a2Z?] 
L,.(011) = —pE-'L,,(0—10) 
L,,(111) = 0 
L,»(0—-10) = —aZ[(ut+y,)— (v+7,)] 
Ly »(1—10) oe —aZE[pt+ytet+y,+2] 
L,,(0—11) = pE~'L,,(010) 


L,,(1-11) = 0. 
(616) My»(010) = [(u+1—y) 0+ 1-7.) +022] | 

M,,(110) = E™ [(utl—yy) (v+1—y,) —0?2?] 

M,,(011) = —pE-'M,,(0—10) 

M,,(111) = 0 


M,,(0—10) = aZ[(u—v,) —(v—7,)] 
M,.(1—10) = aZE“[p—y,+7—7,4+2] 
M,,(0—11) = pE-'M,,(010) 
M,,(1—11) = 0. 


(61c) N,»(010) = —aZ[p—yytv+y+2] 
N,,(110) = aZE“[u—y,—-v—-7,] 
N,,(011) = —pE-1N,,(0—10) 
N,,(111) = 0 
N,»(0—10) = —[@t+1—vy,)(v+1+y7,)—a?Z?] 
N,.(1—10) = E“[(u+1—y,) (v+14+y7,)+a2Z?] 
N,.(0—11) = pE-N,,(010) 


~, r 


N,,.(l—11) = 0. 
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The combinations of electron wave-functions of the forms (54) most 
commonly needed are those for which yu = v. To simplify the notation a 
little for use in the complicated formulae of succeeding sections let us write 
L,(¢or) for L,,(fer), etc.; also L,,(010) will be abbreviated to L,, etc. In 
the zeroth-order approximation one finds 


(62) L,(110) = y,(v+1)'E“'L, 
L.(O—11) = pE“L, 
L,(1—-10) = —aZy,-'L,(110) 


L,(011) = Z,(111) = L,(0—10) = L,(1—11) = 0. 


With one exception relations analogous to (62) hold also for the W- and 
.\-functions; the exception is a negative sign in the formula relating 7,(110) 
to M,. The formulae for L,, M,...R, obtained from (60) and (61) are in 
agreement with the first terms of the expansions given by Greuling and by 
Pursey. 

The W-coefficients occurring in (58) can usually be calculated from the 
algebraic formulae listed by Biedenharn et a/. (1952). It may sometimes be 
more convenient to use tables, references to which are given by Rose (1957). 
The X-coefficient, almost identical with the A-coefficient, is now also well 
tabulated (Kennedy et al. (1954); Smith and Stevenson (1957)). However, 
the values of the particular coefficients needed here are usually easily obtained 
without the aid of the tables. If the source is unoriented or if the neutrino is 
unobserved then either & or g’ is zero and 


yf Dh 
Al je je’ Je 
k’ g' k 


reduces to a Racah coefficient multiplied by a trivial factor. Formulae for 
other reducible cases are given in the Appendix, where also are given simple 
formulae for evaluating the 49's. U(JLE&) and A,, are easily computed from 
their definitions, relations (21) and (47) respectively. The only coefficients in 
(58) not yet discussed are the vector-addition coefficients of the forms 
(l'00|L0) and (jj7’ 3 —3/J0). References to papers containing algebraic 
formulae for certain special cases of vector-addition coefficients and to tables 
of values are given by Rose (1957). Rose also quotes the well-known formula 
for (/l'00'L0), a table of which has been prepared by Sharp et al. (1953). 
De-Shalit (1953) has published a table of values of (j7’ 4 —}|J0), in which 
we have found two errors. The corrected values are 


991 1 
222 2 
4. SPECTRA FROM UNORIENTED SOURCES 
The correction-factor for the shapes of spectra emitted by unoriented sources 


with no neutrino detection is obtained by integrating the So(6, 6’, ¢—¢’) of 
(48) over all electron and all neutrino directions, noting that Z(k) = 6x0; 


991 1 
as ” 4 ee 
30) = —(12/715)* and ee 5 





40) = (81/715)?. 
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thus [{SdQd2’ = R(00000). Considerable simplification occurs in R(gg’kk’t) of 
(58) when all the parameters take the value zero; o is restricted to +1, and, 
because of (57a), only + = 0 occurs. In the discussions of this section, and in 
fact of all further sections, we shall drop the terms of degree of forbiddenness 
higher than the dominant. This is done by omitting the summations over 
Ly, Le, 1)’, and 7,’ in (58). L is assumed to be fixed by the type of transition 
and by the values of — and 7 associated with K,; the value of /’ used is L—/. 
For some choices of Z and K,, two values of J are to be ‘used. To simplify the 
notation a little we shall write (K,; JZ) for (K,; JZ; 0). The dominant terms 
in the correction-factors for all five interactions are contained in the formula 


(63). 
(63) R(00000) = >> (—1)*2(1+8)inz)” 'Dnons 


nigne2 


D> 83158 52gU (IL ake) U (IL ibs) (Kaz; JeLk2)(Kui3 JsL1)* 
J 


| yv 4 j yp 4 fj 
> 14a Ler 4 jf’) AdDim—r 34 7 1,010¢""*” 
on | Ilo & J ty: fe ZZ 
v+1 b J v+1 2 J 
+Aol Le—v—1 3 7] Ad Li—v—1 3 7’ JM,(¢10)g"**?-”~* 
Le f& J Ly & J 
v+1 2 Jj ee ae 
+(—1)?*"] Ad La—v—1 4 Oj’ ] Ag Lime 3} fj’ 
L» f& J i &. & 
v3 fj mi 64 J 
+(—1)'Aol Le—» § 7’ |] Ad Li—v—1 § 7’ J IN(G10)g! tT 
Le & J Ly; a FF 
S = mbitnbothitveotmtnmtiiténc 
6 = $it6e j = vt}. 


Formulae for the shapes of spectra in transitions of arbitrary forbiddenness 
have already been given, by Greuling (1942) for pure interactions and by 
Pursey (1951) for mixtures. In these works the nuclear matrix elements are 
defined in terms of certain Cartesian tensors which are not readily related, 
in general, to the spherical tensors used here. The formulae derived from 
(63), with the aid of values of 4 as supplied in the Appendix, are all identical 
with the earlier work, if certain sums of contractions of the Cartesian tensors 
are assumed to be equal to products of pairs of reduced matrix elements. 
These equalities may be grouped into six types, one of which, (64), is quoted 
as an example. 


(64) O1(K,y xr) -Q0,* (Ke Xr) = (Ki; EE (Ko; Ly. 


The symbol Q;(K Xr), in the notation of Pursey, stands for the matrix 
element of a tensor of rank Z constructed from K X r andr, and is equivalent 





I ES LL I IN TT LN IE LT IIT LTE TITEL OT A AE LENO TO OIL TT LT a RE, AT A RET TY EIRENE A EPR RET RET AS LE 





Seg ty seeegen tne Sin ert ne SER ST ERT Te NRE REE EDS ee atte nee jo ere taes 


OES ER Rae ries get eres oe Ere 


SREP GO VASE SRE EG LORE ME ABS OR A RT LO LITT SIRO ATENEO tt CERT OE 
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to Greuling’s notation (Z!)~! Q,;(K Xr,r). The dot between the two factors 
on the left side of (64) implies contraction of the tensors plus a summation 
over M,. Each of the other five conditions equates a left side similar to that 
in (64) to a product of reduced matrix elements, except for the possible 
presence of a factor of modulus unity. The six relations between the matrix 
elements of Greuling and Pursey and the reduced matrix elements of spherical 
tensors used in this paper can all be obtained from the following set of trans- 
formations: 


(65a) Qz(K) > (K; L, L—1), 
(655) i Q,(Kr) > (K; LL), 
(65c) Q.(K Xr) > (K;LL). 


Though a direct proof of the six equations implied by (65) has not been found, 
they can easily be shown to be true when Q,(a) is a scalar or a vector. Only 
one type of scalar occurs, namely Qo(Ar) = f W*K @dr. From the definitions 
of Section 2 it follows that 
Qo(Kr) = (K; 00) 65,4,; 
thus 
Qo(Kir)-Qo* (Kor) = (Ki; 00) (Ko; 00)*, 

the same result that is obtained by the use of (658). In the formulae of Greuling 
and Pursey there are three different kinds of vectors, namely Q;(a) with a 


equal to K,, K,r, or K, X r. One can show that in each case a4 can be written 
in the form g7,#(1, Z), g being a unimodular factor; Z is 0 for K,, 1 for K,r 


and K, X r. Hence 
3 2 a * 
2 ( J v'a,0dr)( | ¥*6 4dr) , 
Ms fal 


2, ([vretae ) (fvorecr)" 


= Lake Chas LE) (Ky; RL,)*. 


(66) Q:(a)-Q:*(b) 


Now because g = 1 when a is K or K Xr and g = —7 when a is Ar, the 
six relations contained in (66) are identical with those generated from (65). 
In deriving the value of g fora = K Xr we have used the fact that 


(67) A® = —i2' (11 yy—p|1u)B'C””, ifA=BxC. 


As is well known, the works of Greuling and Pursey omit some first-forbidden 
terms; these terms, containing matrix elements of tensors of rank zero, are 
among those discussed in Section 2 under transitions of the U class with 
J; = J;. The terms omitted by Greuling had already been given by Konopinski 
and Uhlenbeck (1941); some of the missing mixture terms have been supplied 
by Smith (1951). The cross-terms for P mixed with V, 7, or A given by 
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Pursey, though evaluated correctly, are never of dominant order of forbidden- 
ness if either 7 or A is present in the 8-interaction; Greuling’s term containing 
Q,-1([@ Xr],r) in his C,y is in the same category. Our results for all first- 
forbidden transitions are exhibited explicitly in R(00000) of Section 9. There 
is complete agreement with Konopinski and Uhlenbeck; to get agreement with 
Smith it is necessary to assume that he has replaced 8 by 1 throughout. In 
papers on allowed or first-forbidden decay the notation on the left of (68) has 
usually been used for the scalar and vector Q, already discussed. 


(68a) f A = Qo(A), 
(685) JB = Q,(B). 


For the first-forbidden transitions with J; equal to J, there is another type 
of scalar of importance. In (69) the notation is defined aid the matrix element 
related to the reduced matrix element used in this paper; K is any one of the 
four vector Dirac operators. 


(69). fa = JwK- rear = i(K;01). 


The second-rank tensor occurring in the first-forbidden correction-factors of 
Konopinski and Uhlenbeck is related to our notation by (70). 


(70) td |Bisl’ = |Q2(6)|* = |(6; 21)|* 

tJ 
The 7A-mixture terms with J = 0 included in R(00000) of Section 9 do not 
seem to have been published before; these terms have been included for the 
sake of completeness, though there is considerable evidence for the absence 
of 7A mixing in the correction-factor. 

The correction-factor (63) contains the results derived for special cases 
by Spiers and Blin-Stoyle (1952a). If an apparent misprint is corrected by 
including a factor 2” within the square root in their definition (18) of N(Z), 
then the expressions i to iv of Banerjee and Saha (1954) seem to agree with 
(63). 

Because most of the properties of (63) are already well known, discussion 
will be limited to minor features. By inspection one can see that the co- 
efficient of a matrix-element product for which Z; = Ly and & = & contains 
no .V,(¢10) term if ¢ = 1, i.e. if the two matrix elements are drawn one from 
STP and the other from V4. Less obvious is the absence of the V,(¢10) term 
for ¢ = 0 when L,; = Lz and &, ¥ &s. In the proof of this statement use is 
made of the identities 


v+1 3 vt+3 v+1 4 v+} 
(71a) Ad L—v—1 § L—v—}] = —[L71(L4+1)]' Ad L—v—1 4 L—v—-3 J, 
E  | = a a 
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v v+2 vy & wth 
(716) Ag L—v 4 L—v—}] = [L“'(L+1)}' Aq L—v 4 L—v—}], 
L L L@é@ 


the justification of which depends upon formulae given in the Appendix. Note 
also that, because of the “‘stretched"’ orbital angular momenta, /+/' = L, 
at most two values of j’ occur for each value of ». 

It is perhaps not inappropriate to give some defence of the introduction of 
yet another notation for nuclear matrix elements. The use of (K;JL) as 
defined by (23) and (20) yields correction-factors in precisely the customary 
form, and also gives a concise description, through the values of J and Z and 
through the evident scalar or vector character of K, of the total, orbital, and 
spin angular momentum carried off by the emitted particles. From the symbol 
K one can tell immediately which 8-interaction is involved. No cumbersome 
complex quantities occur in the correction-factor. An advantage of the use 
of spherical tensors is that the values of the nuclear matrix-element ratios 
found by fitting the theoretical to the experimental correction-factor are more 
easily compared with the shell model. 


5. POLARIZATION OF ELECTRONS FROM AN UNORIENTED SOURCE 
Let us consider the emission of electrons in the direction of a unit-vector 
n from an unoriented source, the neutrinos being unobserved. With the non- 
conservation of parity the longitudinal polarization of these electrons need 
not be zero. Let V+ and N~- be the numbers of electrons observed with spins 
parallel and antiparallel to the direction of motion. Then the polarization, 
P, is defined by 


(72) P = (N+—N-)(Nt+4N-)-1. 
The vector polarization { has been defined in such a way that 
(73) P=C€-n. 


Because only one direction has been defined in the experiment under dis- 
cussion, { must be parallel or antiparallel to n; i.e. |{| = |P. 
Using the results of Section 3 we have 


(74) pf So(6, 0’, ¢—¢')dn’ = > (-1)'nt f Si (8, 0’, — $')d’. 


It is easily demonstrated that 
75) n* = (4/3)'V(8, $). 
Substituting (75) in (74) and using (48) we find that 
(76) P = R(10001)/R(00000). 
We shall call R(10001) the modified correction-factor; the polarization is 


then the ratio of the modified to the normal correction-factor. A formula for 
R(10001) similar to (63) for R(00000) may be obtained from (58). 
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(77) R(10001) = >) (—1)8*2(1+-8:an:)— > D> (Sr1Gnani+5r0l'nens) 


nigne 


a 6 5.96.7, U (J L2é2) U(JLié1) Rus: JoL2) (Kas; JiL1)* 


| . 2 J v ; fj : 
dy Ad Lo-» 4 JAA Li—v 4 fj’ JLA(¢—1r)g’ 
v,j | yi és J z. fy yj 
tl} j ahd. dea 
+Aq Lo—-v—-1 § 7 JAh Li—-v—-1 3 7’ JMG te)g?tt i * 
Le £2 Ji Li & ae 
vtl 3 y £ oS 
+(—1)**! Ag Le—v—1 5 is Ao Li-v i - 
LL» £2 a ye & J 
vi 3 j vl 64 Oj a 
+(—1)""4y Lo—v 3 £ Ag Ly—v—1 3 J’ Nite Oy. 
Le & J Li J | 








Comparison of (77) and (63) shows that it is possible to formulate a simple 
set of rules for generating the modified correction-factor from the normal 
correction-factor. Because of (576) L,(¢—17r) and M,(¢—I17r) occur only if 
¢+7 is odd. The terms in the modified correction-factor containing L,(¢— Ir) 
may be obtained from the normal correction-factor by the following substi- 
tutions: 

(i) If ¢ = 0, Das, > (—1)"Gaen,, and L,(010) > L,(0—11). 
8) Gi) IEF = 1, Dens > (—1)" Ferny, and ZL, (110) > Z,(1—10). 


Rules precisely analogous to (78) hold for M,({—1r). The same is true for 
the terms containing V,(¢—17) for which ¢+7 is odd. However, because of 
the absence of a restriction of the form (576) on the \-functions, there are 
additional terms in the correction-factor, namely those containing V,(1—11) 
and N,(0—10). As will be seen later these additional terms are not very 
important. They do not occur in allowed or U class transitions, except in the 
rather special AJ = 0, first-forbidden type. V,(1—11) is present only if the 
8-interaction is a mixture of contributions from STP and VA. The coefficient 
of N,(0—10) vanishes if the @-interaction is invariant under time-reversal. 
It is thus unlikely, with present ideas about the 8-interaction, that either type 
of additional term actually exists. Also, the zeroth-order values of V,(0—10) 
and N,(1—11) are zero; the first nonvanishing terms in the (pp)-expansion 
for a point charge are: 


(79a) N,(0—10) = —aZE~"y,(2y7,+1)"'p-”T,’, 


(796) N,(1—11) = pE(v+1)7,(2y->+1)—'p-”’T,?. 





Because of their probable unimportance, no mention of the additional terms 
was made in the preliminary report on this work (Lee-Whiting 1958). Neverthe- 
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less, it is possible to give rules for obtaining the additional terms from the 
normal correction factor. Thus the coefficient of 


Gain: (Kunz; J2L2)(Kny; JiLi1)*N (1-11) 

is equal to (—1)':t1+!(1+6,,,,/) times the coefficient of 
Dayan’ (Kn2; J2L2)(Kny JiL1)*N (010) ; 

the index m,’ refers to K,,, = 8K,,. Similarly the coefficient of 
Tom (Knz; J2L2)(Kni; JiLi1)*N (0-10) 

is equal to (—1)*:+*! times the coefficient of 
Dias? as: JoL2)(Kny'; J,L;)*N,(110). 


Though each of the additional terms in the modified correction-factor can 
be generated from two terms of the correction-factor, only one is to be used. 

To get a rough idea of the magnitude of the longitudinal polarization we 
turn again to the zeroth-order approximation discussed in Section 3. Then the 
troublesome additional terms of the last paragraph disappear. The remaining 
combinations of electron wave-functions can be estimated with the aid 
of (62). Note that G,,,, is always accompanied by the factor p/E = v/c, 
Pion, by —aZ/E. Thus the terms which would show a violation of time- 
reversal invariance, I’,,,,, are always accompanied by the factor aZ/E, as 
shown by Jackson et al. (1957) for the allowed transitions. Hence to look 
for the failure of time-reversal invariance one should use a large value of Z 
and a small value of E; the T,,,,, terms will not be present (in the zeroth-order 
approximation) unless the 8-interaction contains contributions from both 
STP and VA. 

If in addition to the zeroth-order approximation we assume no mixing of 
STP with VA and the validity of the two-component theory of the neutrino 
(Lee and Yang 1957), the expression for the polarization is particularly 
simple. No terms with ¢ = 1 appear, all T,,,, are zero, and Grsn, = +Daponi, 
the choice of sign being independent of m2. and m;. Every remaining combination 
of electron wave-functions is equal to p/E times the corresponding com- 
bination in the normal correction-factor. Thus the polarization reduces to 
+v/c. This result is consistent with the now well-known fact (see Curtis 
and Lewis 1957) that the magnitude of the longitudinal polarization must 
approach unity as v approaches c, regardless of the order of forbiddenness. 

One can show that the polarization is not always +v/c by working to the 
second term in the expansion in powers of (pp). Let us define the reduced 
polarization, P’, as the ratio of |P| to v/c. Then for an allowed transition in 
the field of a point charge we find that 


(80) P’ = [Lo(0—11)/L (010) ] (w/c) 
= 1+2aZpylE (vot) (2y04+) 73; 
it is here assumed that ST is not mixed with VA. Note that P’ is greater 


than unity for negatrons, less than unity for positrons. Dr. L. G. Elliott has 
pointed out that this effect is explained qualitatively by the assumption 
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that the 8-particle is created with a polarization equal to the value of v/c 
at its point of creation in the electrostatic field. The expansion of the type (80) 
in powers of (pp) is probably of little use because of the relatively large effect 
of the finite distribution of nuclear charge. To provide evidence for this 
statement, P’ for an allowed transition has been evaluated using the wave- 
functions of Sliv and Volchok (1956) for the case of the nuclear charge dis- 
tributed uniformly throughout the volume of a sphere of radius 1.2 K A'® 
fermis. The quantities tabulated are essentially our F*-4(p, 0); hence we 
have had to put p equal to zero. The behavior of the P’ so computed is very 
similar to the P’ of expression (80). In each case P’ approaches unity mono- 
tonically as E increases or as Z decreases. With the finite charge distribution, 
however, the departure from unity can be considerably greater than with 
the point charge; for example, at Z = 95 and E = 1.1, P’ is 1.021 for the 
distributed charge, 1.005 for the point charge. 

Expressions for the longitudinal polarizations of electrons emitted in 
allowed and in certain first-forbidden transitions have already been pub- 
lished. Our formula (77) gives, in the zeroth-order approximation, results in 
complete agreement with those of Jackson et al. (1957) for allowed transitions. 
Curtis and Lewis (1957) have calculated the polarization for a (AJ = 1, yes)- 
transition with the assumption that the §-interaction is ST, but without 
assuming time-reversal invariance for nuclear forces. When this latter assump- 
tion is made, the formula (6) of Curtis and Lewis is found to be contained in 
the first-forbidden terms of (77). It should be noted that though the L,’, 
M,’, and N,’ defined by these authors are exactly equal to our L,(0—11), 
M,(O—11),and V,(0—11), their R,’ is equal to — V,(0—10). In their approxi- 
mate expression for R,_1’, valid for aZ «K 1, a factor k~'! has been omitted. 
Formulae for all allowed and first-forbidden polarization terms occurring with 
an interaction ST have been derived by Alder et a/. (1957) with the restriction 
(aZ)?<<1. The (pp)-expansions of the wave-functions have been cut off in such 
a way that all terms in p~?, p~', and p® have been retained. The use of our 
zeroth-order approximation with (77) only enables us to check that we agree 
with Alder et a/. in the dominant power of p. 


6. EMISSIONS FROM ORIENTED SOURCES 


Four directions were considered in our discussion of Section 3, viz., those 
of electron momentum, p, of electron polarization, {, of neutrino momentum, 
q, and of nuclear orientation; the last will be indicated by the unit-vector j. 
Following the precedent of Jackson et al. (1957) we shall limit our discussion 
of the more complicated experiments to those in which at most three of the 
four vectors p, ¢, q, and j appear. In this section we consider only those 
cases in which one of the three is j. 

The distribution of electrons and neutrinos from an oriented source is 
obtained by putting ¢ = 0 in (48). 

(81) 44S0(0, 6", @— 6") = DY (—1)P*Z(k)[(2R+1) 2p’ + 1)1*R(ga’kg0) 
go 


x Do (g’ gm —m\k0)¥,"(0', $') VY, (6, 4). 


m 
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Because the complicated expression (81) is unlikely to be required for transi- 
tions other than the allowed, for which a simplified form has already been 
given by Jackson et a/. (1957), no further attention will be given to it. Con- 
siderable simplification of (81) results when either the electron or neutrino 
is unobserved. Taking the less important case first, we find that the distri- 
bution of neutrinos emitted by an oriented source is given by (82). 


(82) G(@’) = ir f Sv, 6’, ¢—¢’)d2 = > (2k+1)1Z(k)R(ORKOO) P; (cos 0’). 
k 


The definition of the Legendre functions P,”(x) used is that of Jahnke and 
Emde (1945). Though it is unlikely that it will ever be possible to measure 
G(6’), it is still of interest to point out that the combination (Grr+G,44) 
always occurs in the Gamow-Teller part of R(ORR00), instead of the com- 
bination (Gzr—Ga4,4) which is measured in electron polarization or distri- 
bution experiments. For experiments in which the neutrino is unobserved we 
transform (81) into (83); F(@) is, of course, the distribution of electrons from 
an oriented source. 


(83) F(@) 1x | So(8, 0’, o— ¢’)dQ’, 
>, (—1)*Z(k)R(RORKO) P; (cos 8). 


k 

If { is one of the three vectors to be considered, p must be another, and 
the third, in this section, is j. We shall give formulae for computing the 
components of { in the directions of the three members of a right-handed, 
orthogonal triad of unit-vectors, viz., n, t, t. The direction of n is that of the 
electron momentum vector p. The vector t lies in the plane of n and j and is 
perpendicular to n; the sense of t is fixed by the condition that t be parallel 
to j when n is perpendicular to j. The remaining unit-vector, t, is perpendi- 
cular to the plane of n and j; its sense is that of the increasing values of ¢. 
An orthogonal projection of all the vectors of interest is shown in Fig. 1. 
Note that 
(84) t=nxX-. 


The three unit-vectors chosen correspond to the directions for which polariza- 
tion measurements could be made most conveniently. It also happens that 
the expressions for the components of { in these directions are considerably 
simpler than those for an arbitrary set of orthogonal directions. 

Letting k represent any one of the three unit-vectors we can compute the 
component of { in the direction of k, ¢;, from the relation (85). 


—1 2 
(85) h%=k- {= | f sat | > (-1)'e" | S, “dQ’. 
e » 
The spherical components of the vector n have already been supplied—see 
(75). Those of + may be seen by inspection of Fig. 1 to be 
(86) pt! = —j2-tetto ta 


For t we use the formula (67), keeping the relation (84) in mind. The substi- 
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Z 


X 


Fic. 1. The relative positions of the right-handed, orthogonal triad of unit-vectors (n, %, t), 
the electron momentum vector p, the direction of nuclear orientation j, and the right-handed 
system of Cartesian axes (X, Y, Z) are illustrated. 


tution into (85) of the expression (48) for S;-* and of the appropriate formulae 
for spherical components leads to: 


(87a) F(0)t, = 2» (—1)*Z(k)B,(k) Px (cos 8), 
(876) F(6)o, = — » (—1)*Z(k)B.(k)[k(k+1)] *P;' (cos 6), 


(87c) F(@)t, = —i >> (—1)*Z(k)R(RORR1)[Rk(R+1)]-*P,' (cos 6). 


(88a) B,(k) = >.> (1k00|g0)R(g0kR1), 
(88b) B,(k) = [6(2k+1)]! >> (1k00|g0)W(11kk,1g)R(gOkR1). 
g 


It follows from (59) that R(RkORR1) is a pure imaginary, and hence that ¢, is 
real. i ortunately the summations over g in (88) can be carried out explicitly. 
The essence of the g-sum in (88a) is 


(89) >> {X,,.(¢or)(1k00|g0)’ —A,,’Xy»(¢+ 107) (1k00|g0) (1k1 —1|g0)}; 


the symbol X stands for L, M, or N. We have used the abbreviation A,,’ 
for A,,(k, g)[((1k1—1 g0)]-'; A,,’ is independent of g. The last term of the 
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last factor of (58) reduces to the form (89) because, for g = k+1 with k’ = k, 
we have 


(90) oAy, = o(—1)2t2t* AL, = (—1)***"A,, = Ay. 


The sums over g in (89) can be carried out with the aid of one of the unitarity 
relations for vector-addition coefficients, namely 
(91) > Gi’ m M—m|J M)(jj' m'’ M—m'|J M) = dnm- 

J 
The first term in (89) yields simply X,,({o7r), the second gives zero. In per- 
forming the analogous sum over g in (884) one encounters summations of the 
form (92). 


(92) [6(2k+1)]! >> W(11kk, 1g) (1k00|g0) 
X {X,»(¢or) (1k00|g0) — A,,’X,(¢+1lor) (1 k 1-1 |g 0)} 


With the aid of the elementary identity (93) one can show that the sum over 
g of the first term in (92) vanishes; the sum over the second term gives 


—X,,(tor) Ay, (R, R). 
(93) [6(2k+1)]* >> W(11kk, 1g) (1k00|g0) (1 & vy —v|g 0) 
g 


= 6|,|.1(1k1 —1|k 0). 


Thus both B,(R) and B,(k) can be calculated from formula (58) for R(g0kR1), 
if the quantities of the form 


[(1k00| 20) X,, (for) — Ay (Rk, g)Xy(E+1e7)] 


are appropriately replaced. For B,(k) the replacement is simply X,,(for); 
for B,(k) it is —A,,(k, k)X,,(¢+1lor). For g one may use either of the two 
values k+1, except, of course, when k = 0; B,(0) = R(10001). 

It is possible to formulate rules, similar to those of Section 6, for obtaining 
the coefficients in the Legendre expansion of the components of polarization 
from the coefficients in the expansion of the electron distribution, i.e. from 
R(RORRO); however, the rules for the .V-functions are so complicated that it 
is probably easier to work from (58) directly. 

Note that the formula for ¢, does not contain any of the quantities A,,(R, g) 
which appear in (58), while every term in the expressions for the transverse 
components ¢, and ¢, contains such a factor. If instead of { one calculates the 
expectation value of 6, formulae analogous to (48) and (58) hold with A,, 
replaced by £~'A,,. Thus the component in the direction of motion of the 
expectation value of ¢ is equal to £,, while the transverse components of the 
expectation value of 6 are equal to (1—v?/c?)? times the corresponding trans- 
verse components of {. Since ¢ is the expectation value of the electron spin 
calculated in the rest frame of the electron (Tolhoek 1956), it should not be 
surprising that the relation between { and ¢ just announced is a Lorentz 
transformation. 

The simplification of formula (58) resulting from the specialization to 
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R(RORRO) and R(gO0kRk1) is not sufficient to merit the rewriting of (58) in 
simplified form. Explicit formulae for allowed transitions are given in Section 
8. Certain first-forbidden transitions are treated in Section 9; it is felt that 
because so much space would be required for a full tabulation, this should be 
delayed until a better knowledge of the 8-interaction permits the leaving out 
of unnecessary terms. 

A general discussion of angular distribution of electrons from an oriented 
nucleus was given by Spiers and Blin-Stoyle (19526), though calculations 
were performed for only a few types of matrix elements. In the post parity- 
nonconservation era, Alder et al. (1957) have published a formulation of the 
problem, but have not carried out enough of the summations to make com- 
parison with our results profitable. No earlier work (except for allowed 
transitions) on the polarization of electrons from oriented sources seems to 
be available. Comparison with previous work on allowed and first-forbidden 
transitions is discussed in Sections 8 and 9. 


7. NEUTRINO CORRELATIONS 


Following the program outlined in the first paragraph of the preceding 
section we next examine experiments in which one of the three directions 
involved is that of the neutrino momentum vector. The angular distributions 
of neutrinos and of electrons emitted by an oriented source was given brief 
consideration in Section 6; see relation (81). When the source is randomly 
oriented, i.e. when Z(k) = 6,0, considerable simplification results. We shall 
refer to the function W(®) defined in (95) as the electron—neutrino angular- 
correlation function; 0 is the angle between the directions of emission of the 


electron and the neutrino. 


(94) W(0) = (4r)*So(8, 0’, 6—¢’), 
> R(gg0g0)P,(cos ®) (for an unoriented nucleus). 
9g 








The relation of W(®) to the coefficients R(gg0g0), given by formula (58), is 
thus particularly simple. Because g’ = g, parity nonconserving terms con- 
taining Grn, Or Iy.2, do not occur. However, combinations of coupling con- 
stants of the form A,,,,, which must vanish if the 6-interaction is invariant 
under time-reversal, do occur. The existence of such terms in the W(@) for 
allowed transitions seems to have been pointed out first by Morita (1953); 
see also Morita and Morita (1957). Our formulae for allowed transitions, 
given in Section 8, contain the results of Morita (1953); they are also in agree- 
ment with the zeroth-order results of Morita and Morita (1957) and of 
Jackson et a/. (1957). Formulae for certain cases of first-forbidden transitions 
have been given by Hamilton (1947), for Z = 0, and by Greuling and Meeks 
(1951). The relation of the present results to those of the latter authors is 
discussed in Section 9. 

The electron polarization will be treated almost exactly as in Section 6. 
Because the source is now randomly oriented, the Z-axis of Fig. 1 has no 
special significance. We are free, therefore, to assume that a neutrino is 
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observed travelling upward along the Z-axis; the unit-vector j is now parallel 
to q, the neutrino momentum vector. The derivation of the formulae for the 
components in the directions of n, t, and t of the vector polarization is very 
similar to that outlined in Section 6. 


(95a) W(0)t: = D> Qn(g’)P,’ (cos 8), 
: 
(95) WO). = — DY Qilg’)[e’(g’+ 1) Py‘ (cos 8), 


(95c) W(O)t, = —i DY R(g’g’Og’1)[g’(g’ +1) *P 2 (cos 8). 


(962) Qnp(g’) = D2 (1g'00|g0)R(gg’0g'1), 


g 
(965) Qu(g’) = [6(2g’+1) > SS (1g'00|g0) W(11g'e', 1g)R(ge’0g'1). 
g 
Just as in Section 6 the g-sums in (96) can be carried out. Both Q,(g’) and 
Q,(g’) can be calculated from formula (58) for R(gg’0g’1), if the quantities 
of the form 


(97) [(1g’00 g0)X,,(fo7r) — Ay, (g’, g)Xy»(F+1er)] 


are suitably replaced. For Q,(g’) the replacement is simply X,,(¢o7); for Q,(g’) 
it is —A,,(g’, g’)Xy.(¢+1lor); the value of g in R(gg’0g’1) can be either of 
g’+1. It is easily seen that Q,(0) = R(10001). 


8. FORMULAE FOR ALLOWED TRANSITIONS 


Within the limits set in the first paragraph of Section 6 formulae corres- 
ponding to all feasible experiments for allowed transitions have been given 
by Jackson et al. (1957). The electron wave-functions are treated by these 
authors in a manner equivalent to our zeroth-order approximation. In addition 
the nonrelativistic approximation for nuclear matrix elements was made, 


i.e., 8 was set equal to —1. If |My? and |Mg7z\* are used for |(J00)? and 
(610)|* respectively, this approximation gives 


(98) (800)? = | Mr, (800){I00)* = —|M,r|’. 
(99) (8610)? = Mer’, (8610)(610)* = —| Mer’. 
Using the definition (A17) of Jackson et al. one can show that 
(100) (100){610)* = —MrMer; 


analogous results follow for matrix elements containing 6. When the same 
approximations are made in our formulae, complete agreement is obtained 
with all the results of Jackson et al. 

For most purposes the more approximate treatment is probably adequate. 
However, it is possible that, when more accurate experiments have been 
performed, it will be desirable to allow for the effects of the finite distribution 
of nuclear charge. To facilitate the making of such an allowance formulae 
in which the zeroth-order approximation has not been made are quoted for 
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the more important cases. First we give the correction-factor for electrons 
from an unoriented source, R(00000) in the present notation. 













R(00000) = 
[D ss| (800)|?+-Dyy| (100)|?+ Drz| (8610)|?+D4 4) (610)|?]LZo 
—2[D sy (800) (100 )* + D4 (8610) (610 )*]Zo(110). 


(101) 






To compute the longitudinal polarization of electrons from an unoriented 
source we need R(10001). 


(102) R(10001) = 
[G s s| (800)? -—Gyy| (100 )|?+G rz| (8610 )|? — G4 4| (610)|?]LZ0(0— 11) 
+2[T sy (800) (100 )* + I 74 (8610) (610)*]Zo(1—10). 












The angular distribution from an oriented source requires R(10110). 





(103) —R(10110) = 
23(2J,4+1)?W (11S ,J,, 1F,) [Grr] (8610)|? —G4.4| (610)|?]Z0(0— 11) 
+20 74 (8610) (610 )*Zo(1 = 10) } 
+2 X 3-2{[G4y(610)(100)* —G75(8610) (800 )*]Z.(0—11) 


+ [Ty 7r(8610)(J00)* — T 54 (610) (800 )*]Zo(1 — 10) }. 










Recall that Z(1) is given in relation (50). To find the three components of 
the polarization of electrons from an oriented source we require: 






(104a) B,(1) = 
—23(27,4+1)2W (SJ; 1F,) {[Drr| (B610)|?+D44| (610)|2]Lo 
—2D 74 (8610) (610)*L (110) } 
+2 & 3-3 {[D,4 (610) (100 )*+ D7 5(8610) (800 )*]Lo 
— [Dy 7(B610)(100)*+D 54 (610) (800 )*]Lo(110) }. 









(104d) —2-3B,(1) = B,(1) with Ly and Lo(110) interchanged. 









(104c) 72> R(10111) = 
23(27, +1) W(11S Ji, 17) [Grr] (8610)|?—G4 4! (610)/2]L0(1 — 10) 

—21 74(8610)(610)*L.(0—11) } 

—2 X 3°? {[G75(8610) (B00 )* — G4 y (610) (100 )*]Zo(1 — 10) 

+[T, 5(610)(800)* — T' py (8610) (700)*]Zo(0—11) }. 













Now, since Lo(1—10) = —aZE-'Ly and Lo(0—11) = pE~'Lo, the terms in 
¢, with relatively large coefficients, at least for nuclei with small values of Z, 
are those containing I'yy, which are zero unless time-reversal invariance 







fails. 
The electron-neutrino angular correlation depends upon R(11010). 










(105) R(11010) = 
[Dyy| (100)/?—D 5 5|(800)|?—4D 4 4| (610) ?+4D 77! (8610) ?]Lo(0—11) 


+2[Ays(100) (800 )* — $A, 7 (610) (8610 )*]Zo(1 — 10). 






LEE-WHITING: FORBIDDEN £-DECAY 1235 


Note that the dominant terms in (105) contain Zo(0—11), the same com- 
bination of electron wave-functions that is important in the longitudinal 
polarization of electrons from an oriented source. The formula of Morita 
(1953) after correction according to Morita and Morita (1957) agrees with 
(105). As suggested by Morita (1953), the different dependences of Zo(0—11) 
and Lo(1—10) upon electron energies makes it possible in principle to test 
for failure of time-reversal invariance by measuring allowed electron—neutrino 
angular correlations. 


9. FORMULAE FOR FIRST-FORBIDDEN TRANSITIONS 


Expression (106) contains the complete first-forbidden correction-factor; 
it is divided into paragraphs corresponding to J = 0, 1, and 2. The cross- 
terms between matrix elements belonging to 7 and A with J = 0 are believed 
to be new. The correction-factor is given in full partly because the original 
publications are scattered, but mainly to facilitate the description of the 
calculation of the polarization. As shown in Section 5 the longitudinal polariza- 
tion is the ratio of R(10001) to R(00000). Each term in (106) produces, under 
the application of the rules (78) of Section 5, a term of R(10001). Because 
these rules are very simple the resulting terms of R(10001) will not be given. 
The rules for the additional terms of R(10001) are sufficiently complicated 
that we deem it useful to list these terms; they are supplied in (107). 


(106) 
R(00000) = [D4a|(ys00)|\2+ Dpp! (8-7s00)|?]Lo+2D 4p (ys00) (Bys00 )*Lo(110) 
+ D4 | (601)|*[Mo—3qNot+$q°Lo] +D rr| (8601 )|?[ Mo+ FqNot+$q°Lo] 
—2D 4 7(8601) (601 )*[ 19 (110) +$92L0(110)] —2D 4 4 (vs00) (601 )* 
[No—3qLo] 
—2D 4 p(Bys00) (601 )*[No(110) — $gZ0(110)] —2D 4 r(ys00) (8601 )* 
[No(110)+3qLo(110)] 
—2Dpr(Bys00) (8601 )*[No+ 3qLo] 


+Dss| (81 l )|?{ Mo+3qNot 4@?Lo+2Li1]+Drvy' (71 l )|*[Mo— 3qNot+ 3q°Lo 
+21] 
—2D sy(I11)(811)*[.Mo(110) +39?L0(110) +22,(110)] —2D4 7(611) 
(8611 )*[Mo(110) + §g?Z0(110) +321 (110)] 
+Drr| (861 1) 2[Mo— qNot tQ? Lot 3Lilt+Da Al (61 1 )|?{ Mo+ FqNo 
+$§q*Lot3Li] 
+[Drrz|(8al10)\?+Dyy (a10)/?]Lo+2D ry(Bal0)(a10)*L5(110) 
+2[D75(8611)(811)*+D, y(611)(J11)*][Mo—Li] 
—2D 4 5(611)(811 )*[Mo(110) — 3gN (110) —Z1(110)] —2D py (8611) {711 )* 
[Mo(110) + 3qNo(110) —Z1(110)] 
+2[Dyy(I11)(a10)*+D 77(8611)(Ba10)*][No—4qLo] 
+2[D 57(B11)(8al10)*+D4 (611) {a10)*][No+3qLo] 
+2Dy7[(111)(8a10)*+ (8611) (a10)*}[No(110) — $gLo(110)] 
+2[D sy(B11)(a10)*+ D4 7(8a10) {611 )*][No(110) + 4¢LZ0(110)] 


+[D rr| (8621 )|?-+ D4 4| (621)|2)[4q?L04+3L1] 
—2D 4 7(621) (8621 )*[4q2L0(110) +3Z1(110)]. 
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(107) R(10001) =... 
+2No(1—11)[$9G 14 (8601 ) (601 )* +G 7, (3601 ) (500 )* 

—Gp4(601)(Bys00)* 
— 3qG ra(8611)(611)*+3qG sy (811 )(J11)* 
+Gry(Bal0)(J11)*+Gr,4(Bal0) (611 )* 
—Gys(al0)(B11)* —Gyr(al0)(B611 )*] 
+2N,(0—10)[T'p7(Bys00) (8601 )* 
—$qV ya (I11)(611)* —3gT 57(B11) (8611 )* 
+ Pys(al0)(611)*— P'75(Bal0)(811 )*]. 


The angular distribution of electrons from an oriented source is calculable 
in terms of the coefficients R(RORRO); in first-forbidden transitions k takes the 
values 0, 1, 2, and 3. Though a complete list of formulae could be generated 
from (58) with little difficulty, the space required to print it would be con- 
siderable. Also, since it is certainly possible that all five types of 8-interaction 
do not occur, many terms in the formulae would probably never be used. 
Current ideas about the @-interaction suggest S7P (see Konopinski in 
Siegbahn 1955), or, more probably, VA (Goldhaber et a/. 1958; Herrmanns- 
feldt et al. 1957). In the remainder of this section results will be presented 
for STP and for VA, the cross-terms between the two groups being ignored. 
Furthermore, if an experiment on a first-forbidden decay is to be interpreted 
with reasonable ease, a transition should be chosen for which cross-terms 
containing matrix elements with different values of J do not occur. Thus in 
considering the emissions from oriented nuclei we exclude (AJ = 0)-trans- 
itions and also those (AJ = 1)-transitions in which the (J = 2)-terms make 
an appreciable contribution. A convenient abbreviation is defined in (108). 


(108) Wi(J) = (—1)¥*-24(2F, 41)! WSIS J, RS). 


(109) R(10110) =... 
+2!w,.(1) { —Gss/(811)|2[ Mo(O—11) +3qNo(0—11) +L,(0—11) 
+2.NVo:(011)+3gLoi(011)] 
—Grr|(8611)|2[My(O—11) —3qgNo(O—11) + Feq2L0(O— 11) +3L,(0—11) 
— Noi (011) +43¢L01(011)] 
—2G 7 s(811) (8611 )*[4)(0—11) —§g?L0(0— 11) —$£,(0—11) — $gZoi (011) 
+3No1(011)] 
—2T rs (B11)(8611 )*[2¢.No(0—10) — gL: (010) +3.N91(010) ] 
—Grr| (Bal0)/?L9(0—11) —2G 77(8611)(Bal10)*[Vo(0—11) — 3gLo(O—11) 
—}L(011)] 
—2G75(811)(Ba10)*[No(O—11) + 3qL0(0—11) +20: (011)]+20 rs 
(811)(Ba10)*[Vo(0—10) — Lo, (010) | 
+Gyy (711)/*{A7,(0—11) — §¥qgNo(O— 11) + 2,(0—11)+2.N 9: (011) 
— 3gL(011)] 
+G44) (611) ?[Mo(O— 11) +3¢No(0—11) + feq?L.(0—11)+4L1(0—11) 
— Noi (011) —4gLoi1 (011) ] 
+2G 41 (611)(711)*[ 17) (0—11) —bg?L0(0—11) — 32, (0—11) 
+4gLo1(011) +3.N01(011)] 
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—2T 4y(611)(I11)*[2gNo(0—10) — gZ01(010) —$N01(010) ] 
+Gyy!|(a10)|?Z5(0—11)+2Gyy(111)(a10)*[No(O—11) —3qL0(0—11) 
+L:(011)] 





+2G 4 (611) (a10)*[No(O—11) + $¢Z0(0—11) —$Lo1(011)] 

+21 4y(611)(a10)*[No(O—10)+3L0:(010)]+ ... } 
+2-4 x 3 X 5!w:(2)[Grr| (8621 )|?—Ga | (621)|?][$q2L0(0—11) 
4-85,(0<10)}+.. .. 















(110) R(20220) =... 
+6}w.(1) iDss| (all )?[Li+2Noi(0— 10) + §qLu(0— 10)] 
+Drr| (8611 )|?[4L1— Noi (0— 10) +4¢L01(0—10)] 
— (8611) (811)* (Dr s[Li— Noi(O0—10) + ¢Loi1(0—10)] 
+Ars(3Noi(0—11) —$qL01(0—11)]) 
— Dr7(8611)(8a10)*Lo:(0—10) 
+2(811)(8e10)*[D s7Lo1(0—10) + A s7rLoi1(0—11)] 
+Dyy|(I11)\?[Zi+2Noi(0— 10) — 3gL01(0—10)] 
+D44\(611)\*{4L1— Noi (0—10) —$qL01(0—10)] 
— (611) (711)* (D4 y[Li— Noi (0-10) —gLo1(0— 10) ]+ Aa y[3.V01(0—11) 
+ $qL0:(0—11)]) 
— (611)(@10)*[D4 yLo1(0—-10) +A, yLoi(0—11)] 
+2Dyy(I11)(a10)*Loi(0—10)+ ...} 
—2-? x 3 X 7$w2(2)[Drr| (8621)|?-+D44|(621)/2]Zi+.... 


(111) R(30330) = —2? x 32 & 5-4w3(2)[Gr7| (8621) |? —G, 4| (621 )\?] 
Li(0—11)+.... 






















The ST parts of (109), (110), and (111) may be compared with the approxi- 
mate results of Alder et al. (1957); the approximations involve the retention 
of only the first three terms in an expansion in powers of p and the neglect of 
(aZ)? in comparison with unity. With the use of our zeroth-order approximation 
it is possible to show that, with one exception, the p-dominant term in the 
coefficient of each power of g in our results agrees with the analogous term 
of Alder et al. The exception is in the coefficient of (38611)(311)*Ars in 
R (20220); instead of the factor 3/4 in the expression for S: at the top of the 
continuation of Alder’s Table IV on page 733, we find 3/8. 














(112). 4<ROUOIMYE) =... 
+2wi(1){Gss\ (B11)/?[Mo(1— 10) + 3¢No(1 — 10) +22, (1 — 10) — Voi (110) 
— $qLo1(110)] 
+G rr (8611) ?[Mo(1—10) —3gNo(1—10)+7eq7Lo0(1 — 10) +3L1(1 —10) 
+4.N (110) — §gLo1(110)] 
+ 2G 7 5(B611)(811)*[ My (1 —10) — 3g?Lo(1 — 10) —LZ1(1 — 10) — }. Noi (110) 
+ tgLoi(110)] 








—Tp5(8611)(811)*[$g.NVo(1 —11) —$ Noi (111) +3 qLn(111)] 
+Grri(Bal0)|2Lo(1—10) 

+Grr(8611)(8a10)*(2.NV'9(1 — 10) — ZgZo(1 — 10) +43L01(110)] 
+G s5r(811)(8a10)*[2.Vo(1 — 10) + 2gL0(1 — 10) — Lo, (110)] 
— TP g7(811)(8a10)*[2No(1—11) + Loi (111)] 
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—Gyy!| (111)|2[Mo(1—10) — 3gNo(1— 10) +2Z1(1—10) — Noi (110) + 
+ 39L01(110)] 
—Gaa\(6ll) *|Mo(1—10)+3qNo(1—10) +7sq?L0(1—10)+4Zi(1 — 10) 
+3No1(110)+§gLo1(110)] 
—2G 4y(611)(J11)*[Mo(1—10) —4g?Z0(1 —10) —Li1(1—10) — No: (110) 
—}$qLo1(110)] 
—T4y(611)(I11)*[$gNo(1 —11) +$No1(111) +43¢Z01(111)] 
—Gyy!|(al0)|?Zo(1—10) 
—G4y(611)(a10)*(2No(1— 10) + 3gLo0(1 —10) +3L£01(110)] 
+TI,4y(611)(@10)*(2No(1—11) —$£01(111)] 
—Gyy(I11)(@10)*[2No(1 —10) — $gZo(1 — 10) —L1(110)]}+.... 


Of the components of polarization of electrons emitted by an oriented 
source, only ¢, is likely to be of interest. One might hope to find that in the 
(AJ = 1, yes)-transitions ¢, could be measured in order to estimate the 
magnitude of I's7 or I'y4, and hence to study the possible failure of time- 
reversal invariance. In (112) we give the (J = 1)-terms in 7R(10111) for ST 
and for VA. From (112) one can calculate, using (106) and (110) also, the 
value of ¢, for electrons emitted perpendicularly to the direction of nuclear 
orientation. The coefficients of !'7s5 and I'4y are, however, never appreciably 
greater than those of Grs and G,y respectively. Hence a clear-cut experiment 
on time-reversal invariance is not feasible in general. For those beta-emitters 
for which cancellation of matrix elements seems to occur, e.g. RaE, it is 
possible that a measurement of ¢, would be informative. 


(113) (11010) = —|(8601)|2D pp[Mo(0—11) +2qNo(0—11) +492L0(0—11)! 
— | (By;00)!2D ppLo(0—11) +2(8601 ) (By,00)* { Drp[No(0—11) 
+ 4qL0(0—11)]+ArpNo(0—10) } 
+ | (601)|?D44{[Mo(0—11) —3gNo(0—11) +4$q2L0(0—11)]+| (ys00)|? 
Daalo(O—11) 


—2D 444601) (y,00)*[No(O—11) —4qZ.(0—11)]+... 
+4D 55) (811) ?{Mo(O—11) —2gNo(0—11) —g?Zo(0—11) —2L,(0—11) 
+8No:(011)+4gLoi(011)] 
+4Dyr\ (8611) *[My(O—11) —2gNo(O—11) +492L0(0—11) —$L,(0—11) 
—4No:(011) + 9Lo1(011)] 
+3 (8611 )(811)* {Dy s{Mo(O—11) —2gNo(0—11) + L,(0—11) +2N 01 (011) 
—2gLo1(011)]+6475No1(010) } 
+4Drr\ (Bal0)|?Lo(0—11) + 3D r7(A611)(Ba10)*[No(0—11) —gLo(0—11) 
—2120,(011)] 
2 (811)(Bal0)* { D g7[No(O—11) —gLo(O—11) +4L01(011)] 
+A gr[No(0—10) —4Lo(010)]} 
—4Dyy| (111))?[Mo(O— 11) +2qNo(0—11) —q?L0(0—11) —2LZ,(0—11) 
+8No,(01 1) —4qLo,(01 1 )] 
—4D 44) (611)|?{Mo(O—11) +2qNo(0—11) +492L0(0—11) —$L,(0—11) 
—4No:(011) —gLo1 (011)] 
— (611) (J11)* {Da y[Mo(O—11) +2qNo(0—11) 4+L,(0—11)+2N01(011) 
+2qL0(01 ] )]4+6A,4 vNoi (010) } 
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—4Dyy|(a10)|2Zo(0—11) —2D yy (111) (a10)*[No(0—11) +gZ0(0—11) 
+4Lo:(011)] 
—3(611)(@10)* { D4 y[No(O—11) + qLo(0—11) — 2Z5:(011)] + A4v[Vo(0—10) 
+2Z,),(010)]}+ acne 
+D rz| (8621 )|?[ P5q?L)(0— 11) +2¢Zo1 (011) +$L1(0—11)] 
— Daa) (621)|2[P5q2L0(0 — 11) —2gZ>, (011) +#ZL,(0—11)]+... 


—4A75(8611)(811)*gZo1(0—11) 
—4[4Dyy|(Z11)|?—Dy 4! (611)/?]gZ01(0— 10) 
+$A,y(611)(Z11)*gZo,(0—11)+... 
—3[Dr7| (8621)|?—D 4 4| (621)|2]gZ01(0—10). 


Formulae (113) and (114) for R(11010) and R(22020) enable one to calculate 
the electron—neutrino angular correlation for STP or for VA. Note that 
R(33030) contains no first-forbidden terms and therefore that P3(cos @) 
does not occur in the angular correlation. Formulae for pure interactions, 
omitting all cross-terms, have been given by Greuling and Meeks (1951); 
again only the first three terms in the expansion in powers of p are retained 
and (aZ)? is neglected in comparison with unity. When the zeroth-order 
approximation is incorporated into (113) and (114), agreement with Greuling 
and Meeks for at least the dominant term in the coefficient of each power of 
q is obtained. 

The coefficients of the coupling-constant combinations Ayy, whose presence 
would reveal the failure of time-reversal invariance, all contain, in the zeroth- 
order approximation, the “small” factor (eZ). Beta—neutrino angular correla- 
tions are not, therefore, a very promising phenomenon for investigating time- 
reversal invariance. 

If time-reversal invariance holds, then the angular correlation for (0 — 0, 
yes)-transitions reduces in the zeroth-order approximation to a particularly 
simple form, namely 


(115) 1+ (v/c) cos; 


the negative sign holds for 7’P, the positive sign for A. Our formula (115) is 
more generally valid than the similar result of Greuling and Meeks, because 
we have considered all cross-terms and because we have not neglected terms 


in (aZ)?, 


APPENDIX 

The first purpose of the Appendix is to define the various nine-parameter 
recoupling coefficients used in this paper and to state, without proof, those 
properties which have been employed in the manipulations leading to the 
formula (58). The definition of the X-coefficient is that of Fano (1951); its 
relationship to other definitions is discussed by Smith and Stevenson (1957) 
in the introduction to their numerical tables. In writing the first part of the 
Appendix much use has been made of the report of Kennedy e¢ al. (1954). 
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The generalized LS-to-jj7 transformation coefficient 





he 5. og 
at et Fs 
a ee 
is defined by (Al). In the definition, magnetic quantum numbers have not 
been written in the vector-addition coefficients; an angular-momentum quan- 
tum number is to be accompanied by the same magnetic quantum number 
each time it occurs. 
















oes iad 

(Al) (Is|j) (l's'\7’) G7 |J) = D> ALU ss’ 7’ JUL) (ss’|S)(LS|J). 
ey Lo 

The nine parameters occurring in the A-coefficient may take any integral or 

half-integral values consistent with the existence of the vector-addition 


coefficients. 
To discuss the properties of the 4-coefficient it is convenient to introduce 


the closely related but more symmetric X-coefficient. 














- £3 . £2 
(A2) ALI’ 8’ 7?) = [(QJ4DQ/+ENQLENQSHEDEX YU 8’ 7’ 
en. ee » 8 2 


The X-coefficient can be expressed as a sum of triple products of Racah 
coefficients. 





















be s 
(A383) XUV s’ 7’ |} = DS (2k+1)WsJj’, jk) Wj’ Ss, s’k) W(ll' TS, Lk). 
a . 


This expansion can be used to compute X-coefficients when simpler methods 
fail. Reflection of the array of parameters in the main diagonal does not 
change the value of X. Interchange of any two rows or of any two columns 
multiplies the X by a factor (—1)'*’+#ti+/'+J+s+8'+8) An X-coefficient with 


a parameter equal to zero may always be evaluated by means of (A4). 


















BIS oy 
(A4) XU sf j' P= (-1)*7"[(274+1) (QL +1) W(isl’s’, jL)8 55x. 
‘= 

1 


When any two parameters in one row or column are equal to 3, the 4- 
coefficient can be calculated from reasonably simple algebraic formulae, 
(A5); in this case A is equivalent to the ordinary LS-to-jj transformation 
matrix for particles of spin $. The formulae (A5), which we derived directly 
from the definition (Al) by making use of the relative simplicity of the vector- 
addition coefficients when M, = J, has already appeared in a report by 
Kennedy and Cliff (1955), who also give tables of values in both fractional 
and decimal forms; Racah (1950) and Matsunobu and Takebe (1955) have 


published formulae for closely related coefficients. 
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The third nine-parameter coefficient used in this paper, Ao, is merely an 
A modified by multiplication by a collection of miscellaneous factors, the 
repetition of which would be very tedious in writing formulae like (58). 


ew 
(A6) Adl’ 4 7’) = [(214+1)(2’+1)(2L+1)~]}*(U'00|L0) 
iC 
i .&2 
XAll’ 3 7’ Gv (gp). 
£.3:-2 


Here, of course, / and /’' are integers. The presence of the special vector-addi- 
tion coefficient ensures that Ao vanishes if /+/’+ZL is odd. Except when exam- 
ing spectra which have been measured with the highest precision, it suffices 
to retain only the first term in the MacLaurin series for the spherical Bessel 
function contained in the definition of Gy (gp); thus 


(A7) Gy(qo) = (qo) "jv (qo) = [(2l’+1)!!* = 1 X3X 5... (24D) 


Values of A» can therefore be computed using (A7), (A5), and the well- 
known formula for the special vector-addition coefficient appearing in (A6) 
(e.g. Rose 1957, page 47). Except when small corrections of order of forbidden- 
ness higher than the dominant order are being considered, it happens, as was 
seen in Section 2, that /+/’ is equal to L. In this case a small amount of 
simplification, which we shall not give, is possible in the formulae (A5); note, 
however, that the simplified formulae printed by Banerjee and Saha (1954) 
contain misprints in 


1 4 14k 1 4 i+h 
All’ 4 I'+4] and All’ 3 I'+3 
| ae | L Lev Gi 
A result valid for /+/' = L and useful in discussing the U class transitions 
is 
- ay 
ar 47 = | 
L 1 543 
when j = /+3 and j’ = l’+4, = 0 otherwise. Also for ]+/' = L we find 
(A8) (U'00|LO) = Liid'1}~*{ (22)! (20) 1[(2L) 7744. 


Next we shall derive an identity between vector-addition coefficients and 
a special class of X-coefficients. From the definitions (Al) and (A2) we find 


(A9) (34 paltpta)(j 7’ —p —ol|k’—p—o) (tk pto —p—oalg 0) 


At 
= DS [(2t+1)(2k’+1) (2/41) (2 +1) X17’ 
iu’ l P g 


K (Aj p —pll0)(47’ o —oll' 0) (ll'00'¢ 0). 
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Multiply both sides of (A9) by (47 » —p|/0)(4j7’ ¢ —oll’ 0) and sum over p 
and o. Then the use of the relation 


(A10) (474-310) = 2 forl = j4}, 


the correctness of which is made evident by inspection of the algebraic for- 
mulae for vector-addition coefficients, permits us to transform (A9) into 
(All). 


(AIL) [(2t+1) (2k’+1) (2/-+1) (2/’+1)]!(1'00|g0) X 


= (-— ] yt? -*"'114-(— 1) ren 
x {2-4 —1) *7" Gj"b— 4R'0) (th'00|g0) 
+5n(Ji'h4|k'1) (1 k’ 1 —1|g 0)}. 


Now from a recurrence relation for vector-addition coefficients given by 
Condon and Shortley (1935, page 74) it follows that 


(A12)  (97’33|R'1) = 
[R’(k’ +1) 3 +H) + (— 1)" (G+ DIGI’ 4 — 43) R’ 0). 


Using (A12) one can rewrite (A11) in the form 


- £8 
(A13)  [(2/+1) (2/’+1) }'(W/’00|g0) a(i ) = 
eis 


(—1)'2°*?[(2t4+1) (2k +1) {1+ (—1) FPG 74 — 3k’ 0) 
[(—1)7-**" (tk'00|g0) —8 Ary, 9-4 (R', g)]- 

(Ald) Ay4.5-4(R', 8) = 
2k’ (k’+1)]*(1k'1—1|g0)[(7’ +3) +(—1) 72 4-9). 


Of course Ay_4,;-4(0, g) = 0. The equation (46) used in Section 3 is equiva- 
lent to (A138). By putting ¢ = 0 in (A13) one obtains 


(AI5)  [(27+1) (27’+1)]}*(W/’00|g0) W (ijl’j’, 3g) = 
a[1+(—1)**'*"](57’ § — 41g 0), 


of which relation (25) of Section 2 is a special case. Note also that (A15) is 


equivalent to 
(A16)  Z(Ujl’j’, 3g) = 2 "(2741 (27+1))'GI 4 —4)g 9), 


valid for even values of /+/'+g. The Z-coefficient in (A16) is that discussed 
by Biedenharn et al. (1952); relation (A16) may be used to calculate 
(jj'3—4\|g0) from the tables of Z. 
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NOTES 


THE EFFECT OF ELONGATION ON THE ELECTRICAL RESISTIVITY OF 
POLYCRYSTALLINE ZINC 


H. R. PEIFFER* 


In recent years there have been many investigations of the changes in elec- 
trical resistivity introduced upon cold-working face-centered cubic metals 
(Broom 1954). These changes have been correlated to changes in the defect 
structures of the material. Very recently experiments have been performed 
on a body-centered cubic metal such as molybdenum (Martin 1957; Peiffer 
1958). A further experiment concerning these changes in a hexagonal metal 
such as zinc would be most interesting. The research described in this note 
is concerned with changes in the electrical resistivity of zinc produced by 
tensile elongations. 

The zinc used in these experiments was 99.999% pure. The specimens were 
first annealed at 150° C for 1 hour and then were elongated at room tempera- 
ture up to approximately nine per cent. The specimens were transferred 


6 
% ELONGATION 


Fic. 1. 


*Work performed at RIAS Inc., 7212 Bellona Ave., Baltimore 12, Md. 
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immediately after elongation to a Dewar flask containing liquid nitrogen. 
The resistivity of the zinc was measured in this container by means of an 
L & N Precision Kelvin Double Bridge. The specimens were periodically 
warmed to room temperature and were allowed to remain at this temperature 
for long periods. They were then again immersed in liquid nitrogen and the 
resistivity was subsequently measured to determine whether any recovery 
had occurred. None was observed even after the zinc was held at room tem- 
perature a few days. 

The resistivity increase as a function of strain is shown in Fig. 1. The rate 
of change of the resistivity with elongation decreases as the elongation in- 
creases. Furthermore, the total resistivity increase anneals out in the region 
of 90° C to 120° C, which is around the temperature range of recrystallization 
for zinc. It appears, then, that the total increase in the resistivity is caused 
by dislocations. 

Since no recovery of the resistivity is noted at room temperature, one is 
led to believe that few or no point defects are formed. However, there is a 
likely possibility that the point defects are formed and diffuse so very rapidly 
that their contribution disappears before the specimen is immersed in liquid 
nitrogen. 


















1. ‘Broom, T. 1954. Advances in Physics, 3, 9. 
2. Martin, D.G. 1957. Acta Met. 5, 371. 
3. PEIFFER, H. R. 1958. J. Appl. Phys. In press. 
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LETTERS TO THE EDITOR 


Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not later than six weeks 
previous te the first day of the month of issue. No proof will be sent to the authors. 


Decay of Heavy Hyperfragment 


The study of heavier hyperfragments can give very useful results regarding the interaction 
of the A® particle with different nuclei. Looking through the literature one finds hardly any 
report on hyperfragments heavier than Be. The reason is that the study of such events is 
difficult because the decay mode for elements heavier than helium is predominantly non- 
mesonic, which leads to a high nuclear excitation. From such a disintegration, emission of 
more than one neutron is quite common, which does not help to analyze the event completely. 
But even an incomplete analysis of such events can help in the understanding of the inter- 
action of the A® hyperon with nuclear matter by providing information about the stimulated 
decay processes of the A° hyperon in different nuclei, about which we know very little at the 
present time. 

We are presenting here a partial analysis of an event which is interpreted as a non-mesonic 
decay of either a boron or carbon hyperfragment. This event was observed during the normal 
scanning of G-5 6-in. X4-in. stripped emulsions of 400 » thickness which were exposed near 
Guam (Marian Island 6° N. geomagnetic latitude) at an altitude of 102,000 ft. 

A slow, unstable, and heavily ionizing nuclear fragment (Track 1) is ejected from a primary 
star A, of type 7+0n, comes to rest and disintegrates, producing a secondary star B with four 
prongs which are not coplanar. (The projection drawing of the event is shown in Fig. 1.) In 
the primary star two prongs were due to doubly charged particles and the rest were due to 
singly charged particles. All these tracks were followed to the end of their range or until they 
left the stack. No evidence for the coproduction of a heavy meson or hyperon was found. 


Fic. 1. A nuclear fragment 1 from a primary star A stops in emulsion and produces a secondary star B which 
has four prongs. Both tracks 3 and 4 were produced by a proton, and track 5 by a proton or an He nucleus. 
Track 2 was produced by an alpha particle. 


Data relating to the secondary star are given in Table I. Tracks 2 and 3, when compared 
with other similar tracks in the emulsion, appear to be due to doubly and singly charged 
particles, respectively. Both of them stop in the emulsion and on account of their short length 
no direct measurement can be made for mass determination. The mean width of track 2 is 
about 0.76 uv and this agrees fairly well with 0.75 » obtained from the width of number of alpha 
particle tracks from thorium stars contained in the emulsion, which confirms that the charge 
of particle 2 is two. Track 4 is steep and leaves the stack after traveling a total distance of 550 u 
in two plates. It appears to arise from a singly charged particle. The delta-ray count on this 
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TABLE I 
MEASUREMENT ON BORON OR CARBON HYPERFRAGMENT 





Range, Energy, Momentum, 
Track — Identity M Mev Mev/c 
a a 805 
b a 202 
1 5B or §C 18.5 
2 3He or *He 15.4 4.0 (*He) 172.7 
3 p.d,t 21.9 1.3 (p) 49.89 
4 p* 98.0 439.81 
5 pa 3.5 0.35 (p) 25.6 (p) 
1.0 (a) 86.4 (a) 


*Energy was calculated from ionization density. 


track was 0.42 delta rays per 1004, which indicates that the particle may be a proton. Its 
ionization density was 3.8 times minimum so its energy calculated from ionization density is 
98 Mev. Since these tracks are not coplanar, the identity of the short track 5 cannot be inferred 
from momentum imbalance. The general appearance of this track indicates a charge of one or 
two. 

Since the secondary star has four nuclear tracks, one of which is due to charge two and 
none of which is due to a negative meson, the total charge of the second star must be at least 
five. So charge conservation requires that the primary track 1 may be due to a hyperfragment 
of an isotope of boron or carbon, depending upon whether the recoil track 5 is due to a singly 
charged or a doubly charged particle. 

In order to conserve momentum, one or more neutrons may also be emitted, so that the 
visible energy released from the secondary star is greater than 104 Mev. The energy released 
from non-mesonic decay of A® hyperon is 175 Mev while for mesonic decay it is 37 Mev. As 
the energy released in the present event (neglecting even the energy taken up by neutrons) 
is much greater than 37 Mev, we consider the present decay of the hyperfragment as non- 
mesonic; the large energy of the proton track 4 indicates that the initial stimulation of the A° 
hyperon was caused by a proton. Unfortunately the shortness of the track 4 puts limitation on 
its identification. One cannot rule out the alternative possibility of track 4 being a 7-meson. 
In that case its energy calc ulated from ionization density i is about 16 Mev. Thus an alternative 
interpretation of the event is that it may be due to mesonic decay of Li or Be hyperfragment. 
The total energy released in the latter case is about 22 Mev. The non- -coplanarity of four 
tracks indicates that one or more neutrons are still involved in the mesonic decay of Li or Be 
hyperfragment. As we do not know, in either case, the energy of the neutron, we can only put 
a lower limit on Q and so it is not possible to say anything about the binding energy of the 
particle. 


We would like to thank Drs. K. E. Davis, N. M. Duller, and C. A. Randall for sharing 
with us their plates which were exposed near Guam. We would also like to thank Professor 
Marcel Schein for our stimulating introduction ‘o this field in his laboratory during the summer 
1956. 


RECEIVED May 27, 1958. P. L. Jain 
Puysics DEPARTMENT, 

UNIVERSITY OF BUFFALO, 

BUFFALO, NEW YORK. 


The Influence of the Martensitic Transformation on 
the Electrical Resistance of Sodium* 


Although Barrett (1948, 1955, 1956) has reported the discovery of a martensitic transfor- 
mation in sodium, there has hitherto been no direct evidence of the effects of this transfor- 
mation on the transport properties of sodium. Recent experiments in this laboratory have 
shown that the martensitic transformation can influence many of the thermodynamic,f 
mechanical, ft and electrical properties of sodium and lithium. Consequently a more careful 
study of the low-temperature electrical resistance of sodium was made and some results of this 
study are reported here. 


*Issued as N.R.C. No. 4855. 
tSee Martin (1958). 
¢L.Verdini (private communication). 
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The electrical resistance of three specimens of very pure sodium was measured between 
about 100° K and 4° K as a function of both decreasing and increasing temperature. Two of 
the specimens (Nos. 1 and 3) were bare sodium wires about | meter in length and about 0.5 
mm in diameter. The third specimen (No. 2) consisted of a glass capillary tube about 40 cm 
long and about 0.2 mm diameter$ filled with sodium. The residual resistance ratios (R4.2./R273) 
of the three specimens were as follows: No. 1, 4X10~*; No. 2, 2X 1074; and No. 3, 3X10™. 

The phase of sodium stable at high temperature is body-centered cubic and, according to 
Barrett, no spontaneous change to the low-temperature hexagonal close-packed phase occurs 
until the specimen is cooled to 36° K. By 4° K, about 5% of the specimen may transform 
spontaneously to the h.c.p. phase, the exact amount depending on the purity, grain size, and 
surface condition of the specimen. On the warming up of the specimen, spontaneous recon- 
version should begin at about 45° K and be complete in the temperature region between 80° 
and 110° K, according to the state of the specimen. By temperature-cycling each specimen 
and measuring its resistance at a suitable fixed temperature (40° K) we were able to show that 
all the specimens began to transform and to revert at approximately these temperatures. We 
had no method of determining the actual fraction transformed at any temperature but by 
assuming that the measured resistance at 40° K was linearly related to the proportion of each 
phase present we could estimate the approximate relative amounts transformed. A trans- 
formation curve obtained in this way is shown for specimen No. 3 in Fig. 1. The excellent 
definition of this curve suggests that this would be a useful method of investigating other 
properties of the martensitic transition in sodium. 
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Fic. 1. The transformation curve for sodium specimen No. 3. 


The measurements showed that at a given temperature conversion to the low-temperature 
phase reduced the resistance of the specimen. The maximum reduction amounted to about 
11%; this was found in specimen 1 and occurred at about 28° K, The capillary specimen, 
No. 2, showed the least change, the maximum reduction being approximately 5%, which also 
occurred at about 28° K. The percentage difference between the resistance of specimen No. 3 
measured with dec reasing and increasing temperatures is shown as a function of the tempera- 
ture in Fig. 2. The specimen was annealed at 170° K between the two runs shown in the figure; 
the tail at the high temperature end of the diagram appeared to extend up to about 120° K. 


§We are grateful to Dr. S. B. Woods for the loan of this specimen. 
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Fic. 2. The percentage difference between the resistance of specimen No. 3 measured with decreasing and 
increasing temperatures plotted as a function of the temperature. 


We draw the following conclusions from the results. 

Firstly, the fact that the transition is mot inhibited in specimens cast in glass capillary 
tubes suggests that many of the existing measurements of the transport properties of sodium 
at low temperatures should be used with caution in any comparison with theoretical pre- 
dictions. 

Secondly, if we accept, following Barrett, that the maximum spontaneous conversion to the 
h.c.p. phase even at the lowest temperature is only about 5%, then the resistivity of this 
phase at 28° K must be at least three times smaller than that of the b.c.c. phase at the same 
temperature. Indeed, a rough estimate of the temperature dependence of the resistance of 
the h.c.p. phase in the range from 25 to 60° K can be made from these results, and it appears 
that this can be approximately represented by a Bloch-Griineisen function with a characteristic 
temperature 6g ~ 280° K as compared with 6g ~ 200° K for the b.c.c. phase. 

We hope to test this conclusion further by measuring the resistance of specimens whose 
composition is more completely changed by suitable straining. 





BARRETT, C.S. 1948. Am. Mineralogist, 33, 749. 
—— 1955. J. Inst. Metals, 84, 43. 

——- 1956. Acta Cryst. 9, 671. 

Martin, D. L. 1958. Phys. Rev. Letters, 1, 4. 
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Spontaneous Magnetization of a Thin Film Measured by Electron Diffraction 


An electron beam is deflected by a magnetic field (the Lorentz effect). This effect is observable 
in a diffraction pattern obtained from a ferromagnetic substance. In this letter, we describe how 
the spontaneous magnetization of a thin film of nickel and nickel oxide was estimated by 
means of electron diffraction. A thin crystalline film of thickness 50-500 A can give a distinct 
interference diagram; this is convenient for magnetic analysis. Magnetic analysis of a thin 
film is usually carried out by means of the Faraday effect (Daniels and Wesemeyer 1958). 

A double-exposure technique was used to measure the deflection of the incident beam by a 
magnetic field. The diffraction pattern of a non-ferromagnetic film was first photographed, 
then that of the ferromagnetic specimen was photographed on the same film. It was noticed 
that the two diffraction patterns were not concentric, on account of the Lorentz effect. 

A thin film of nickel was prepared by vacuum evaporation and then was heated in air at 
250° C for 10 hours. This film was used for the diffraction experiment. Figures 1 and 2 are 
double exposures of this film and of a gold foil. In each case the incident beam was almost 
perpendicular to the film. In Fig. 1 the incident beam grazed the boundary of the film; in 
Fig. 2 it penetrated the center of the film. We notice that in Fig. 1 the rings from the specimen 
and the gold are eccentric, while in Fig. 2 they are concentric. In both figures we can detect 
only nickel oxide (NiO, NaCl type lattice, lattice constant 4.17 A). No diffraction pattern 
from nickel is visible. Since pure nickel oxide is not ferromagnetic, the observed Lorentz 
effect means that the specimen contains amorphous nickel. 

From Fig. 1 we can estimate the spontaneous magnetization H of the film. We have an 
equation 


H = (h/elLy)AZ, 


where e = electron charge (1.6 X_10-?° e.m.u.), Z = camera length (495 mm), A = the 
wavelength of the electrons (0.0304 A), h = Planck’s constant (6.6 X 107%? erg-sec), AZ = the 
deflection of the incident beam measurable from the ring eccentricity found in Fig. 1 (0.18 
mm), and / = the effective range of the magnetic field (about 0.1 mm). Therefore, we obtain 
H = 50 gauss. 


DaNIELs, J. M. and WESEMEYER, H. 1958. Can. J. Phys. 36, 405. 
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Can, J. Phys. Vol. 36 (1958) 


Note: Figures 1 and 2 follow. 





Fic.1, Diffraction pattern of a NiO-Ni film and that of gold are superposed. Incident beam grazed the boundary 


of the film. There is recognizable the ring eccentricity as the result of magnetic effect. Wavelength, 0.0304 A. 
Camera length, 495 mm. Positive enlarged 3.3 times. 


Fic, 2. Incident beam passed through the center of the film. Magnetic effect is quite faint. 


Yamaguchi—Can, J. Physics 
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